sn 

.ssr 


AO  (10  463 
$4.00 


SRDS  Report  Numbor 
RD64.160 

FINAL  REPORT 

Centfoct  Nvinbor  FA*WA*4717 
Profoct  Numbor  206«003*R 

METEOROLOGICAL  ASPECTS 
OF  THE  SONIC  BOOM  , 


THE 


COMPANY 


’’NOTICES'' 


The  information  contained  herein  is  the  result  of  a  Federal 
Aviation  Agency  sponsored  Contract  No.  FA-WA-4717. 

When  Government  drawings,  specifications,  or  other  data  are  used 
for  any  purpose  other  than  in  connection  with  a  definitely  related 
Government  procurement  operation,  the  United  States  Government 
thereby  incurs  no  responsibility  nor  any  obligation  whatsoever;  and 
the  fact  that  the  Government  may  have  formulated,  furnished,  or  in 
any  way  supplied  the  said  drawings,  specifications,  or  other  data, 
is  not  to  be  regarded  by  implication  or  otherwise  as  in  any  manner 
licensing  the  holder  or  any  other  person  or  corporation,  or  con¬ 
veying  any  rights  or  permission  to  manufacture,  use,  or  sell  any 
patented  invention  that  may  in  any  way  be  related  thereto. 


RELEASED  TO  THE  CLEARINGHOUSE  FOR 
FEDERAL  SCIENTIFIC  AND  TECHNICAL 
INFORMATION  FOR  SALE  TO  THE  GENERAL 
PUBLIC. 

4 

* 


For  Solo  by  Clooringhouio  for  Fedorol  Sciontific  ond  Toehnicol  '.nformotion 
Springfield,  VA.  22151  •  Price  $4.00 


FINAL  REPORT 


Contract  Number  FA-WA-4717 

Project  Number  206-003-R 

SRDS  Report  Number  RD64-160 


METEOROLOGICAL  ASPECTS 
OF  THE  SONIC  BOOM 


September  1964 


Prepared  by 
Edward  J.  Kane 
Tnomas  7.  Palmer 


"This  report  has  been  prepared  by  THE  BOEING  COMPANY 
under  research  and  development  sponsorship  of  the  Office  of 
Supersonic  Transport  Development  and  the  Research  Division, 
Systems  Research  and  Development  Service,  Federal  Aviation 
Agency  by  Contract  No.  FA-WA-4717,  Arthur  Hilsenrod,  SRDS 
Project  Manager.  The  contents  of  this  report  reflect  the  views 
of  the  contractor,  who  is  responsible  for  the  accuracy  of  the 
data  presented  herein,  and  do  not  necessarily  reflect  the  official 
views  or  policies  of  the  FAA.  This  report  does  not  constitute  a 
standard  specification  or  regulation." 


THE  M  M  £  SIf  ^  COMPANY 

AIRPLANE  DIVISION 
RENTON,  WASHINGTON 


ABSTRACT 


Tnls  report  is  a  study  of  the  effect  of  changing  meteorological  con¬ 
ditions  on  the  sonic  hocm  produced  during  steady  level  fll^t.  The 
Influence  of  variations  In  atmospheric  teo^jerature,  pressure,  and  wind 
on  this  noise  are  Investigated.  Simplified  methods  are  established  for 
estimating  the  effect  of  these  variations.  Ccmblnatlons  of  meteorological 
conditions  which  can  produce  anomalous  propagation  such  as  coo^lete  cut¬ 
off,  focusing,  and  extreme  lateral  spread  are  discussed.  The  effect  of 
air  txirbulence  near  the  ground  Is  considered.  A  number  of  comparisons 
with  test  data  measured  at  Oklahoma  City  (196^+)  are  presented,  and  recom¬ 
mendations  for  additional  experimental  and  theoretical  work  are  outlined. 
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SECTION  I  INTRODUCTION 


It  has  been  recognized  for  a  number  of  years  that  the  sonic  bocoi  gen¬ 
erated  during  supersonic  flight  would  be  an  iioportant  factor  In  the  design 
and  operation  of  a  commercial  supersonic  transport  airplane.  Accordingly, 
research  Into  the  factors  which  effect  the  sonic  boom  have  been  actively 
pursued  by  Indxistry  and  govemioent  agencies.  Methods  for  predicting  the 
Influence  of  the  airplane  configuration  have  been  established,  and  are 
well  substantiated  with  both  flight  test  and  wind  tunnel  data.  Ihls  ap¬ 
proach  however.  Is  valid  only  In  a  homogeneous  atmosphere  with  constant 
meterologlcad  properties  between  the  airplane  and  the  ground. 

Until  now,  considerably  less  research  has  been  devoted  to  developing 
an  understanding  of  how  an  <-..zplane'8  shock  waves  propagate  through  non- 
uniform  atmospheric  conditions.  Some  methods  of  analysis  have  been 
based  on  acoustic  proi>agatlon  through  nonuniform  tenperaturt.  regions. 
Althou^  this  Is  acceptable  for  predicting  shock-wave  locations  and  pat¬ 
terns  on  the  ground,  this  approach  yields  little  useful  Information  about 
the  shock  wave  strength  under,  and  to  the  sldi  of,  the  fll^t  track. 
Various  correction  factors  such  as  the  square  root  of  the  ratio  of  the 
ambient  pressure  at  the  ground  to  that  at  the  airplane,  ^ Pg/ Pa,  have  been 
used  to  account  for  the  effect  of  the  variation  In  atmospheric  properties 
between  the  airplane  and  the  ground.  These  aom  based  on  approximate 
analysis  of  acoustic  waves. 

A  more  detailed  approach  to  the  problms  was  taken  In  Ref.  6.  In  this 
work  the  shock  waves  were  assumed  to  propagate  at  velocities  dictated  by 
their  strength,  and  the  effect  of  pressure,  temperature,  and  wind  shear 
along  the  path  of  propagation  was  taken  Into  account.  An  additional 
value  of  this  approach  Is  that  It  allows  solutions  of  the  shock  wave 
strength  In  regions  of  focusing  where  the  acoustic  theory  predicts  totally 
unreal  values.  This  work  Is  expanded  In  Appendix  II.  It  was  also  pro¬ 
grammed  for  a  digital  ccnputer  (Refs.  6  and  83).  The  lacter  computer 
program  Is  available  from  NASA. 

Ihfr  theory  and  method  of  Refs.  6  througd^  8  have  been  used  In  this 
report  to  predict  variations  In  overpressure  which  would  occur  with 
variations  In  the  atmospheric  properties  between  the  airplane  and  the 
ground.  The  purpose  of  this  effort  was  to  determine  If  variation  In 
atmospheric  properties  could  significantly  Influence  the  boom  as  It 
propagates  between  the  airplane  and  the  ground.  Simplified  methods  for 
estimating  the  effect  of  such  variations  are  established,  and  are  sum¬ 
marized  In  Appendix  VIII.  Compeo'lsons  with  experimental  data  are  shown 
and  recommendations  for  future  work  are  also  outlined. 

Ihe  study  was  divided  Into  Investigations  of  the  effects  of  variations 
of  meteorological  conditions  In  the  atmosphere,  and  the  effects  of  local 
turbulence.  It  has  been  found  that; 
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Veirlatlons  in  temperature,  wind,  and  press\ire  can  influence  the 
boom  strength  on  the  ground. 

Variations  in  temperature  and  wind  can  influence  the  lateral  distri¬ 
bution  of  the  boom  to  the  side  of  the  fli^t  track. 

Variatioxis  in  tenperature  and  winds  can  cause  ancmalcxis  propagation 
such  as  conplete  cut-off  (no  boom  heard  on  the  groiond),  focusing 
(local  intensification  of  boom  strength),  or  extreme  lateral  spread 
(no  cut-off  to  the  side  of  the  flight  track). 

!Ilhe  above  effects  are  significant  only  at  Mach  numbers  below  1.3. 
Local  turbulence  may  change  the  shape  of  the  pressure  wave  on  its 
way  to  the  ground. 


A  number  of  horizontally  stratified  model  atmospheres  have  been  in¬ 
vestigated.  !Ibe  results  of  the  investigation  have  indicated  that  for 
flight  at  Mach  numbers  above  I.3  the  largest  influence  of  changing 
meteorological  conditions  on  the  sonic  bocm  overpressure  is  generally  no 
more  than  abouv  1^  percent  frcm  that  generated  in  the  still  (no  wind) 
standard  atmosphere.  For  fligiit  at  Mach  numbers  between  1.0  euid  I.3  the 
meteorological  conditions  between  the  airplane  emd  the  ground  may  result 
in  more  significant  variations  in  the  overpressure.  Ube  effect  of  tem- 
perattire  variations  from  the  standard  temperature-height  curve,  and  of 
various  winds  is  shown  in  Fig.  1.  These  results  eoe  typical  of  the  type 
found  throughout  the  investigation. 


EFFECT  OF  TEMPERATURE! 


MACH  NUMKR 


CUT-OFF 
IN  STANDARD 
ATMOSPHERE 


MACH  NUMBER 

Fig.  1  EH»ct  of  Tomporatun  and  ^indt  on  Sonic  Boom  Oyorprotturo. 


The  figure  shows  that  for  Mach  numbers  less  than  1.3>  temperatures 
lower  than  standard  at  the  ground  generally  reduce  the  overpressure, 
while  higher  teiQ>eratures  generally  Increase  the  boom.  For  physically 
real  conditions,  this  variation  may  be  as  much  as  ±15  percent  at  Mach 
1*2.  For  Mach  nimbers  less  than  1.3>  headwinds  generally  Increase  the 
boon  while  tailwinds  and  slde^/lnds  decrease  It.  Winds  may  ca\ise  vari¬ 
ation  In  the  oveirressure  from  that  In  a  still  atmosphere  (no  wind)  of 
as  much  as  ±20  percent  at  Mach  1.3. 

Winds  may  also  cause  the  ovex^ressure,  at  the  lateral  cut-off,  to  be 
higher  than  that  under  the  fll^t  track  for  these  low  Mach  nvu&bers. 
However,  the  situation  In  this  case  Is  not  fully  understood  becaus'*  this 
phenomena  accompanies  cut-off  where  the  shock  front  Is  locally  normal  to 
the  ground,  nils  precludes  a  reflection  of  the  wave  from  the  ground  and 
the  nonnal  doubling  of  the  free  air  overpressure  assumed  In  most  calcula¬ 
tions  (see  Section  III.D.l).  Finally,  winds  may  cause  lateral  distri¬ 
bution  of  the  boom  over  much  wider  areas  than  noxmally  predicted.  This 
phenomena  may  occur  at  all  Mach  numbers.  However,  it  need  ca\ise  little 
concern  for  two  reasons.  First,  the  oveipressure  In  the  extended  region 
drops  off  quite  rapidly  with  distance.  Second,  It  cannot  occur  If  the 
airplane  files  at  altitudes  above  those  vdiere  maximum  winds  exist.  The 
contemplated  supersonic  flight  altitudes  for  the  supersonic  transport 
are  generally  above  these  maximum  winds. 

The  Influence  of  local  turbulence  seems  to  be  that  of  a  distorting 
mechanism  idilch  deforms  the  Inccmlng  pressure  signature  on  Its  way  to 
the  ground.  Initial  efforts  to  describe  the  deformation  process  were 
not  successful,  but  the  development  of  a  more  sophisticated  approach  Is 
continuing.  Ihls  work  Is  presented,  and  the  proposed  steps  for  the 
completion  of  the  theory  are  outlined.  Present  Indications  are  that 
Interactions  of  the  shock  waves  with  certain  turb\ilent  "eddies"  result 
in  a  scattering  of  small  portions  of  the  Incident  wave  energy  to  other 
parts  of  the  wave.  This  process  would  lead  to  rounded  signatures  at 
some  points  on  the  ground  rnd  spiked  (or  very  sharp  peaks)  signatures 
at  others. 

A  limited  amount  of  the  measured  data  from  the  Oklahoma  City  flight 
test  series  was  analyzed  statistically.  To  avoid  normalizing  the  data, 
the  analysis  considered  data  measured  using  the  F-lOl+A  airplane,  flying 
at  Mach  1.5  at  an  altitude  of  28,000  feet,  at  times  when  the  cloud  cover 
was  less  than  3/l0*  observations  were  grouped  for  times  near  O7OO, 

0900,  1100  and  I3OO  Central  Standard  Time.  Results  from  this  work  in¬ 
dicated  that  the  Important  scattering  parameters  are  the  angle  of  the 
path  of  propagation  of  the  shock  wave,  and  the  time  of  day  (as  related 
to  the  turbxilent  intensity).  The  data  indicates  that  the  upper  and 
lower  bounds  of  the  overpressure  of  the  front  shock  of  a  deformed  pres¬ 
sure  signature  are  respectively  2.0  and  0.3  times  that  for  the  imdeformed 
measured  signature.  (The  shapes  of  these  signatures  are  Indicated  in 
Fig.  42.) 
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SECTION  II  EFFECTS  OF  VARIATIONS  IN  ATMOSPHERIC  PROPERTIES 


Concidered  in  1  liin  section,  are  the  cffoctc  or  vnrious  horizontally 
stratified  atmospheric  models  on  the  sonic  bcom  produced  in  stead.y  level 
flight.  The  results  are  presented  in  two  ’categories;  i  .e. , ■'effects  of 
the  standard  ntnk'sphere,  and  effects  of  variations  from  the  standard  at¬ 
mosphere.  Consideration  of  these  effects  leads  to  simplified  prediction 
o^  u.e  influence  of  variations  in  atmospheric  properties  on  the  sonic 
uoom  and,  in  most  cases,  will  make  use  of  a  computer  program  unnecessary 
for  routine  estimates  of  the  boom  strength.  Special  cases  such  as  cut¬ 
off,  focusing,  etc.  are  considered  in  Section  III. 

(A)  STANDARD  ATMOSPHERE  -  The  U.G.  Standard  Atmosphere,  1962  (Ref. 
1),  which  forms  the  basis  for  calculations  of  the  performance  character¬ 
istics  of  any  airplane  configuration,  has  been  used  to  establish  the 
basis  of  the  meteorological  effects  on  the  sonic  boom.  This  model  is 
representative  of  the  mean  atmospheric  properties  prevalent  in  the  mid¬ 
latitudes.  Boom  strength,  distribution,  and  extent  on  the  ground  has 
been  established  for  .he  standard  model. 

(1)  Sonic  Booii  Under  the  Flight  Track  -  To  predict  the  sonic  boom 
strength  in  a  nor.iicmogeneou£  atmosphere  the  Whitham  theory  (Ref.  2)  has 
been  modified  to  account  for  propagation  of  the  shock  wave  throu^  a 
region  of  '.rarying  density.  Ihe  correction  factor  used,  which  yielded 
relatlvi.'ly  close  agreement  with  test  data,  vac  a  geometric  mean  pressure 
given  by  yPa  Pg  (Refs.  3  through  5)*  This  factor  was  used,  to  re¬ 
place  the  homogeneous  ambient  pressure  used  in  the  Whitham  eqiiation  for 
the  boom  strength  under  the  airplane.  The  Whitham  equation  may  be 
written  as  shown  below: 


h3/4 


[jT^"  F(Y.tf)dY] 


Eq  (1) 


where 

K 

P 

.n 

h 

M 


7 

F(Y.  /')dY 

•'O 


Croiind  reflectivity  factor 

Ambient  pressure  at  airplane 

Airplane  distance  above  ground 

Airplane  Mach  number 

Ratio  of  specific  heats  (l.4  for  air) 

I(Y  ,$)  -  functioi  of  airplane  geometry  and  lateral 
”  location  of  the  observer  (Ref.  2) 


To  account  for  sonic  laocm  propagation  through  the  atmo- 
Bpbtre,  a  factor  aay  he  applied  to  this  equation  vhlch  Is  a  function  of 
the  ataospherlc  properties,  naisely  ambient  priissure,  P  ,  and  tenqperature, 
T  ,  alxplane  Mach  number,  and  height  above  tlie  ground.  For  f  =  -*90*  , 
under  the  fll^t  track,  Eq.  (l)  becomes: 


AP_ 


wdtr  Bt.  tnek 


=  K(P.T.Mfl)  K.P 


a 


H3/*  (y+1)» 


21/4  y 


1/1 


=  K(P.T.M»AP^ 


Specializing  this  for  the  prt^rtles  In  the  standard 

atmosphere : 


AP 


ni.  track  ^^^'^^^^Whltham  ^^Whithaa 


Eq  (2) 


Ihe  factor  Ka'1>  ft  function  of  airplane  hel{^t  above  the 
ground  and  Mach  number  In  the  Standard  Atmosphere.  Ihe  variation  of  this 
factor  vas  calculated  by  the  method  given  in  Refs.  6-8  and  is  shown  in 
Fig.  2  for  the  ground  located  at  sea  level.  It  Is  compara^^ylth  the  usu- 
ed  correction  factor  yPg/Pa«  (This  Is  equivalent  to  yPaPg  when  multi¬ 
plied  by  the  Pa  from  Eq.  (l).) 

It  can  be  seen  that  and  )/Pg/Pa  are  very  close  at  the 
lowe.'  altitudes  which  accounts  for  the  agreement  with  the  early  flight 
test  uata. 


Location  of  the  ground  above  meem  sea  level  will  affect 
the  variation  of  IC^  with  height  above  the  ground.  This  Is  primarily 
due  to  the  change  In  ambient  pressure  at  the  ground  from  that  at  sea  level. 
Curves,  similar  to  those  In  Fig.  2,  have  been  prepared  for  the  standard 
atmosphere  with  the  ground  located  at  2000,  kOOO,  and  6000  feet  above 
mean  sea  level.  These  curves  are  presented  In  Appendix  I  and  may  be 
used  \dien  the  boom  Is  being  calculated  for  areas  which  are  above  mean 
sea  level. 


Each  time  the  atmospheric  properties  vary  from  those  In 
the  standard  atmosphere,  the  variation  of  K.  (  P,  T,  M.H),  with  Mach 
number  and  hel^t  above  the  ground  changes.  Generation  of  a  number  of 
charts  similar  to  Fig.  2  for  every  conceivable  variation  would  be  a 
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AIRPLANE  HEIGHT  ABOVE  GROUND  -  1000  FT 

Fig.  2  Atmospheric  Correction  Factor  hr  U.  S.  Standard  Atmosphere,  1962. 

hiige  task  of  little  value,  lu  that  nuiterous  possible  combinations  exist. 
The  problem  of  accounting  for  small  variations  from  the  standard  atmos¬ 
phere  is  covered  in  Section  II. B. 

(2)  Lateral  Distribution  of  Boom  Strength  -  The  lateral  distribution  of 
sonic  boom  strength,  as  given  by  Whitham  (Ref.  2)  varies  inversely  as 
the  3/^  power  of  the  distance  from  the  airplane  to  the  point  of  the  ob¬ 
server  on  the  ground.  In  many  cases  the  airplane  configuration  is  such 
that  the  lateral  variation  of  boom  strength  is  also  dependent  on  the 
variation  ofI(Yo^  ^  )  with  6  .  The  combination  of  these  effects  is  given 
below; 


AP=  AP 


rHi3/4r  i(Y„, 0  11/2 

ldJ  Lhy.-90“)J 


under  fit.  track 


T 


where 


H  ■  airplane  altitude 


D  ■  distance  from  airplane  to  observer  at  lateral  distance  Y  from 
flight  track  -  + 

0  m  -  tan"^  (H/Y) 


Rearranging  this  expression  so  that  the  ri^t  side  is  independent  of  the 
airplane  geometry; 


AE.  r -  fij.  PLVl " 

AP«.a.rflt.tmkL  I(Y„.«)  J  “ 


Eq.  (3) 


The  variation  of  boom  strength  with  lateral  distance,  as 
computed  by  the  method  of  Refs.  6-8.  is  shown  in  Pig.  3  compared  to  the 
prediction  obtained  by  using  F.q,  (3). 


This  comparison  shows  that  agreement  with  Eq.  (3)  is  quite 
close.  The  maximum  differences  are  less  than  ±10  percent  and  are  con¬ 
fined  to  the  very  low  Mach  number  points  (M^  1.2). 

The  above  results  indicate  that  in  the  standard  atmosphere 
the  variation  of  boom  strength  with  lateral  distance  may  be  very  closely 
approximated  by  Eq.  (3)  with  APuadM  m.  trmek  being  obtained  by  the 
methods  of  Section  II. A. 1.  The  lateral  distribution  of  sonic  boom 
strength  predicted  by  this  method  should  be  suitable  for  most  routine 
calculations . 


Fig.  3  Variation  of  Boom  StrangtU  with  Lataral  Distance  for  U.  S,  Standard  Atmosphere,  1962. 


(3)  Uitvral  Extent  of  Boom  Strongth  -  Ten^raiu/  v^x'latlons  In  the  at¬ 
mosphere  will  cause  the  rays,  describing  the  path  .-1  the  shock  wave  from 
the  airplane  to  the  ground,  to  be  distorted.  Under  circumstances 
certain  of  these  rays  will  not  reach  the  ground  but  will  b?  refracted 
back  Into  the  atmosphere.  It  Is  i>osslble  to  determine  mathcun^ittcally 
the  lateral  extent  of  the  boom  by  detemlnlng  the  lateral  location  uf 
the  last  ray  to  reach  the  ground.  This  has  been  done  by  specifying  bhat 
the  direction  cosine  of  the  shock  front  at  the  ground  where  lateral  cut¬ 
off  occuzs  is  equal  to  unity.  Ihe  location  of  the  lateral  cut-off  In 
the  standard  atmosphere  for  the  ground  at  mean  sea  level  Is  shown  In 
Fig.  4  for  various  airplane  Mach  numbers  and  altitudes. 


Fig.  4  Lalwal  Location  of  Sonic  Boom  Cvf-dff  In  U.  5.  StonJorJ  Atnoiphoro,  1962. 

Details  of  the  calculations  and  a  general  equation  for  the 
determination  o’f  the  lateral  location  of  cut-off  on  the  fxound  for  a 
general  atmosphere  are  given  In  Section  II.B.6. 

The  location  of  the  lateral  cut-off  point  defines  the  point 
at  which  the  shock  front  begins  to  degenerate.  Beyond  this  point  noise 
may  be  heard  as  a  low  rumble.  Data  from  the  Oklahoma  City  flight  tests, 
for  Instance,  Indicate  that  the  pressure-time  trace  recor^d  beyond  the 
lateral  cut-off  location  Is  similar  to  a  sine  wave.  There  Is  no  evidence 
of  the  sharp  pressure  rise  which  Is  produced  by  a  shock.  Thus,  beyond 
the  location  of  lateral  cut-off  the  shock  front  has  degenerated  Into 
something  similar  to  an  acoustic  noise  front.  This  Is  Illustrated  by 
comparing  the  pressure -time  records  reproduced  In  Fig. 
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BEFORE  UTERAL  CUT-OFF 


BEYOND  LATERAL  CUT-OFF 


Fig.  5  Pnttun-Tim  Traces  Near  Lofsra/  Cirt-OH. 

At  the  present  time  there  are  no  data  which  Indicate  the 
rate  at  which  the  acoustic  noise  decays  to  zero,  but  the  maximum  pres¬ 
sure  difference  above  the  local  ambient  pressure  is  quite  small  and  the 
pressure  rise  time  is  quite  largje.  This  should  cause  little  annoyance. 
Pbr  practical  purposes,  then,  the  lateral  cut-off  location  would  define 
the  boundary  of  shock  wave  noise  ^rtiich  would  be  produced  by  a  supersonic 
airplane . 


(4)  Routins  Calculation  of  Sonic  Boom  in  tho  Standard  Atmosphort  -  Once  the 
airplane  geometry  inputs,  r(Yo»fl)r  (see  Eq.  (1))  have  been  established 
for  each  altitude  and  Mach  number  of  interest,  routine  calculations  of 
sonic  boom  strength,  lateral  distribution,  and  lateral  extent  for  steady 
level  flight  in  the  U.  S.  Standard  Atmosphere,  1962  may  be  obtained  in 
the  following  manner: 

(a)  Cempute  APunter  t,^|froni  Eq.  (2)  and  Fig.  2 
for  each  Mach  number  and  altitude. 


(b)  Cempute  lateral  distribution  of  the  boom  strength  frem 
Eq.  (3)  for  each  Mach  number  and  altitude. 

(c)  Obtain  the  location  of  lateral  cut-off  from  the  ciirves 
in  Fig.  4  for  each  Mach  number  and  altitude,  and  tenninate  the  lateral 
distribution  of  bocm  strength  at  this  point. 

The  above  method  Is  represented  schematically  in  Fig.  6,  and  may  be 
used  for  M  >  1.2  with  relatively  good  accuracy. 


(B)  VARIATIONS  FROM  STANDARD  ATMOSPHERE  —  Variations  in  temperature, 
pressure,  and  wind  in  the  atmosphere  will  produce  variations  in  the  sonic 
boom  strength  received  at  the  ground.  The  effect  of  some  typical  changes 
in  atmospheric  proi)erties  frean  those  in  the  standard  atmosphere  model 
are  considered  in  this  section.  Methods  of  approximately  accounting  for 
these  influences  are  presented,  and  where  approximations  are  not  possible 
the  importance  of  the  resulting  effects  are  discussed.  In  most  cases 


Involviiig  routine  calculations  it  is  possible  to  make  estimates  of  the 
boom  strength,  lateral  distribution,  and  lateral  extent  without  vising 
complex  coorputer  calculations. 


a. 


LATERAL  DISTANCE 


Fig,  6  Routint  Caleulatlen  of  Sonic  Boom  in  U.  S.  Standard  Atmoophoro,  1962. 


(1)  Tomporaturo  -  Atmosphere  temperature  variations  between  the 
airplane  and  the  ground  will  cause  the  ray  path,  which  describes  the 
motion  of  the  shock  wave  through  the  atmosphere,  to  be  deformed.  In 
general,  increasing  temperatures  along  the  path  of  propagation  (negative 
lapse  rates)  will  cause  the  ray  to  bend  up  while  decreasing  temperatures 
(positive  lapse  rates)  will  cause  it  to  bend  down.  This  distortion  of 
the  ray  path  will  Influence  the  boom  -strength  received  on  the  ground. 

The  effect  of  temperature  variations,  on  the  boom  strength  under  the  air¬ 
plane,  have  been  studied  for  a  number  of  model  atmospheres.  These  models 
(Refs.  9-16)  characterize  the  following  meteorological  conditions; 

•  Various  tropopause  heights  and  associated  temperature  gradients. 

•  Temperature  Inversions  near  the  ground  caused  by  nocturnal  radiation, 
snow  cover,  and  coastal  stratus. 

•  Multiple  temperature  inversions  due  to  mixing  and  advection. 

•  Frontal  temperature  inversions. 

•  Combinations  of  the  above. 

They  are  summarized  in. Fig.  T»  which  shows  the  standard 
atmosphere,  for  reference. 
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Flg>  7  Summary  and  Kay  of  Atmospharic  Modals. 


The  effect  of  these  variations  on  hoom  strength  was  studied 
in  each  temperature  model.  It  was  found  that  the  temperature  effect  was 
primarily  a  function  of  airplane  Mach  number  and  the  ratio  of  the  abso¬ 
lute  temperature  at  the  airplane  to  that  at  the  ground.  This  is  illus¬ 
trated  in  the  summary  of  the  calculated  data  shown  in  Fig.  8.  To  better 
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Fig.  8  Effect  af  Trnnpwafur*  Variafions  on  Sonic  Boom  Produced  Under  fhe  Flight  Track. 


display  the  variation  in  boom  strength,  the  overpressure  produced  under 
the  airplane  in  each  atmospheric  model  has  been  divided  by  the  overpres¬ 
sure  generated  in  the  standard  atmosphere  for  the  corresponding  Mach 
number  and  altitude. 

The  data  presented  in  Fig.  8  indicate  that  significant 
changes  in  the  overpressure,  caused  by  temperature  variations  from  the 
standeu^i  atmosphere,  are  mainly  confined  to  the  Mach  number  range  between 
that  required,  for  ccMnplete  cut-off  (no  boom  hecurd  on  the  ground)  and 
about  Mach  1.2.  For  Mach  numbers  greater  than  1.2  the  effect  is  rela¬ 
tively  small.  The  influence  of  temperature  variations  for  Mach  1.2  is 
shown  in  Fig.  9»  This  figure  presents  the  data  shown  in  Fig.  8  for  Mach 


Fig,  ?  Approx/maft  Variation  of  Ovarprassvro  with  Ttmpcrafvr*  for  Maeh  1,2. 


XT  fra  ««c«yw  yo  wthtw  imiTMnur.^ 


When  the  airplane  Mach  number  is  less  than  or  approximately 
equed  to  1.2  extreme  care  must  be  used  In  estimating  the  boom  strength 
generated  under  the  airplane.  Because  all  the  meteorological  factors  In¬ 
volved  may  significantly  affect  the  boom  strength,  simple  approximations 
are  not  possible  In  most  cases.  A  more  complete  discussion  of  these 
effects  In  the  low  Mach  number  range  Is  given  In  Section  III.  However, 
for  routine  analysis  It  may  be  possible  to  approximately  estimate  the 
effect  of  temperature  variations  on  the  boom  strength  by  referring  to 
Figs.  8  and  9*  The  change  at  Mach  1.2  would  be  estimated  from  Fig.  9 
and  the  variations  for  the  rest  of  the  Mach  number  range  would  be  esti¬ 
mated  from  Fig.  8.  As  the  variations  In  boom  strength  for  Mach  l.$  do 
not  usually  exceed  ±2  percent,  the  effect  of  temperature  In  this  range 
might  be  neglected  for  routine  calculations. 

(2)  PrcMur*  -  pressure  vcurlatlons  between  the  airplane  and 
the  ground  will  Influence  the  booni  t^tenflty  received  on  the  ground. 
Investigation  of  the  equations  describing  the  propagation  of  a  shock  wave 
through  a  nonhomogeneous  horizontally  stratified  media  (Refs.  6  and  7) 
Indicate  that  the  Influence  of  small  variations  from  an  established  pres¬ 
sure-height  curve  can  be  described  by  a  relationship  between  the  ambient 
pressures  at  the  airplane  and  the  ground.  If  the  established  pressure- 
height  curve  Is  taken  as  the  one  for  the  standard  atmosphere  the  Influence 
of  variations  from  It  can  be  estimated  from  Eq.  (k).  (This  equation  Is 
developed  In  Appendix  III. ) 


'  uncter 


Pa 


Eq.  (4) 


und«r 
lit.  track 


Where 


Pg 

(Pg).«. 

Pa 


Ambient  pressure  at  ground  In  model 

Ambient  pressure  at  ground  In  standard  atmosphere 

Ambient  pressure  at  airplane  (tapeline)  altitude 

Ambient  pressure  lu  standard  atmosphere  at  airplane 
(tapeline)  altitude 


The  Influence  of  pressure  on  sonic  boom  strength  Is  indepen¬ 
dent  of  airplane  Mach  number.  This  is  because  the  path  of  propagation 
(i.e.  the  retj'"  path)  is  primarily  a  function  of  the  temperature  variation 
and  airplane  Mach  number,  and  Is  generally  Independent  of  the  pressure 
variation  In  the  atmosphere. 


This  expression  has  been  checked  by  making  a  number  of  com¬ 
parisons  with  results  obtained  using  the  method  of  Refs.  6-8.  The 
atmospheric  temperature  models  shown  In  Fig.  7  were  used  for  the  compeu:!- 
son.  First,  the  sonic  boom  was  computed  In  each  model  assuming  the  U.S. 
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Standard  Atmosphere,  1962  (Ref.  1)  pressure -height  curve.  Boon  strength 
was  then  recomputed  assuming  a  pressure-height  curve  which  was  developed 
by  using  the  hypsometric  equation  for  each  temperature  model  (Ref.  17). 
iniese  two  results  were  then  ratioed  to  obtain  the  effect  of  pressure 
variations  from  the  standard  pressure-height  curve  on  the  sonic  boom, 

AP/  \  ^  comparison  between  Eq.  (U)  and 

••"P-  I  BU  irark 


i.e., 


the  computed  results  is  shown  in  Pig.  10,  for  each  group  of  temperature 
models. 


These  data  show  that  agreement  between  Eq.  (U)  and  the  com¬ 
puted  results  is  quite  close.  This  equation  should  be  adequate  to  allow 
estimation  of  the  effect  of  normal  pressure  variations,  from  the  standard 
atmosphere,  on  the  sonic  boom  produced  under  the  flight  path. 


(3)  CembiBBd  TcmpBraturB  and  Prattur#  -  Veuriations  of  both 
temperature  and  pressure,  in  the  atmosphere,  will  affect  the  strength  of 
the  sonic  boom  received  on  the  ground.  The  influence  of  these  variations 
on  the  sonic  boom  is  a  function  of  airplane  Mach  number,  altitude,  and 
atmosphere  thermodynamic  properties.  In  Section  II.B.l  it  was  observed 
that  the  influence  of  temperature  variations  from  the  stsuidard  teiipera- 
ture-height  prof^e  was  a  function  of  airplane  Mach  number,  altitude, 
and  the  absolute  temperature  profile.  This  could  be  put  in  the  function¬ 
al  form: 


=  K(M.T.H) 

umler 
(It.  track 


Furthermore,  it  was  observed,  in  Section  II. B. 2,  that  the 
effect  of  pressure  variations  from  the  standard  pressure-height  profile 
was  primarily  a  function  of  that  profile  and  airplane  altitude.  This 
could  be  put  in  the  functional  form: 


■Cd 

press 


)  *  =  K(P.H) 

/under 

n*  . _ .. 


Hi,  track 


The  combined  effects  of  pressure  and  temperature  variations 
from  the  standard  atmosphere  can  be  obtained  by  the  product  of  the  abovs 
two  functional  quantities. 


Specifically,  if  the  first  of  these  functional  quantities  is 
noted  as Kt (temperature  correction  factor)  and  the  second  is  put  in  the 
form  given  by  Eq.  (14- ),  the  product  of  these  is: 


under 
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Fig.  ,0  Efftct  of  Prozsuro  Variations  on  Sonic  Boom  Produced  Under  the  Flight  Track. 
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Thus,  for  a  nonstandard  atmosphere  the  overpressure  under 
the  flight  track  would  be  given  by  Eq.  (5)* 

[(P^J  [(]^J  Eq.  C 

where 


Kt 


(Pg) 

Pa 


■t  Temperature  correction  factor 
B  Ambient  pressure  at  ground 

B  Ambient  pressure  at  ground  in  standard  atmosphere 

■  Ambient  pressure  at  airplane  tapeline  altitude  above 
ground 


(Pa)j^  ■  Ambient  pressure  at  airplane  tapeline  altitude  above 
ground  in  standard  atmosphere 


«  Atmospheric  correction  factor  for  stemdard  atmosphere 
(see  Section  II.A.l,  and  Appendix  l) 

APu,u..v  ^  Eq.  (l)  (See  Section  II. A. 1) 

Wnittuifn 


For  most  routine  calculations,  the  temperature  correction 
factor,  Kt,  may  be  determined  by  the  method  outlined  in  Section  II.B.l. 

In  some  sx>ecial  cases,  such  as  for  Mach  numbers  less  than  1.2,  special 
methods  may  be  required  in  order  to  make  estimates.  However,  these  situ* 
ations  would  feill  into  the  specied  class  of  problems  which  EU'e  discussed 
in  Section  III.  Thus,  use  of  the  correction  factors  outlined  in  Eq.  (5) 
eliminates  the  need  to  generate  a  number  of  curves  (Fig.  2),  to 
allow  estimation  of  the  effect  of  normal  atmospheric  variations  from 
standard  conditions. 


(4)  Wind*  -  Variation  in  wind  speed  and  direction  (i.e. 
wind  shear)  between  the  airplane  and  the  ground  will  tend  to  distort  the 
ray  path  in  much  the  same  manner  a=  variations  in  temperature.  In  gen- 
er^,  headwinds  cause  the  rays  to  bend  up  away  from  the  ground,  while 
tailwinds  cause  them  to  bend  down.  Distortion  of  the  ray  paths  will 
cause  some  variation  in  sonic  boom  strength.  This  variation  was  investi¬ 
gated  by  constructing  a  set  of  model  wind  profiles  (Refs.  18-35)  and 
calculating  the  resulting  overpressures  on  the  ground.  The  wind  models 
were  selected  to  be  characteristic  of: 

•  Gradients  in  zonal  and  merldonal  wind  components 

•  High-speed  Jet  streams  near  the  tfopopause 

•  Low-level  Jet  streams  over  the  great  plains 


Each  model  was  assumed  to  be  omnidirectional,  i.e.  no  lat¬ 
eral  shear.  The  wind  models  are  summarized  in  Fig.  11  for  each  of  the 
above  categories. 

The  effect  of  winds  on  the  overpressures  received  at  the 
ground  was  studied  by  assuming  that  the  wind  velocities  were  aligned 
parallel  and  perpendicular  to  the  airplane  flight  path.  In  this  manner, 
the  influence  of  headwinds,  tailwinds  and  sidewinds  were  studied.  The 
a/erpressures  computed  by  the  method  of  Refs.  6-8  under  the  flight  track 
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with  wind  were  divided  by  the  overpressures  ih  the  same  model  atmosphere 
with  no  wind,  to  better  indicate  the  influence  of  each  profile.  These 
results  Eire  shown  in  Fig.  12  for  the  various  assumed  wind  directions. 
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Fig,  12  Effect  of  Wind  on  Sonic  Boom  Under  the  Flight  Track. 


This  figure  illustrates  that  the  headwinds  and  tailwinds 
exert  the  most  powerful  influence  on  the  overpressure  under  the  flight 
path  for  Mach  numbers  less  than  about  1.3<  Sidewinds  seem  to  have  rela¬ 
tively  little  influence  on  the  sonic  boom  under  the  airplane  regardless 
of  Mach  number.  The  results  fo^  winds  parallel  to  the  flight  track  are 
similar  to  those  obtained  by  varying  the  temperature  profile  (Section 
II.B.l).  For  Mach  numbers  above  1.5>  the  effect  of  winds  is  relatively 
small  (of  the  order  of  ±2  percent).  In  the  Mach  number  range  between 
1.2  and  1.5  the  wind  components  peirallel  to  the  flight  path  way  have  a 
significant  effect.  In  some  cases,  especially  at  low  Mach  numbers,  winds 
may  cause  focusing  of  the  sonic  boom  under  or  to  the  side  of  the  flight 
track,  These  situations  represent  a  special  set  of  cases  which  are  con¬ 
sidered  in  detail  in  Section  III.B. 

Care  must  be  taken  in  estimating  the  influence  of  wind  on 
the  boom  under  the  flight  track  for  Mach  numbers  less  than  about  1.3, 
but  for  Mach  I.3  and  moderate  winds  the  effect  may  be  estimated  from  a 
cross  plot  of  some  of  the  data  in  Fig.'  12.  This  is  shown  in  Fig.  13 
where  the  data  at  Mach  I.3  has  been  plotted  against  the  sum  of  the  speed 
of  sound  ratio  (SLg/a.^)  and  the  relative  wind  component  (IP/a, )  parall.el 
to  the  flight  path. 


J^ig.  13  Approximate  Variation  of  Ovorprasswo  with  ^ind  Along  Flight  Path  for  Moch  1.3. 
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The  data  from  this  figure  may  he  used  to  estimate  the  in¬ 
fluence  of  headwinds  and  tailwinds  at  Mach  1.3-  Figure  12  data  may  he 
used  to  estimate  the  effect  a/er  the  remainder  of  the  Mach  number  range, 
for  most  routine  cedculations.  In  routine  cedculations,  the  effect  of 
sidewind  components  could  prohahly  he  ignored  since  the  most  extreme 
winds  cause  less  than  a  5  percent  change  in  the  overpressure  under  the 
flight  path  at  Mach  1.2,  and  less  than  ±2  percent  at  higher  Mach  numbers. 

(5)  Lateral  Distribution  of  Boom  Str«n9iii  -  It  was  shown  in  Section  II.A.2 
that  according  to  theory  the  lateral  distribution  of  the  sonic  boom  strength 
varied  inversely  as  the  3/^  power  of  the  disteuice  from  the  airplane  to  any 
point  on  the  groiuid.  Hie  expression  developed  in  that  section,  Eq.  (3), 
was  compared  to  data  computed  using  the  method  of  Refs.  6  through  8  for  each 
of  the  nonstandard  atmospheric  models  shown  in  Fig.  7»  This  comparison  is 
shown  in  Fig.  l4  for  several  Mach  numbers.  No  wind  shear  is  assumed. 

The  figure  shows  that  the  maximum  deviation  from  the  approx¬ 
imate  equation,  Eq.  (3),  is  about  8  percent  for  Mach  1.2.  The  majority 
of  the  computed  points  show  very  close  agreement,  especially  at  the  Mach 
numbers  above  1.2.  The  results  of  this  comparison  would  indicate  that 
Eq.  (3)  should  yield  a  relatively  good  approximation  for  both  standard 
and  nonstandard  atmospheric  models. 

The  presence  of  wind  shears  will  also  influence  the  latered 
distribution  of  the  boom  strength.  In  some  cases  the  approximation 
obtained  by  using  Eq.  (3)  nay  not  be  sufficient  and  more  complicated 
methods  must  be  used.  The  validity  of  the  simple  approximation  was 
checked  in  the  standard  atmosphere  with  the  mean  zonal  and  high-speed 
Jet  wind  profiles  (Fig.  11 ).  The  wind  directions  were  taken  individually 
as  tailwinds,  headwinds,  and  sidewinds  and  the  distribution  of  boom 
strength  under  the  airplane  was  computed  for  each  case.  The  computed 
results  are  compared  to  Eq.  (3)  in  Fig.  15(a)  for  the  mean  zonal  wind 
profile  and  in  Fig.  15(b)  for  the  high-speed  Jet  wind  profile. 

rf 

The  figure  shows  that  for  the  mean  zonal  wind  profile  (Fig. 
15(a)),  the  agreement  with  Eq.  (3)  is  quite  good  at  all  Mach  numbers 
above  1.2.  In  the  case  of  high-speed  Jet  profile  (Fig.  15(b))  agreement 
is  not  quite  as  good,  especially  at  low  Mach  numbers.  From  this  compari¬ 
son,  it  appears  that  estimates  computed  using  the  approximate  equation, 

Eq.  (3),  are  sufficient  in  cases  of  moderate  winds  (with  maximum 
wind  speeds  less  than  about  100  feet  per  second)  at  all  Mach  numbers, 
and  for  high  winds  at  Mach  numbers  greater  than  about  1.5.  The 
results  at  low  Mach  numbers  for  the  high-speed  Jet  profile  illustrate 
some  of  the  special  cases  discussed  in  Section  III.B.  For  instance,  at 
Mach  1.2  both  the  headwind  and  tailwind  cause  complete  cut-off  so  that 
no  boom  would  be  heard  on  the  ground.  However,  as  a  sidewind  the  high¬ 
speed  Jet  would  cause  lateral  focusing  off  the  flight  track  (i.e.  over¬ 
pressures  to  the  side  of  the  flight  track  higher  than  those  under  it). 

In  the  cases  of  high  wind  speeds  in  the  vicinity  of  the  airplane  altitude, 
extreme  care  must  be  exercised  when  estimating  the  boom  strength  and 
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Fig.  15  Efhet  of  Vfind  on  Lotorol  Distribution  of  Sonic  Boom  Strongth. 


distribution  on  the  grounds  This  situation  represents  a  set  of  special 
cases  which  are  covered  In  more  detail  In  Section  III.B.  Criteria  axe 
established  In  that  section  so  that  these  special  situations  may  be 
Isolated  and  Investigated  separately  If  necessary.  For  ordinary  meteor¬ 
ological  conditions  the  approximation  established  previously  should  be 
sufficient  In  making  routine  estimates  of  the  lateral  distribution  of 
the  boom  strength  on  the  ground. 

(6)  Lataral  Extant  nf  Bnon  Strnngth  The  lateral  extent  of  the  boom 
strength  to  the  side  of  the  flight  track  In  a  still  atmosphere  Is  deter¬ 
mined  primarily  by  the  temperatxire  distribution  between  the  airplane  and 
the  groimd,  the  flight  altitude,  and  airplane  Mach  number.  The  location 
of  the  lateral  cut-off  point  does  not  necessarily  define  the  location 
beyond  which  no  noise  wl’l  be  heard  as  noted  In  Section  II. A. 3.  It  is 
possible  to  theoretically  determine  the  lateral  extent  of  noise  by  In¬ 
vestigating  the  eqtiatlons  describing  the  path  of  the  shock  front  from 
the  airplane  to  the  groiind.  These  paths  are  commonly  called  ray  paths 
and  their  lateral  extent  Is  found  by  determining  the  location  of  the 
last  ray  to  reach  the  ground. 

Consider  a  general  temperature  profile  such  as  that  shown  In 
Fig.  l6,  where  the  temperature  Is  assumed  to  vary  linearly  between  signif¬ 
icant  levels. 


The  lateral  location  of  the  last  ray  to  reach  the  ground  for 
this  profile  Is  given  by  (this  equation  Is  developed  in  Appendix  IV): 


In  some  cases  the  temperature  between  two  levels  Is  constant, 
such  as  In  the  stratosphere  of  the  standard  atmosphere.  Assuming  this 
occurs  between  levels  z.^  and  Zn,^. ,  (l.e. 
involving  (z  ,  -z  )  becomes: 

in+ 1  m 


®m+l 


(Z  -  Z  ) 

m-fl  nm' 


Eq.  (7) 


This  Is  further  illustrated  in  Appendix  IV  where  Eq.  (6) 
is  expanded  for  the  standard  atmosphere. 

The  assumptions  Involved  in  deriving  Eq.  (6)  are  that  the 
speed  of  sound  varies  linearly  between  levels,  and  that  the  shock  front 
travels  at  speeds  nearly  equal  to  the  local  speed  of  sound.  The  influ¬ 
ence  of  these  assumptions  has  been  checked  by  comparing  results  computed 
by  the  method  of  Refs.  6-8  with  those  computed  using  the  above  equations. 
A  typical  example  of  this  ccmpetrison  is  shown  in  Fig.  17  for  the  model 
temperatvre  profile  B-5  (Fig.  j),  at  several  Mach  numbers  and  altitudes. 
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The  figiu'e  shows  that  the  effect  of  the  assuniptionB  is  not 
very  significant,  as  the  results  computed  by  both  methods  are  quite  close. 

Variations  in  the  temperature  profile  will  cause  variations 
in  the  location  of  the  lateral  cut-off.  The  extent  of  the  latereQ.  dis¬ 
tribution  of  sonic  boom  strength  was  computed  for  each  of  the  atmospheric 
models  shown  in  Fig.  7* 

The  variation  in  lateral  extent  is  shown  in  Fig.  l8  for  the 
two  models  which  produced  the  widest  deviations  from  the  values  in  the 
standard  atmosphere.  The  extent  in  the  stelndard  atmosphere  is  shown  for 
reference. 


The  figure  shows  that,  in  the  extreme  cases,  vGurlatlons 
from  the  standard  atmosphere  locations  3  to  10  miles  may  occur.  In 
genered,  temperatures  lower  than  standard  on  the  ground  will  increase 
the  lateral  extent, Ymax>  while  ground  temperatures  higher  them  standard 
will  decrease Yinax* 

Wind  shears  will  also  affect  the  location  of  latered  cut-off 
off*  The  magnitude  of  the  variation  caused  by  moderate  winds  was 
studied  in  the  standard  atmosphere  with  the  mean  zonal  wind  profile 
(Fig.  U-).  The  value  of  Ymax  was  computed  by  the  method  of  Refs.  6-8, 
for  a  headwind,  tailwind,  and  sidewind.  For  the  purposes  of  compeurison 
these  values  were  divided  by  the  value  of Ymax  with  no  wind  for  the  same 
Mach  number  and  altitude  and  are  shown  plotted  in  Fig.  19  against  air¬ 
plane  Mach  number. 

The  figure  shows  that  for  moderate  winds  the  maximum  varia¬ 
tion  occurs  at  Mach  numbers  near  1.2.  For  Mach  numbers  greater  than 
about  1.5  the  variation  is  of  the  order  of  ±5  percent.  In  general,  tail¬ 
winds  and  sidewinds  (on  the  downwind  side)  increase  the  value  of  Ymax, 
while  headwinds  and  sidewinds  (on  the  upwind  side)  decrease  Ymax* 
some  cases  strong  winds  may  substantially  Increase  the  magnitude  of  Ymax* 
Ihe  conditions  required  for  this  to  occur  are  described  more  fully  in 
Section  III.B.4.  For  purposes  of  routine  calculation  the  influence  of 
moderate  winds  on  the  location  of  the  lateral  cut-off  point  might  be 
ignored  or  estimated  from  Fig.  19. 


(7)  Routine  Colculotiona  of  Sonic  Boom  in  Gonoral  Atmotphoio  -  Ulie  method  for 
calculating  sonic  boom  distribution  on  the  ground  was  outlined  in  Section 
II. A. 4,  for  the  U.S.  Standard  Atmosphere,  19^2,  with  no  wind.  A  similar 
t  procedure  would  be  used  in  making  routine  estl^tes  for  a  general  at¬ 

mosphere  with  wind.  The  following  method  should  be  used  with  extreme 
care  for  Mach  numbers  between  1.0  and  1,3,  especially  in  cases  when' 
wind  shears  are  to  be  considered.  Criteria  for  meteorological  conditions 
which  may  cause  anomalies  in  the  overpressure  and  distribution  of  the 
sonic  beam  are  developed  in  Section  Til,  These  should  be  checked  when 
makir^  calculations  in  the  above  Mach  number  range. 
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Fig.  18  Variation  of  Lataral  Extant  With  Tamparotvra  Profila. 
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Fig.  19  Variation  af  Location  of  Latoral  Cut-Off  9lith  9lind 


Once  the  airplane  geometry  influence,  KY^  ,0),  (see  Eq.  (l)) 
has  been  determined  for  each  Mach  number  and  altitude  of  interest  the 
calculation  of  sonic  boom  strength  and  distribution  may  proceed  as  follows ; 

(a)  Calculate  A  P  nt.  track  from  Eq.  (5)  and  Fig.  2 
(or  Appendix  Figs.  I-l,  1-2,  or  1-3  for  the  ground  located  above  0  MSL) 
for  each  altitude  and  Mach  number 

•  The  variation  of  the  temperature  correction  factor,  ,  with  Mach 

number  may  be  estimated  from  Fig.  8,  with  the  value  at  Mach  1.2 
estimated  from  Fig.  9' 

•  The  variation  of  the  correction  for  wind  shears  with  Mach  number  may 

be  estimated  from  Fig.  12  with  the  value  at  Mach  I.3  taken  from  Fig.  13. 

(b)  Compute  the  lateral  distribution  of  sonic  boom  strength 
from  Eq.  (3)  for  each  altitude  and  Mach  number.  (Caution  must  be  exer¬ 
cised  when  estimating  the  lateral  distribution  for  low  Mach  numbers  and 
high  winds.  Fig.  I5.) 


(c)  Obtain  the  location  of  lateral  ait-off  from  Eq.  (6) 
and  Eq.  (7)  for  each  Mach  number  and  altitude,  and  terminate  the  lateral 
distribution  at  that  point.  (Effect  of  moderate  winds  on  the  lateral 
cut-off  location  may  be  estimated  from  Fig.  19.) 

The  above  procedure  was  generally  Illustrated  in  Fig.  6, 
which  may  be  referred  to  as  a  guide  for  performing  routine  calculations 
of  the  sonic  boom  Intensity  and  distribution  on  the  ground  In  a  nonstand¬ 
ard  atmosphere. 

(C)  SUMMARY  OF  RESULTS  -The  effect  of  variations  in  temperature  and  pres¬ 
sure  on  the  sonic  boon  Intensity  and  lateral  distribution  were  established 
In  this  section  for  the  U.  S.  Standard  Atmosphere,  19&.  Variations  in 
these  properties  from  the  standard  conditions  and  the  presence  of  winds 
were  also  Investigated.  Slnqpllfled  methods  for  rredlctlng  the  influence 
of  these  changes  were  developed.  Further,  these  investigations  Indicated 
that  the  effect  of  varying  meteorological  conditions  for  flight  at  Ifech 
numbers  above  I.3  on  the  sonic  boom  intensity  Is  generally  no  more  than 
^5  percent  from  that  generated  In  the  still,  standard  atmosphere.  Flight 
at  Mach  ntanbers  below  I.3  may  result  In  more  significant  variations  In 
the  overpressure. 
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SECTION  m  A  REVIEW  OF  CONDITIONS  WHICH  MAY  CAUSE  ANOMALOUS 
PROPAGATION 

Various  laeteorologlcal  conditions  may  exist  between  the  airplane  and 
the  ground  which  will  cause  unusual  Influences  on  the  sonic  boom  Intensity 
and  Its  lateral  distribution  at  the  ground.  It  is  the  purpose  of  this 
section  to  discuss  these  conditions,  their  effect  on  the  oveipressure, 
and  when  possible,  to  present  criteria  for  determining  the  conditions 
necessary  to  produce  the  anomaly.  ISie  presence  of  the  required  condi¬ 
tions  may  result  In  fooislng  (local  Intensification  of  the  boom),  com¬ 
plete  cut-off  (no  boom  heard  on  the  ground),  extreme  lateral  spread 
(boom  heard  at  extremely  large  distances  to  the  side  of  the  fli^t  tr~..ck), 
ud  distortion  of  the  wave  form  due  to  Interaction  with  turbulence. 
Variations  In  temperature  alone  can  cause  focusing  or  cut-ol^,  while 
variations  In  both  wind  and  teeperature  can  caxise  focusing,  cut-off,  or 
extreme  lateral  spread.  For  purposes  of  simplicity  It  Is  convenient  to 
consider  these  In  terms  of  temperature  and  wind  produced  anomalies.  Hie 
effects  of  low  altitude  turbulence,  the  shock  wave  history  near  cut-off, 
and  the  effect  of  high  altitude  turbulence  and  shower  clouds  are  also 
discussed. 

(A)  TEMPERATURE  PRODUCED  ANOMALIES  -  The  primary  Influence  of  tem¬ 
perature  variations  in  the  atmosphere  is  to  distort  the  ray  path,  which 
describes  the  path  of  the  shock  front  from  the  airplane  to  the  ground. 

In  some  situations,  the  temperature  variation  can  be  such  that  the  rays 
are  refracted  so  completely  that  the  boom  will  not  reach  the  ground. 

This  condition  Is  known  as  "complete  cut-off"  of  the  sonic  boom.  On  the 
other  hand  the  conditions  may  be  such  that  the  boom  will  be  focused  (l.e. 
locally  Intensified)  at  a  point  on  the  ground.  Both  of  these  phenomena 
occur  at  very  low  Mach  numbers  where  the  speed  of  the  airplane  relative 
to  the  ground  Is  nearly  equal  to  the  local  sound  speed  at  the  ground. 

(1)  Compittt  Cut-oH  -  Investigation  of  the  ray  path  equations  (Eqs. 
II. la)  through  (II. If)  In  the  Appendix)  Indicates  that  for  conditions  o'^ 
no  wind  the  sonic  boom  will  be  prevented  from  reaching  the  ground  if  the 
velocity  of  the  airplane  relative  to  the  groiuid  Is  less  than  the  maximum 
speed  of  sound  et  any  altitude  between  the  airplane  and  the  ground. 

Stated  in  terms  of  the  airplane  Mach  number,  M  <  3 mix/ limi¬ 
ting  value  of  Mach  number,  ,  is  determined  by  putting  this  in  the 

form  of  an  equality  which  Is  expressed  in  Eq.  (8)  below. 


where 


M 


cut-off 


a. 


Eq.  (8) 


M  „  «  Largest  airplane  Mach  number  at  which  complete  cut-off 

°  will  occur. 


a„.,  =  Largest  value  of  sound  speed  between  the  airplane  and 

the  ground. 


sound  speed  at  the  airplane. 
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This  relationship  also  defines  the  meteorological  conditions 
required  for  no  boom  to  reach  the  ground.  It  has  been  applied  in  Fig.  20 
to  several  examples  to  Illustrate  the  application  to  the  standard  and 
nonstandard  atmospheres.  The  model  atmospheres  considered  are  illustrated 
in  Fig.  7*  The  figure  shows  that  for  altitudes  above  36>000  feet,  the 
limiting  Mach  number  lies  generally  between  1.1  and  1.2.  No  boom  would 
be  beard  on  the  ground  for  flights  in  a  still  atmosphere  at  Mach  numbers 
less  than  the  at  the  alrplan'  altitude. 


(2)  Fecutini  on  tli«  Ground  -  Focusing  or  intensification  of  the  sonic 
boco  C.V..  ressure  will  occur  if  a  set  of  adjacent  ray  paths  describing 
the  prcpag.ition  of  a  portion  of  the  shock  front  tend  to  cane  together, 
nils  leads  to  cusping  or  folding  over  of  the  shock  front  at  a  point  where 
the  area  between  the  adjacent  rays  (i.e.  ray  tube  area)  tends  to  go  to 
zero.  This  phenooena  is  discussed  in  more  detail  in  Ref.  36.  A  general 
expression  for  the  variation  of  the  area  between  the  adjacent  rays  is 
given  in  the  Appendix  (Eq.  (II-6)).  nils  expression  has  been  investi¬ 
gated  to  determine  the  meteorological  conditions  required  for  the  ray 
tube  area  to  approach  zero.  It  was  found  that  if  focusing  takes  place, 
it  can  occur  only  at  the  location  of  a  cut-off  of  the  sonic  boom  (either 
lateral  or  under  the  flight  track).  For  a  still  atmosphere  focusing 
can  occur  only  under  the  airplane  and  takes  place  simultaneously  with  the 
complete  cut-off  at  that  point.  The  Mach  number  at  vdilch  focusing  may 
occur  on  the  ground  is  given  by  Eq.  (9)  which  is  derived  in  Appendix  V. 

a, 

M_  =  ^  if  and  only  if  a  =  a  (Eq.  9) 

Focua  a,  g  max  '  ^ 


sound  speed  at  the  ground 
sound  speed  at  airplane 

largest  value  of  sound  speed  between  the  airplane  and  the 
ground 

If  a  temperature  inversion  exists  near  the  ground  such  that 

_  the  boom  cannot  be  focused  at  the  ground  because  the  cut-off, 

if  it  exists,  will  occur  at  the  top  of  the  inversion.  Thus,  in  an  atmos¬ 
phere  such  as  B-5  (Fig.  7)  the  boom  will  not  be  focused  at  the  ground, 
regardless  of  Mach  number.  Referring  to  Fig.  20,  the  line  shown 

for  the  Standard  Atmosphere,  and  model  A-3  would  also  represent  the 
^ focua  effect  of  approaching  the  Mach  number  at  which  focusing 

may  occur  is  shown  in  Fig.  21  where  the  variation  in  overpressure  pro¬ 
duced  at  a  fixed  altitude  by  varying  the  Mach  number  is  shown  for  the 
Standard  Atmosphere  and  Model  A-3.  The  data  for  these  curves  were  com¬ 
puted  using  the  method  of  Refs.  6-8. 
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Fig.  21  EHmcf  of  Approaching  on  Sonic  Boom. 


The  figure  Bhows  thr.t  the  boom  magnitude  may  eubstantlelly 
increase  as  the  cut-off  Mach  number  is  approached.  However,  two  factors 
should  be  noted  in  considering  these  results.  First,  except  for  continued 
flight  at  this  Mach  number  the  focusing  occurs  at  only  one  point  on  the 
ground  imder  the  fll^t.  track  and  thus  would  affect  only  a  very  small 
area.  Second,  and  perhaps  more  important,  is  the  interpretation  of  these 
results.  At  the  present  time,  a  factor  is  applied  to  all  sonic  boom 


estimates  to  accoxint  for  the  reflection  of  the  oblique  shock  vrave  frcan 
the  ground.  Hils  factor,  ,  (Eq.  (l))  is  us\ially  taken  as  approxi¬ 
mately  2.0.  However,  at  cut-off  the  shock  front  is  normal  to  the  ground 
and  there  is  no  reflected  wave  from  the  gro\md.  The  value  of  Kh  Is 
not  clearly  defined  under  these  circ\mistances.  This  aspect  of  the  pro¬ 
blem  is  considered  in  more  detail  in  Section  III.D.l  but,  to  summarize 
here,  the  conclusion  is  that  KaHiust  \rary  between  approximately  2.0  for 
the  oblique  shock  wave  to  approximately  1.0  for  the  normal  wave  at  cut¬ 
off.  In  this  respect  the  overpressure  experienced  at  the  ground  near  an 
atmospheric  focus  may  be  only  slightly  larger  than  that  predicted  for 
Mach  numbers  higher. than  that  required  for  focusing,  where  Kr«*  2.0.  As 
will  be  shown  in  the  following  sections,  the  seme  argument  would  hold  for 
wind  Induced  focusing  as  this  too  occurs  only  at  a  cut-off  where  the 
shock  front  is  normal  to  the  ground. 

In  making  routine  calculations,  Eqs.  (8)  and  (9)  should  be 
checked  to  determine  the  Mach  number  at  which  cut-off  or  focusing  nay 
occur.  If  an  inversion  exists  near  the  ground  Eq.  (9)  may  be  ignored. 
When  the  Mach  number  of  interest  lies  very  close  to  Mach  number  required 
for  focusing,  (the  proximity  may  be  estimated  from  Fig.  21)  the  approxi¬ 
mate  methods  outlined  in  Section  II  may  not  be  sufficient  to  adequately 
determine  the  boom  strength  on  the  ground.  If  an  estimate  is  required 
for  these  Mach  numbers  more  sophisticated  methods,  such  as  those  detailed 
in  Refs.  6-8,  must  be  applied. 

(B)  WIND  PRODUCED  ANOMALIES  -  Wind  sheeirs  between  the  airplane  and 
the  ground  will  distort  the  shock  front  as  it  travels  through  the  atmos¬ 
phere.  In  some  situations  these  shears  may  either  prevent  the  boom  from 
reaching  the  ground  or  significantly  increase  the  boom  lateral  extent. 

In  others  they  may  intensify  or  focus  the  boom  locally.  The  meteorologt 
call  conditions  necesseiry  to  produce  these  phenomena  are  presented  in  this 
section.  These  results  should  be  used  in  conjunction  with  routine  esti¬ 
mates  to  isolate  those  situations  where  unusual  propagation  may  occur. 

It  may  be  noted  that  with  the  exception  of  leurge  increases  in  location  of 
lateral  cut-off,  wind  produced  anomailieB  sure  generally  confined  to  the 
low  Mach  number  flight  regime.  At  higher  Mach  numbers  (usually  above 
1.3)  the  wind  shears  required  to  produce  the  unusual  effects  are  too 
large  to  be  realistic. 

(1)  Comp/tt«  Cuf-o// -  Investigation  of  the  ray  propagation  equa¬ 
tions,  (Appendix  Equations  (ll-la)  through  (ll-lf))  indicates  that  for  an 
atmosphere  with  wind,  no  boom  will  reach  the  ground  if  the  ray  directly 
under  the  airplane  is  refracted.  Experimental  evidence  of  complete  cut¬ 
off  has  been  noted  in  Ref.  4  for  very  low  Mach  numbers  (near  1.2).  In 
the  simplest  case,  this  requires  that  the  airplane  speed  relative  to  the 
ground  be  less  than  the  speed  of  sound  at  the  ground.  The  maximum  Mach 
nimber  at  wiilch  cut-off  would  occur  is  given  by  Eq.  (lO)  below. 

(a+U)  -U- 

^ - -  Eq.  (10) 
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where 

M  „  ■  largest  airpi’^e  Mach  number  at  which  complete  cut-off 

“  will  occur 

a  a  sound  speed  at  some  level  between  the  airplane  and  the 

ground 

U  a  tailwind  speed  component  at  the  same  level  as  selected 

for  a. ( U  is  negative  if  it  is  a  headwind  component ) 

U,  a  tailwind  comijonent  at  airplane  (U,  is  negative  if  it 

is  a  headwind  component) 

a^  a  souni  speed  at  airplane 

(a+U)  ®  largest  value  of  sound  speed  and  wind  component  speed 

m*x  which  occurs  between  the  airplane  and  the  ground 


The  calculation  of  M  cut-off  is  considerably  simplified  if 
the  combined  sound  speed  and  wind  speed  for  the  atmospheric  model  being 
considered  is  plotted  against  altitude  from  the  ground.  This  step  aids 
in  selection  of  the  maximum  value  of  (a+U)  tcot  lies  between  the  air¬ 
plane  and  the  ground. 

Several  examples  of  M  cut-off  been  calculated  using  the 

atmospheric  models  from  Fig.  7  ond  the  wind  models  from  Fig.  11.  These 
results  are  shown  in  Fig.  22  for  both  headwind  and  tailwind.  The*  cut-off 
Mach  number  with  no  wind  is  shown  for  reference. 

The  figure  shows  that,  in  general,  headwind  components  in¬ 
crease  the  cut-off  Mach  number  above  that  for  the  no  wind  case,  while 
tailwind  components  decrease  the  cut-off  Mach  number.  However,  in  the 
case  of  very  high  speed  winds,  the  tailwind  components  may  increase 
the  ;ut-off  Mach  number  for  altitudes  above  the  maximum  wind  speed.  This 
is  Illustrated  in  the  case  of  the  high-speed  jet.  It  should  also  be 
noted  that  the  altitude  at  which  complete  cut-off  occurs  is  indicated  by 
M  cut-off  “  In  "the  case  of  the  high-speed  jet  tailwind  in  the  Stand¬ 

ard  Atmosphere  for  flight  at  altitudes  above  28,000  feet,  the  cut-off 
would  occur  at  about  or  above  28,000  feet. 

(2)  Focusing  (/n</*r  fh#  F/igfcf  Troci -Under  some  conditions,  winds  may 
cause  focusing  of  the  sonic  boom  UD.der  the  flight  track.  The  mechanism 
is  similar  to  that  described  in  Section  III. A. 2.  It  has  been  found,  up¬ 
on  investigation  of  the  ray  tube  area  expression,  (Appendix  Eq.(lI-6)) 
that  focusing  under  the  fll^t  track  and  complete  cut-off  must  occur 
simultaneously.  The  boom  at  the  ground  may  be  intensified  if  the  cut-off 
occurs  at  the  ground.  The  Mach  number  at  which  focusing  may  occur  on  the 
ground  under  the  flight  track,  in  the  presence  of  wind,  is  given  by  Eq. 
(ll)  which  is  developed  in  Appendix  V. 
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Flf,  22  Efhet  of  Wind  on  Cut-Off  Mach  Humbar. 


If  the  quantity  (a-t-U)  is  not  largest  at  the  ground,  cut¬ 
off  and  possible  focusing  will  occtir  at  the  level  where  this  quantity 
has  its  greatest  VEilue.  This  is  slmilsur  to  the  effect  of  an  inversion  in 
the  still  atmosphere.  Referring  to  Fig.  22,  the  Mach  cut-off  lines  for 
the  headwind  and  tailwind  would  also  represent  the  Mpocu*  lines.  An 
interesting  example  of  cut-off  above  the  ground  with  no  focusing  possible 
at  the  ground  is  Illustrated  quite  vividly  in  the  case  of  the  Standard 
Atmosphere  and  the  high-speed  Jet  tailwind  profile.  For  airplanes  flying 
above  28,000  feet  at  Mach  numbers  near  ,  the  cut-off  would  occur 

at  28,000  feet,  thus  precluding  the  possibility  of  focusing  at  the  ground. 

The  effect  of  wind  on  the  overpressure  under  the  airplane 
was  investigated  by  the  method  of  Refs.  6-8, in  che  Standard  Atmosphere 
for  winds  increasing  lineeurly  from  zero  at  the  ground  to  the  maximum 
value  at  the  airplane.  The  wind  magnitude  at  the  airplane  was  increased 
until  cut-off  and  focusing  at  tht'  ground  were  achieved.  These  results 
sure  shown  in  Fig.  23  where  the  overpressure  with  wind  has  been  divided 
by  that  obtained  without  wind  for  the  airplane  at  a  fixed  Mach  number 
and  altitude. 


The  figure  shows  that  the  focusing  effect  occurs  at  wind 
speeds  very  near  those  required  to  produce  complete  cut-off.  The  local 
intensification  may  not  be  as  high  as  it  appears  at  first  sight  because 
of  the  possible  variation  in  ground  reflection  factor  K^with  shock  wave 
angle.  This  consideration  is  discussed  more  fully  in  Section  III.D.l. 

Another  interesting  consideration  is  that  wind  induced  fo¬ 
cusing  can  occur  only  at  the  very  low  Mach  numbers.  This  may  be  seen  if 
Eq.  (11)  is  rearranged  so  that  the  wind  required  for  cut-off  and  focusing 
at  the  ground  may  be  calculated  as  a  function  of  the  airplane  Mach  num¬ 
ber.  An  especially  simple  example  may  be  found  in  the  Standard  Atmos¬ 
phere  with  a  wind  varying  linearly  from  zero  at  the  ground  to  a  maximum 
at  the  airplane  and  by  assuming  that  the  airplane  is  above  36,000  feet. 
The  wind  required  for  cut-off  and  focusing  has  been  determined  for  this 
case  and  is  shown  plotted  in  Fig.  24. 
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The  figure  shows  that  the  headwind  speeds  at  the  airplane 
required  to  produce  cut-off  and  focusing  are  in  excess  of  200  feet  per 
second  for  Mach  numbers  above  1.3»  Wind  velocities  In  the  atmosphere 
rarely  exceed  this  vedue  at  altitudes  above  40,000  feet,  where  airplanes 
such  as  a  large  supersonic  transport  would  fly  at  these  Mach  numbers.  At 
higher  Me,ch  numbers  the  wind  velocity  required  for  cut-off  and  focusing 
would  be  larger  them  emy  which  could  physically  exist  in  the  atmosphere. 

(3)  Focusing  to  th»  Sid*  of  th*  Flight  Track  -  As  In  the  case  of  focusing 
under  the  flight  track,  intensification  of  the  boom  to  the  side  of  the 
fll^t  track  may  occur  simultaneously  with  the  lateral  cut-off.  Investi¬ 
gation  of  the  ray  tube  area  expression  in  the  Appendix  (Eq.(II-6))  leads  to 
em  expression  of  the  wind  required  to  produce  these  phenomena.  Die  ex¬ 
pression  for  the  required  wind  conditions  is  derived  in  Appendix  VI  of 
this  repoil;.  The  result  has  been  applied  to  several  examples  which  are 
presented  in  Fig.  25. 
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ag  =  SOUND  SPEED  AT  GROUND 
aa  =  SOUND  SPEED  AT  AIRPLANE 

NOTE;  WIND  SPEEDS  ARE  MEASURED  RELATIVE 
TO  THOSE  CM  THE  GROUND  (SEE  TEXT) 


The  required  wind  speed  components  are  measured  relative 
to  those  on  the  ground.  For  example,  the  sidewind  speed,  V*  would  be 
the  absolute  value  of  the  algebraic  difference  of  the  sidewind  speed  at 
the  ground  and  that  at  the  airplane,  i.e.  |Vj  -Vg  I-  Wind  is  a  vector 
quantity  so  that  the  direction  must  be  accounted  for  by  the  sign  of  the 
number.  The  following  directions  were  defined  in  developing  the  above 
curves ; 

•  Winds  in  the  direction  of  the  flight  path  are  positive  (i.e.  tailwinds). 

•  Winds  in  the  opposite  direction  of  the  flight  path  are  negative  (i.e. 
headwinds). 

•  Winds  coming  from  the  rl^t  side  nomal  to  the  flle^t  path  are  positive. 

•  Winds  coming  from  the  left  side  nomal  to  the  fligh+  path  are  negative. 

The  figure  indicates  the  wind  components  required  to  cause 
off-track  focusing  at  several  Mach  numbers.  It  is  evident  that  this 
phenomenon  is  restricted  to  the  very  low  Mach  numbers  because  of  the 
magnitudes  of  the  wind  speeds  required  at  the  higher  Mach  numbers. 

To  indicate  the  proximity  to  the  required  value  of  wind 
which  may  cause  an  unusual  effect,  several  cases  have  been  computed  in 
the  Standard  Atmosphere.  The  sonic  boom  strength  at  the  lateral  cut-off 
point  was  computed  by  the  method  of  hefs.  6-8  for  various  wind  speeds 
and  directions.  These  results  were  then  divided  by  the  boom  strength 
under  the  flight  track  for  no  wind  in  order  to  better  illustrate  the  var¬ 
iation.  llie  effect  of  winds  on  the  boom  strengt.h  at  the  latered  cut-off 
point  is  shown  in  Fig.  26. 

The  figure  shows  that  in  some  circumstances  and  at  very  low 
Mach  numbers  the  overpressure  to  the  side  of  the  flight  track  may  exceed 
that  under  the  flight  track.  The  same  phenomenon  may  also  be  seen  in 
Fig.  15  for  the  standard  atmosphere  with  the  high  speed  Jet  sidewind  pro¬ 
file  for  Mach  1.2.  Here  again,  the  variation  of Kr,  the  reflection  fac¬ 
tor,  with  shock  angle  must  be  considered  as  the  focusing  effect  occurs 
simulteuieously  with  cut-off  where  the  shock  is  nesurly  normal  to  the  ground. 
This  is  discussed  in  more  detail  in  Section  III.D.l. 

A  better  indication  of  the  range  of  wind  speeds  which  are 
required  to  produce  significant  off-track  focusing  may  be  obtained  by 
a  cross-plot  of  the  data  in  Fig.  26  for  Mach  1.2  on  the  plot  of  Fig.  25. 

The  winds  which  will  produce  overpressures  at  the  lateral  cut-off  po¬ 
sition  equed  to  these  normally  produced  under  the  flight  track  without 
wind  were  picked  to  form  the  boundary.  Thus,  the  wind  range  which  will 
produce  this  effect  is  shown  as  the  cross  hatched  area  in  Fig.  27. 
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V  ~  (SIDEWIND  AT  GROUND  -  SIDEWIND  AT  AIRPLANE),  FPS. 


F/g.  26  Efftet  of  WinJi  on  Sonic  Boom  at  Lotoraf  Cut-Off, 


For  Mach  numbers  greater  than  1.3>  the  relative  wind  speeds 
required  to  produce  the  off-track  focusing  become  very  large.  Tbis  is 
illustrated  when  comparing  the  computed  data  for  Mach  1.2  and  Mach  I.3 
in  Fig.  15  for  the  high  speed  Jet  sidewind.  The  data  at  Mach  I.3  shows 
that  the  influence  of  the  high  speed  winds  is  practically  negligible. 

For  most  practical  cases  of  transonic  flight  (Mach  numbers  near  1.0)  at 
high  altitudes  (above  h0,000  feet)  the  winds  at  the  airplane  would  not 
be  of  sufficient  magnitude  to  produce  the  off-track  focusing  effect.  Thus, 
the  primary  consideration  would  be  one  of  focusing  under  the  flight  track 
at  or  near  the  cut-off  Mach  number.  The  probability  of  occurrence  of 
the  focusing  phenomenon  either  under  or  to  the  side  of  the  flight  track 
is  considered  in  more  detail  in  Section  V. 
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In  considering  routine  calculations  of  the  sonic  hoom  the 
wind  speed  and  sound  speed  variation  hetween  ti;e  airplane  and  the  ground 
should  he  determined.  The  criterion  developed  in  this  and  the  preceding 
section  should  he  checked  In  the  low  Mach  number  rEuige  to  see  If  the 
possihillty  of  focusing  exists  for  flight  at  the  altitudes  of  Interest. 
If  this  possibility  does  exist  more  sophisticated  methods,  such  as  that 
of  Refs.  6-8,  must  he  used  to  determine  the  hoom  strength  at  the  ground. 
Otherwise,  the  methods  outlined  In  Section  II  should  he  sufficient  to 
determine  the  hoom  strength  on  the  ground. 

(4)  Extreme  Loteral  Spread  -  In  some  c'ses  winds  may  Increase  the 
location  of  the  lateral  cut-off  to  very  laige  distances  from  the  fllg^it 
track.  This  phenomenon,  which  has  been  noted  here  as  extreme  lateral 
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spread,  nay  occur  at  all  Mach  nvonbers.  An  Investigation  ol'  Appendix  Bqs. 
(ll-la)throu^(II-lf )  vhlch  describe  the  -^y  path  propagation,  has  led  to 
an  expression  for  the  winds  required  to  j  ’oduce  the  extrenx!  spread.  The 
relationship  Is  given  In  Eq.  (12)  and  Is  aeveloped  In  detail  In  the  Ap¬ 
pendix  VII. 


where 

U,  ■  Tailwind  speed  at  airplane  (negative  if  headwind) 

V,  •  Sidewind  speed  at  airplane  (positive  if  coming  from  right 

side  of  flight  track 

M  ■  Airpleine  Mach  number 

■  Sound  speed  at  airpleuie 
U  ■  Tailwind  speed  at  height,  H  ,  above  ground 

V  “  Sidewind  speed  at  height, H,  above  ground 

a  ■  Sound  speed  at  height, H,  above  ground 

The  (+)  sign  on  the  term  Vm ^-1  V  is  used  when  investigating 
extreme  lateral  spread  on  the  left  side  of  the  flight  track  emd  r.he  (-)  sign 
is  used  when  investigating  the  possibility  on  the  right  side. 

The  equation  shows  that  for  the  winds  to  produce  extreme 
latersJ.  spread  the  value  of  the  queuitity  A  at  the  airplane,  i.e.  Aa,  must 
be  greater  than  at  any  other  point  between  the  airplane  emd  the  ground. 

The  application  of  Eq.  (12)  to  a  specific  case  will  illustrate  its  use. 

The  exeimple  chosen  was  the  Stemdard  Atmosphere  with  the  high-speed  Jet 
profile.  Airplane  Mach  number  was  chosen  as  2.0,  emd  various  wind  direc¬ 
tions  were  selected.  The  right  side  of  Eq.  (12)  was  computed  for  each 
assumed  direction  emd  is  shown  plotted  in  Fig.  28. 
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Hm  flgor*  ■)  m's  the  eltltude  renge  and  wind  directions  for 
which  extreae  lateral  spn  v.  z:::y  occur.  It  should  be  noted  that,  of  the 
wind  directions  selected,  the  0*  (tailwind),  90*  (sidewind),  and 
60*  wind  directions  could  cause  this  phenomenon,  and  then  only  If  the 
airplane  Is  below  26,000  feet.  The  former  direction  (tailwind)  would 
cause  the  extreae  spread  on  both  the  rlj^t  and  left  sides  of  the  fll^ 
track  while  the  latter  two  dlrectiuns  (60*  and  90*)  would  cause  the  ex> 
treat  spread  mUy  on  the  left  side.  In  this  example,  the  extreae  lateral 
spread  would  not  occur  for  an  airplane  flying  above  28,000  feet.  Tbd 
winds  used  In  the  exaaple  In  Fig.  26  are  characteristic  of  acst  wind 
profiles  In  the  ataosphere.  nuis,  the  problem  of  large  Increases  of  the 
lateral  cut-off  location  would  not  exist  for  flight  at  altitudes  above 
the  level  of  maxi ana  wind. 

The  effect  of  hl^  winds  on  the  location  of  lateral 
cut-off  was  Investigated  for  a  soaewhat  acadealc  case,  but  the  results 
were  characteristic  of  the  phenomenon.  The  overpreesure  on  the  ground 
was  calculated  by  the  method  Refs.  6-8  for  a  wind  profile  which  varied 
linearly  frca  sero  at  the  ground  to  a  naxlmum  at  the  airplane  In  the 
Standard  Ataosphere.  The  asjpltude  of  the  wind  at  the  airplane  was  In¬ 
creased  until  It  exceeded  tha^  requlreu  to  produce  the  extreae  lateral 
spread.  The  results  of  this  calculation  are  shown  In  Fig.  29. 

The  figure  shows  that  althou^  the  wind  Increases  the  lo¬ 
cation  of  the  lateral  cut-off  the  overpressures  In  that  region  are  quite 
low.  For  this  reason  routine  calculations  of  the  overpressure  on  the 
ground  al^  Ignore  this  effect  as  the  overpressure  beyond  the  normal  cut¬ 
off  location  drops  off  qjulte  rapidly  with  distance  froa  the  flight  track. 
In  any  case,  as  was  noted  earlier,  this  phenomenon  will  probably  be  un- 
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(C)  EFFEC' 5  OF  LOW  ALTITUDE  TURBULENCE  -  The  previous  discussion  was 
limited  to  the  study  of  the  effect  of  a  horizontally  stratified  atmos¬ 
phere  on  the  propagation  of  the  sonic  boom.  In  this  section,  it  will  be 
shown  that  any  study  of  the  meteorological  aspects  of  the  sonic  boom  can 
be  separated  into  effects  due  to  the  horizontally  stratified  atmosphere 
and  effects  due  to  low  level  turbulence.  The  turbulent  process  in  the 
atmosphere  is  the  result  of  some  form  of  instability.  This  may  be  either 
a  result  of  mechanical  instabilicy,  such  as  produced  by  wind  shear,  or 
flow  over  obstacles,  or  thermal  instability  such  as  produced  by  solar 
heating  of  the  ground.  These  forms  of  instability  produce  random,  turbu¬ 
lent  fluctuations  in  wind  and  temperature  which  can  only  be  studied  and 
described  in  statistical  terms  such  as  the  correlation  coefficient  and 
spectrum  analysis.  Consequently,  when  the  effects  of  a  turbulent  at¬ 
mosphere  on  the  sonic  boom  are  being  siudie  .  it  is  ncccssarj'  to  ucc 
statistical  parameters  for  the  development  of  a  theory  of  turbulent  ef¬ 
fects  on  the  sonic  boom  vdiich  will  be  introduced  in  this  section.  The 
theory  describes  either  the  turbulent  temperature  or  the  turbulent  wind 
effects  on  the  sonic  boom.  The  effects  of  ccmbinod  temperature -wind 
turbulence  are  currently  being  studied. 


(1)  Atmo*ph«rie  Stiuetur#  -  llie  flow  patterns  Of  the  atmosphere  can  be 
visualized  as  being  ccsnposed  of  oscillations  varying  in  size  from  3  to  k 
thousand  miles  in  wave  length  to  ones  which  are  locally  on  the  order. of  a 
few  feet.  In  order  to  ascertain  the  effect  of  the  wave  lerigth  of  at¬ 
mospheric  properties  on  the  propagation  of  the  sonic  boon,  a  harmonic 
analysis,  or  more  rigorously  a  power  spectrum  analysis  of  data  containing 
these  socillations,  can  be  made. 


The  power  spectxtim  of  wind  shown  in  Fig.  30  ■‘■s  a  typical 
plot  of  kinetic  energy  of  the  horizontal  wind  speed  (which  is  obtained 
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Fig.  30  Horizontal  Wind—Spood  Spoctrum  at  Braakhcven  National  Laboratory  at  abiiut 
100m  H tight  (Rtf.  38). 


from  the  spectral  analysis)  versus  the  frequency  of  the  waves  in  cycles 
per  hour.  The  tes^^eraturc  cpectnaa  cmtvc  is  very  sisiilar  in  fciui  to  the 
horizontal  wind  spectrum.  The  following  results  may  be  deduced  from 
Fig.  30: 

•  There  is  one  maximm  in  the  wave  energy  spectrum  neEu*  100  hours  per 
cycle.  This  represents  waves  due  to  large  scale  weather  systems. 

•  A  secondaiy  maximm  occurs  near  12  hours  per  cycle  which  may  be 
identified  with  the  diurnal  cycle  of  meteorological  conditions. 

•  There  is  a  broad  minimum  centered  near  a  period  of  one  hour  per  cycle. 
This  flat  part  of  the  curve  (the  so-called  spectral  gap)  means  that 
there  are  very  few  eddies  of  a  size  from  1  to  20  miles  in  the  atmos¬ 
phere  when,  for  example,  the  average  wind  speed  is  10  mph. 

•  A  third  maximian  is  found  near  periods  of  one  minute.  This  means 
that,  for  the  same  average  wind  speed  of  10  mph,  there  are  many  turbu¬ 
lent  oscillations  with  a  wave  length  of  about  900  feet  imbedded  in  the 
passing  atmospheric  flow  field. 

These  data  show  that  the  study  of  the  effects  of  the  at¬ 
mosphere  on  the  propagation  of  the  sonic  boom  can  be  separated  into  two 
parts.  One  part  of  the  study  is  the  effect  of  the  stratified  atmosphere, 
in  which  the  horizontal  gradients  on  wind  and  temperature  may  be  neglec¬ 
ted  because  they  are  small  relative  to  the  distances  (20  -  4o  miles)  ef¬ 
fected  by  the  sonic  boom.  As  was  pointed  out  above,  eddies  of  this  size 
are  rare  in  the  atmosphere.  The  second  part  of  the  study  are  the  effects 
of  turbulent  temperature  and  wind  fluctuations  over  small  periods  of  time 
and  space  on  the  sonic  boom.  These  fluctuations  are  superimposed  upon 
the  relatively  slowly  varying  wind  and  temperature  patterns  of  the  large 
weather  systems,  which  for  purposes  of  studying  the  turbulence  are 
treated  as  being  nearly  stationary.  The  magnitude  and  frequency  of  the 
short  period  variations  depend  upon  such  parameters  as  time  to  day,  wind 
speed,  cloud  cover  and  the  lapse  rate  of  temperature  with  altitude. 

The  structure  of  a  horizontally  stratified  atmosphere  is 
obvious,  but  the  description  of  a  tiurbulent  atmosphere  requires  the  in¬ 
troduction  of  statistical  concepts  which  may  be  unfamiliar.  When  a  flow 
is  characterized  as  being  tiirbulent,  it  is  implied  that  fields  of  irregu¬ 
lar  and  random  fluctuations  of  scalar  quantities  (e.g.  temperature)  eind 
vector  quantities  (e.g.  wind)  occur  about  some  mean  value.  The  fluctua¬ 
ting  part  of  the  flow  field  can  be  regarded  as  the  superposition  of  a 
large  number  of  different  sized  oscillations  which  can  be  given  analytical 
form  by  a  three  dimensional  Fourier  analysis  of  the  velocity  field. 

To  find  a  quantity  which  is  easily  measured  and  yet  is  a 
good  statistical  measure  of  the  amount  of  energy  contained  within  each 
particular  eddy  size  of  wave  length,  Taylor  (Kef.  6T)  observed  that  the 
following  Fourier  transformation  pair  exists  between  the  correlation 
coefficient,  R|j  for  scalar  fields  (or  correlation  tensor  for  vector 
fields)  and  the  wave  spectnm  (tensor),  $jj 


Eq.  (13) 


while  ♦  ,  I  represents  the  energy  density  as  a  function  of  wave  number. 

The  turb\ilent  eddies  decay  slowly  and  are  unchanged  over 
the  distances  they  travel  when  carried  by  the  mean  atmospheric  flow  dur¬ 
ing  the  time  Intervals  used  in  correlating  the  measurements.  This  Im¬ 
plies  that  for  homogeneous  turbulence,  the  Eiilerlan  or  time  lag  correla¬ 
tions  may  be  tised  Instead  of  the  difficult  Lagranglan  or  space  correla¬ 
tion.  A  similar  analysis  mlg^t  be  made  of  the  cvoss- spectrum  between  the 
two  variables,  temperature  and  wind,  which  lead  to  a  measure  of  the  rate 
at  which  the  temperature  is  being  advected  or  carried  along  by  the  turbu¬ 
lence,  This  suggests  a  way  In  which  the  effects  of  combined  temperature 
and  wind  fields  may  be  Incorporated  Into  the  theory  of  turbulent  scatter¬ 
ing  of  shock  waves. 

Other  parameters  which  are  useful  In  describing  turbulence 
are  the  mean  wind,  IT  ,  and  the  RMS  wind  speed  (u  2)1/2.  Prom  dimensional 
analysis  (Refs,  Ul  and  42)  It  can  be  argued  that  the  spectrum  of  the  ver¬ 
tical  eddies  varies  with  hel^t  above  the  ground  In  such  a  mannei-  that  If 
the  frequency  of  the  eddies,  o) ,  Is  divided  by  the  mean  wind  u  and  mul¬ 
tiplied  by  the  hel^t,  z  ,  a  reduced  frequency  f  Is  obtained  which  is 
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eBsentlally  the  ratio  of  height  to  wave  length,  f=Y*  reduced  fre¬ 
quency  is  used  to  normalize  data  to  a  comnon  hase.^  The  "scale  of  tur- 
hulence"  vhlch  Is  roughly  l/lO  the  height  of  the  obstacles  to  the  flow, 
and  the  varleuice  of  the  longitudinal,  vertical  and  latitudinal  velocity 
ccBqponents  are  also  meaningful  parameters  vhlch  may  be  used  to  describe 
the  turbulent  variations  of  the  properties. 

(2)  Variation  in  Turbulant  Intantity  -  The  description  of  hOW  the  turbu¬ 
lence  decays  progressively  from  large  eddies  to  heat,  led  Kolmogoroff 
(Refs.  69  and  70)  to  formulate  a  similarity  hypothesis  In  19^1.  This 
hypothesis  states  that  the  Influence  of  the  large  eddies  on  the  smaller 
eddies  as  the  large  eddies  decay  to  smaller  Oues  diminishes  gradually 
and  conseqxiently,  small  eddies  tend  to  have  uniform  properties  for  all 
types  of  turbulence.  Therefore,  the  properties  of  these  small  eddies  are 
solely  determined  by  the  average  rate  of  dissipation  of  energy  per  unit 
mass.  This  hypothesis  leads  to  the -form  of  the  pover  spectrum  In  the 
Inertial  subrange  (e.g.  Fig.  30). 

To  Increase  the  tractablllty  of  equations  such  as  those  for  , 
diffusion  e«d  acoustic  propagation  susd  to  mitigate  the  problem  of  mea¬ 
surement  of  the  atmospheric  parameters.  It  Is  often  assumed  that  the 
tvirbulence  Is  homogeneous  and  Isotropic.  If  the  turbulence  quantita¬ 
tively,  has  the  same  structure  In  all  parts  of  the  flow.  It  Is  homo¬ 
geneous.  It  Is  Isotropic  If  Its  statistical  properties  are  Independent 
of  direction.  Homogeneity  and  Isotropy  are  usually  good  approximations’ 
near  the  ground  (in  the  atmospliere)  when  turbulence  Is  fully  developed, 
when  the  sky  Is  clear  and  the  winds  are  strong.  Because  of  the  mathe¬ 
matical  complexity  Inherent  In  describing  non-hcmogeneous  turbulence, 
all  Investigators  to  date  have  assumed  that  the  turbulence  Is  isotropic 
and  homogeneous  when  deriving  the  scauterlng  equations. 

Examination  of  turbulent  spectra  such  as  those  assembled 
In  lAanley  and  Panofsky  (Ref.  4l,  p.  I6I  ff)  Indicates  that  the  following 
significant  factors  about  atmospheric  txirbi^ence  near  the  ground  should 
be  borne  In  mind  when  studying  the  propagation  of  the  sonic  boom  In  low 
level  turbulence. 

•  Hhs  variance  of  vertical  wind  velocity  near  the  grc\aid  which  Increases 
with  the  square  of  the  wind.  Is  a  function  of  the  surface  roughness, 
and  Is  relatively  Independent  of  the  laps*  rate  of  temperature  with 
altitude. 

•  The  variance  of  the  lateral  component  (l.e.,  the  cross  wind)  Is  very 
sensitive  to  lapse  rate,  but  not  to  surface  roughness  or  wind  speed. 

•  Ihe  variance  of  the  longitudinal  con^onent  of  the  wind  (l.e.,  along 
the  wind)  depends  upon  the  lapse  rate  of  tecg>erature,  surface  rou^- 
nesB  and  mean  wind  speed. 

•  Generally,  the  variance  of  the  horizontal  wind  components  Is  two  to 
three  times  that  of  the  vertical  wind  (l.e.,  the  turbulence  Is  non- 
hcmogeneous  and  non- Isotropic). 

•  It  Is  Indicated  that  the  eddies  tend  to  be  elongated  along  the  mean 
wind  direction  when  the  winds  are  strong  (according  to  Lumley  and 
Panofsky  (Refs.  4l  and  42). 
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•  Examination  of  the  turbulent  spectra  presented  in  Bef.  shows  the 
maximum  value  of  the  energy  density  in  the  spectrum  of  vertical  velo¬ 
city  tends  to  occur  near  a  wave  length  which  is  four  times  the  heigdit 
above  the  ground,  while  the  analogous  teiqperature  spectrum  maximum 
occurs  near  a  wave  length  ten  times  the  height  above  the  ground.  !lhus, 
the  maximum  turbulent  energy  density  tends  to  occur  at  longer  wave 
lengths  with  increasing  altitude. 

Since  the  measurements,  to  date,  of  the  sonic  boom  have 
taken  place  either  over  deserts  or  cities,  it  is  Important  to  note 
that  there  is  little  data  on  the  homogeneity  or  Isotropy  of  tvirbulence 
from  that  kind  of  region.  However,  the  assumptions  of  hemogeneous  and 
isotropic  turbulence  should  yield  an  adequate  preliminary  solution 
which  will  Indicate  the  important  parameters  that  Influence  the  distor¬ 
tion  of  the  sonic  boom  signature. 

(3)  Scattering  of  Acoustic  Energy  by  Turbulence  -  Ihe  principal  investiga¬ 
tors  in  this  field  have  been  Perkeris  (Ref.  4$),  Blockhintzev  (Refs.  46 
and  47)  Lifi^thlll  (Ref.  48)  and  Kraichnan  (Ref.  49).  IRiey  have  been 
concerned  primarily  with  scattering  of  small  amplitude  sound  waves.  The 
problem  of  determining  how  acoustic  energy  is  redirected  (l.e.  scattered) 
by  the  interaction  between  acoustic  waves  and  a  field  of  turbulence  is 
difficult  and  because  of  the  Inherent  ccBq>lexlty,  the  simplest  case, 
that  of  homogeneous  and  isotropic  turbulence  (l.e.,  the  intensity  of  the 
turbulence  has  no  preferred  orientation  and  it  is  uniformly  distributed) 
has  been  the  only  case  studied  to  date. 

The  initial  attenq)t  at  analyzing  the  effects  of  turbulence 
on  the  sonic  boom  by  Palmer  (Ref.  81)  was  to  extend  the  analysis  of  the 
variation  of  the  aaq)lltude  due  to  turbulent  scattering  developed  by 
Tatarski  (Ref.  ^l).  Since  this  was  the  problun  of  greatest  immediate  in¬ 
terest. 


The  extension  of  this  analysis  had  the  advantage  of  having 
been  experimentally  verified  (Ref.  t>4)  in  the  lower  levels  of  the 
atmosphere  for  small  amplitude,  hi^  frequency  sound  waves.  Inherent 
disadvantages  where  the  necessity  of  assuming  the  constancy  of  the  struc¬ 
ture  of  atmospheric  t\irbulence  with  altitude,  of  ignoring  the  non-linear 
effects  associated  with  large  amplitude  waves,  and  of  neglecting  the 
effects  of  variations  in  the  turbulence  power  spectrum  as  a  function  of 
frequency. 


This  analysis  consisted  of  treating  the  "N"  wave  as  being 
composed  of  a  Poxirler  series  sum  of  frequencies  and  considering  one  wave 
from  the  resulting  chain  represented  by 


N(a))  =  ^  V 

"  n«  1 


n*i 


n 


sin  n 


0)  t 
0 


51 


w-rjwvu'j  MV  M'uW ‘j  wui  um  g  UJt  kA  mj^  m  n.  vjuvaj  lu't  luiwj  i 


MMfWV-W-WmV  ■ 


yvjniK  jViWUiQ^^WATVVWVtfVVWtf^  tfTIMTCtflUffljtf 


I 


vhere  is  the  amplitude  of  the  N-wave 

0).  is  the  frequency  of  the  fundamental  harmonic 
of  the  N-wave. 


Tartaraki  (Ref.  51)  has  shown  that  the  amplitude  fluctua¬ 
tions  of  a  sound  wave  that  is  propagating  throu^  a  turbulent  medium  are 
lognoxmally  distributed.  Kie  variance,  o  of  the  amplitude  fluctua¬ 
tions  of  a  monochromatic  acoustic  wave  is  then  given  by 


2 

O 


vdiere  0  is  a  geometrical  constant,  Cn  is  an  atmospheric  structiure 
function,  which  depends  on  tesqierature,  wind  cmd  the  angle  of  inci¬ 
dence  of  the  wave  upon  the  turbulent  region,  S  is  the  path  length  and  1 
is  the  wave  length.  The  N-wave  can  be  represented  as  a  Pourier  series, 
where  the  terns  denote  harmonic  sine  waves  of  varying  amplitude  and  fre¬ 
quency.  Ihe  total  variance  of  the  pressure  amplitude  of  the  N-wave  may 
be  found  as  the  sum  of  the  presstire  variances  of  each  harmonic  providing 
that  each  propeigates  Independently  of  all  others  (i.e.,  uncorrelated) 
throu^  the  turbulent  regions  of  the  atmos^ere.  If  the  harmonics  do 
net  propagate  Independently  the  lognormal  form  of  the  statistical  pres¬ 
sure  amplitude  distribution  is  still  preserved  (Ref.  63)  provided  that 
the  wave  length  of  the  N-wave  lies  within  the  atmospheric  inertial  sub- 
raxige  of  the  power  energy  spectrum  as  is  shown  in  Pig.  30. 

Since  the  parameters  which  enter  into  the  functional  form 
of  the  variance  are  functions  of  the  atmospheric  turbulent  structure, 
the  path  length,  the  wave  length  of  the  N-wave,  and  the  angle  of  propaga~ 
tlon,  it  follows  that  study  of  the  scattering  of  the  sonic  boom  should 
classify  the  data  according  to  a  scheme  based  upon  these  parameters.  In 
particular,  the  atmospheric  turbulent  structure  will  be  related  to  time 
of  day  and  cloud  cover. 

This  analysis  (Ref.  8I)  is  based  upon  the  scattering  of  a 
sinusoidal  train  of  small  amplitude  waves.  As  a  result  of  the  tests  at 
Qklahoma  City,  it  became  apparent  that  the  problem  of  the  interaction 
of  the  sonic  boom  with  turbulence  is  not  one  of  continuous  interaction 
of  a  wave  train  with  turbulence  which  is  unaffected  by  boundary  condi¬ 
tions.  It  is  rather,  one  involving  the  near  field  of  a  large  aji5)litude 
pulse  scattered  by  a  tvirbulent  field  of  temperature  and  wind  near  the 
earth's  svirface.  The  amplitude  statistics  that  were  developed  for  con¬ 
tinuous  waves  do  not  predict  the  angular  energy  distribution  of  the  scat¬ 
tered  energy  from  the  direction  of  propagation  and  does  not  accurately 
depict  the  non-linear  Interactions  which  occvir  under  these  conditions. 
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The  next  step  was  to  examine  the  theory  for  scattering  of 
aji  acoustic  pulse.  Batchelor  (Ref.  has  shown  that  the  following 
derivation  of  the  field  resulting  from  scattering  of  an  acoustic  pulse 
is  equivalent  to  that  developed  in  the  differing  forms  of  Pekeris  (Ref. 
1+5)^  Blockhintzev  (Ref.  k6  and  hj),  Llghthill  (Ref.  48)  and  Kraichnan 
(Ref.  49).  The  procedure  followed  in  formulating  the  theory  was  to  de¬ 
rive  the  wave  equation  from  the  Navler-Stokes  equation  and  the  equation 
of  continuity.  A  forcing  function  is  then  applied  to  this  equation,  and 
the  solution  of  the  resulting  partial  differential  equation  is  then 
accomplished  by  using  the  method  of  perturbations  and  the  Fourier  trans¬ 
form.  This  technique  gives  the  resultant  angular  distribution  of  scat¬ 
tered  energy  in  terms  of  a  differential  cross-section  vdiich  is  defined 
as  the  ratio  between  the  power  scattered  per  unit  solid  angle  in  a  given 
direction  to  the  incident  power  intensity  per  unit  area.  The  use  of 
the  Fourier  transform  in  solving  the  partial  differential  equation,  gives 
a  solution  for  the  angular  distribution  of  energy  which  can  be  expressed 
in  terms  of  a  power  spectrum  of  a  turbulent  quantity  such  as  is  described 
in  the  section  on  turbulence.  Section  III.C.2.  An  Important  concept 
used  in  the  derivation  is  the  scattering  vector,  k  ,  whi(^  is  defined 
as  the  vector  difference  between  the  vector  wave  number,  K  ,  representing 
the  incWent  wave,  and  a  unit  vector,  /  ,  multiplied  by  ^e  absolute  val¬ 
ue  of,  IKI,  vdilch  represents  the  scattered  wave,  K  (i.e.  K  =  |K|  J"). 

For  instance,  if  a  plane  wave  front  such  ps  an  acoustic 
pulse  is  Incident  upon  the_scattering  volume,  the  energy  will  be  re¬ 
directed  along  the  vector  /  ,  At  a  sufficiently  large  distance  from  the 
scattering  volume,  there  will  be  an  interaction  between  the  incident  wave 
and  the  scattered  wave.  The  propagation  of  the  disturbance  which  results 
from  this  interaction  is  represented  by  the  scattering  vector  It  . 
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SCATTERED  SPHERICAL  WAVE 


The  derivation  of  the  equations  governing  tvirhulent  scat¬ 
tering  is  as  follows: 


The  Navier-Stokes  Equation: 


\9 1  9  Xj  /  9  Xj  3  9xj  I  9  Xk 


9  Ui  ,  9  Uk 

Eq.  (17) 


and  the  equation  of  continvilty: 


+  (C  u,) 
9t  9  Xi 


Eq.  (l8) 


are  differentiated  partially  with  respect  to  x,  and  t  in  order  to  elimi¬ 
nate  the  time  deiwndent  velocity  terms,  where 

Jjj  is  Dirac  delta  =  1  when  1  ■=  J  and  Is  zero  otherwise 

Q  is  the  mass  density 

t  is  time 

Uj  is  the  velocity  component  in  the  Xi  - th coordinate,  i  =  1,  2,  3 

p  is  the  pressure 

p  is  the  coefficient  of  viscosity. 

Neglect  of  second  order  terms  as  a  result  of  assianlng  that  the  pressure, 
tenqierat^ire  and  density,  individually  are  the  sum  of  a  mean  plus  a  small 
perturhation  (this  procedure  is  valid  for  small  amplitude  acoustic  waves ) 
and  by  assuming  that  the  wave  follows  adiabatic  processes,  yields  the 
wave  equation: 


v  _  JL  ais  _ 
af  at* 


Eq.  (19) 


where 


is  a  parameter  of  the  wave  such  as  pressure  or  density 
is  the  Laplacean  operator 


is  the  speed  of  sound 
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Inclusion  of  terns  Involving  propagation  parameters  and  second  order 
quantities  can  be  grouped  In  the  operational  calc\ilus  form: 


1  P.(t)  D,(¥) 


Eq,  (20) 


where 

(t)  is  a  function  describing  the  departures  from  the  mean  of 
variables  describing  the  local  properties  of  the  medium 
such  as  temperature 

D,  is  a  linear  differential  or  integral  operator  with  respect 

to  space  or  time 

This  results  in  the  modified  wave  equation. 


^  =  3;  P.  (?)  D.  (V)  Eq.  (21) 

at  o  t  B  " 


When  applying  the  method  of  perturbations,  it  is  now  as¬ 
sumed  that  the  wave  function  after  a  single  scattering  is  composed  of 
the  incident  wave  'I'jeind  a  scattered  wave.  If  it  is  then  fxirther 
assumed  that  the  incident  wave  can  be  represented  by 


where 


A  is  the  amplitude 

•T  is  the  wave  vector 

7  is  the  radius  vector 

u  is  the  angular  frequency, 

the  sec'jnd  approocimatlon  to  the  wave  function  is  that 


Eq.  (22) 


ijf  = 


Eq.  (23) 


i*vwwTCii^iv\Mv%MW«uwuwu»fuirw\sw  ww«ar\nruww«  vtoirwvutfV'VMvuvu  ^Ji#vn#vrywuwJi^vfvwuwuvuv  WAAiV*  AK  ri#>7VJWA 


Substitxrtion  of  Eq.  (22)  and  Eq.  (23)  in  Eq,  (21 )  yields; 


y"m _ ^  ^  —  Api(K*^-“0 

‘  ^  dt^  ~ 


Q  (r) 


Eq.  (24) 


vhere  Q(T)  represents  the  results  of  performing  the  operations  indicated  hy 
Eq.  20. 


Ihe  formal  solution  of  Eq.  (24)  is 


= /^e' Q(70 


vhere  F'ls  a  small  domain  near  rand  Q  (rO  is  a  function  describing  the 
relevant  properties  of  the  medium  in  r*.  IMs  function  represents  an 
acoustic  pulse  originating  at  the  scattering  center. 


The  quantity 


^1  (k.7/-  «t)  +  i  |k  II  r-r'  I 


1 7-?' I 


Eq.  (26) 


can  be  expanded  in  an  infinite  series  to  that  at  large  distances  from 
the  scattering  element  >diere  the  initial  wave  has  not  been  affected  by 
the  redirection  of  energy,  the  second  order  terms  in  the  series  can  be 
neglected.  The  scattered  wave  function  is  then 


_  A  pi  (KT  —  »t) 

Q(T')d®7'  Eq.  (27) 


Tho  quantity  within  the  integral  is  the  Fourier  transform 
of  the  function  tliat  describes  the  distribution  of  a  relevant  property 
of  the  medium  such  as  temperature  fluctuations  or  wind  within  the  vol¬ 
ume  >diich  is  projected  along  the  scattering  vector  k  by  the  "dot"  product. 


A  pKi^’T-at)  _ 


Eq.  (28) 
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Equation (28)  Inc Icates  that,  within  the  small  scattering  region  Q(P),  the 
only  part  of  the  complex  amplitudes  of  the  Fourier  representation  of  the 
turbulence  which  reacts  with  the  incident  wave  so  that  phases  combined 
is  the  component  along  the  scattering  vector  It  (see  Fig.  31 )•  All  other 
amplitudes  will  combine  out  of  phase  to  cancel  each  other. 

The  directional  disti'ibution  of  intensity  of  the  scattered 
energy  can  be  formulated  in  terms_of  the  general  scattering  cross-sec¬ 
tion,  0  (t)  •In  physical  terms,  a(i)ls  defined  as  the  ratio  between  the 
intensity  of  scattered  energy  per  unit  solid  angle  in  the  direction  J  per 
unit  scattering  volume  and  the  Intensity  of  the  incident  wave  per  unit 
area  of  the  wave  front.  The  intensity  of  any  ira.ve  may  be  found  by  taking 
the  mean  of  the  product  of  the  wave  function  V  and  its  complex  conj^ate 

*or  Thus,  the  scattering  cross-section  is  proportional  to  the 

ratio  of  the  intensities  of  the  scattered  wave  to  the  incident  wave,  i.e., 


o(i) 


oc 


UMW* 

^00 


Further,  it  can  be  shewn  that  the  spectrum  of  a  turbulent  quantity  4>(k) 
is  given  by  the  mean  of  the  Fourier  representation  of  the  turbulence  3*  and 

its  complex  conjugate  or  Eq. (28) expresses  the  relation 

between  the  scattered  wave  function  Wj  and  JF  so  that  'Pj'Pj*  in  terms  of 

?y*nay  be  used  to  give  the  scattering  cross-section, 


®  (t)  —  ®(k) 


Eq.  (29) 


where  4>(k)is  the  spectrum  of  the  scattering  quantity  such  as  wind  or 
temperature. 

The  scattering  cross-section  for  txirbulent  temperature 
fluctuations  can  be  found  by  finding  the  Foiirier  transform  of  the  tur¬ 
bulent  temperature  field  and  applying  Eq. (29).  This  gives 


=  COS^  e  «I)^(k)  Eq.  (30) 

for  the  directional  distribution  of  intensity  of  energy  scattered  by 
turbulent  temperature  field.  The  T  subscript  refers  to  temperature.  A 
sJ.mllar  analysis  for  scattering  by  a  turbulent  wind  field,  U  ,  gives 
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for  the  angular  distribution  of  intensity,  where  E  is  the  turbulent 
energy  spectrum  function  along  the  incident  wave  vector  Tf . 


When  Eq, (31) is  integrated  over  a  sphere  surrounding  the 
scattering  region  to  find  the  total  amount  of  energy  scattered  from  u 


weak  shock  wave,  the  integral  is  infinite  (l.e.  since  cot-^  approaches 


Infinity  as  $  approaches  zero,  an  infinite  amount  of  energy  Is  directed 
forward).  In  order  to  resolve  this  discrepancy,  Llghthlll  (Ref.  56) 
treated  the  scattering  of  h  shock  wave  by  turbulence  as  being  due  to  the 
Interaction  of  an  acoustic  inilse,  represented  bj'  a  Fourier  sum,  with  a 
set  of  turbulent  quadrupoles  with  a  velocity  discontinuity  (ea,,)  across 
the  shock.  He  then  assmes  that  the  scattered  energy  propagates  with 
the  speed  (cao)nnd  develops  the  scattering  equation 


«U)  = 


64nr2 


L  2 


ico8  e 

(sin ^4 


sin^  $ 


Eq.  (32) 


where 


,  =  i  - 1 

M, 


M,=  shock  Mach  number. 


This  equation  predicts  a  maximum  of  scattered  energy  in  the 
direction  in  which  the  original  wave  is  propagating  after  a  single  scat¬ 
tering,  and  that  the  intensity  of  tlie  scattered  energy  will  obey  the  in¬ 
verse  square  lav.  Thus  the  scattered  energy  will  remain  finite.  Ihe 
intensity  of  the  scattered  energy  depends  upon  the  square  of  the  wind 
speed  fluctuations. 


lo  has  been  Impllctly  assumed  in  all  studies  of  the  scat- 
terlt^  of  sound  by  turbulence  that  the  interaction  between  the  acoustic 
energy  and  the  ‘■.urbulence  is  weak.  However,  Meecham  and  Ford  (Ref.  59) 
have  shown  that  the  power  spectra  of  the  sound  emitted  by  the  turbiilent 
8\irface  layer  have  two  fozms  depending  on  whether  the  emission  Is  to  the 
left  or  right  of  the  peak  in  the  inertial  subrange  as  defined  in  Fig.  30. 


tlon  given  by 


In  the  inertial  subrange  $  (oj)  is  proportional  to  a  func- 


21 
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in  db 


Eq.  (33) 
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In  the  raxige  of  the  larger  eddies  the  relation  la. 


♦  (•)- 9VaJ  M®  (■^)^  in  db 


Eq.  (3*^) 


vhere  4(tt)  Is  the  spectral  density,  p  Is  the  density,  V  Is  the  volioBe 
being  considered,  ho  Is  the  speed  of  sound,  M  Is  the  Mach  number  of  the 
turbulence  defined  as  the  IOC  wind  speed  divided  by  the  speed  of  sound, 
o  Is  the  angular  frequency  of  emission  cmd  i  Is  a  characteristic 
"length"  of  the  turbulence.  A  schoiatlc  expressing  this  relation  Is 
ehovn  In  Fig.  32. 

Since  most  of  the  turbulent  acoustic  pover  Is  emitted  at 


It  seems  reeisonable  to  assume  that  the  strongest  Interaction  betveen  the 
shock  vave  end  turbulence  will  occur  with  this  particular  size  eddy. 

Oils  ml£^t  imply  excitation  of  the  eddy  to  Wgb  energy  levels  and  sub¬ 
sequent  re-emlsslon  of  a  decaying  pulse  of  energy  or,  resonant  scatter¬ 
ing.  .Alternatively,  this  Interaction  ml^t  Introduce  a  perturbation  on 
the  shock  vave  front  (Ref.  60,  p.  Il4).  This  vave  vould  then  travel 
aloxig  the  shock  front  transversely  to  the  dlrectlm  of  propagation  of 
the  shock.  !Ilhls  effect  vould  be  In  analogy  of  polarization  of  the  tur¬ 
bulent  elements  postulated  by  Ll^thlU  (Ref.  56). 


FREQUENCY  («) 


Fig.  32  Schematic  of  Turbulonf-Acouttlc  Radiation, 
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(4)  Criti^M*  of  Curront  Acowotic  Scottoriof  Tkoory  o>  Applied  to  Shock  Wovot  - 

It  Is  assumed  in  current  scattering  theory,  either  explicitly  or  impli¬ 
citly,  that  the  linearized  wave  equation  holds,  that  the  Interaction  be¬ 
tween  the  turbulence  eind  the  acoustic  energy  is  weak,  and  that  the  me¬ 
thod  of  small  perturbations  is  valid.  Further,  three  dimensional  geo¬ 
metry  (implying  homogeneous  turbulence)  is  assumed  which  (Ref.  57)  oay 
not  be  applicable,  particularly  in  view  of  the  observed  tendency  of 
turbulent  elements  (Ref.  4l,  p.  210)  to  elongate  along  the  wind.  In  two 
dimensions,  Huygshens*  principle  does  not  hold,  and  a  pulse  of  radiation 
will  persist  as  a  reverberation.  In  addition,  it  is  assvmed  that  the 
method  of  peirturbatlons  may  be  used  (in  some  cases  twice)  and  that  the 
oscillations  of  the  scattering  elements  occur  in  a  particular  mode. 

The  degree  of  caiq)lexlty  in  studying  scattering  depends 
on  the  strerigth  of  the  Interaction  of  the  incident  wave  with  the  scat¬ 
tering  elements.  If  the  distance  between  turbulent  elements  is  large, 
it  may  be  assumed  that  the  scattered  waves  from  turbulent  elements  do 
not  Interact  with  each  other.  If  on  the  other  hand  there  is  a  hl^ 
density  of  scattering  elements  the  waves  Interact,  and  any  classification 
into  weak  and  multiple  scattered  waves  loses  its  meaning.  If  the  scat¬ 
tering  is  weak  there  are  two'  types  of  scattering,  coherent  and  inco¬ 
herent  (Ref.  58).  Coherent  scattering  occurs  when  the  scattered  waves 
are  in  phase  in  some  regions  and  strong  reinforcement  occurs.  Inco¬ 
herent  scattering  occurs  when  all  the  waves  are  generally  out  of  phase 
and  destructive  interference  occurs.  The  degree  of  coherence  depends 
upon  the  distribution  of  scatterers  and  the  orientation  of  the  scattering 
wave  vector  Ic  .  If  the  scatterers  can  'be  treated  as  'being  orderly  dis¬ 
tributed  in  horizontal  planes,  and  if  the  scattering  vector  Tc  is  nor¬ 
mal  to  this  horizontal  plane,  the  maximum  interaction  occurs  and  strong 
beams  are  -  produced  at  intervals  of  (ml)  where  m  is  an  Integer  and  X 
is  the  wave  length. 

It  is  fairly  conclusive  from  the  N-wave  data  obtained  at 
Oklahcma  City  that  the  maximum  of  scattered  energy  is  radiated  at  angles 
to  the  direction  of  propagation  other  than  straight  ahead  as  is  pre¬ 
dicted  by  the  theory  which  results  in  Eq. (27)for  high  frequerny  fluctua¬ 
tions  and  Eq. (26) for  shock  wave  scattering.  Ihe  reason  for  this  conclu¬ 
sion  is  derived  from  the  statistical  analysis  of  the  Oklahoma  City  sonic 
boom  test  data  presented  in  Section  IV. B. 4.  IHiis  data  shows  that  Test 
House  4,  at  which  point  the  shock  wave  -was  at  angle  of  approximately 
70  degrees  to  the  ground,  had  consistently  more  spiked  N-v^-ves  of  higher 
relative  amplitude  than  Test  House  1  where  the  shock  wave  was  at  an 
angle  of  about  45  degrees.  If  it  is  assianed  that  the  intensity  of  the 
scattered  wave  falls  off  as  some  inverse  function  of  distance  from  the 
scattering  center  (such  as  inverse  square  law  for  spherical  wa-ves),  the 
observations  mentioned  above  Implies  that  the  scattering  vector  Tc  ,  de¬ 
fined  in  Fig.  31>  is  more  perpendicular  to  the  ground  at  Test  House  4  than 
at  Test  House  1.  This  observation  does  not  agree  with  Eq.  (28)vdiich  pre¬ 
dicts  a  mEixlmum  scattering  in  the  forward  direction.  Inspection  of  N- 
wave  traces  shows  that  the  "spikes"  or  peaks,  quite  often,  nearly  coin¬ 
cide  with  the  leading  edge  of  the  trace.  This  may  indicate  that  the  ef¬ 
fect  of  the  turbulence  introduces  a  perturbation  on  the  shock  front  which 


travels  transversely  along  the  shock  wave  with  a  "group"  velocity  which 
Is  fester  than  the  Mach  number  of  the  shock  wave.  IHie  osclllatoiy  na¬ 
ture  of  the  traces  further  Indicates  that  seme  degree  of  resonance  is 
present.  Some  possible  mean  radiation  patterns  for  a  resonating  scat- 
tcrer  are  shown  in  Fig.  33* 
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Fig.  33  MultIpoU  Radiation  Pattornt, 


An  analytical  ccoparlson  of  the  relative  importance  of  the 
scattering  effects  of  a  turbulent  velocity  field,  the  scattering  of  a 
turbulent  temperature  field,  and  the  scattering  when  both  fields  ai% 
present  has  not  been  made  at  this  time*  It  is  very  likely  that  the  re¬ 
sults  of  the  obsesrved  data  which  are  not  consistent  with  the  theory  are 
a  result  of  the  interaction  of  the  shock  wave  with  a  medium  in  which 
both  the  temperature  and  wind  fields  are  turbulent.  Ibe  characteristics 
of  the  turbulent  atisospbere  dasorlbed  In  Beesion  ZI.C.l  indicates  that 
further  progress  in  quantitatively  describing  the  effects  of  the  earth's 
turbulent  boundary  layer  on  the  sonic  boom  vill  depend  upon  theoretically 
extending  scattering  theory  to  incozrorate  combined  turbulent  tempera¬ 
ture  wind  fields,  a  study  of  the  effects  of  inhomogeneity  and  non-iso¬ 
tropy  in  the  turbplenctf  flow,  and  concomitant  experimental  measurements 
to  seiTTe  as  a  guide  to  the  theo»^tlcal  development.  One  of  the  princi¬ 
pal  experimental  problems  is  to  determine  the  extent  of  the  regions  in¬ 
volved  in  the  single  scattering  process  together  with  the  altitude 'at 
which  it  occurs.  If  this  region  is  small  and  the  altitudes  low,  the 
homogeneous,  isotropic  assumption  will  be  a  good  approximation. 

Work  to  date  (see  Appendix  IX )  has  shown  that  the  log¬ 
normal  statistical  distribution  will  probably  describe  the  areal  distri¬ 
bution  of  overpressures.  Also,  the  pertinent'  parameters  are  probably 
the  direction  of  propagation  of  the  N-wave  and,  when  taken  separately, 
the  spectrum  of  wind  or- the  spectrum  of  temperature  fluctuations.  Some 
progress  has  been  made  in  merging  the  effects  of  the  two  fields,  but  it 
is  not  sufficiently  complete  to  be  Incorporated  in  this  report. 

(D)  MISCELLANEOUS  PHENOMENA  -  Various  considerations  have  been  men¬ 
tioned  in  the  previous  material  which  could  not  easily  be  reviewed  in 
any  of  those  sections.  Ihe  purpose  of  this  section  is  to  discuss  these 
considerations  in  detail  so  that  they  may  be  used  in  evaluating  the  sonic 
boom  produced  on  the  ground  under  various  circumstances. 

(1)  RtfUctlon  Factor  Naor  Cvt-off  -  In  Sections  III. A  and  III.B  meteo- 
rolog.tcal  conditions  were  defined  which  could  cause  local  Intensification 
or  focusing  of  the 'shock  strength.  It  was  noted  in  the  development  of 
these  criteria  that  focusing,  if  it  occurred,  accompanied  cut-off.  At 
cut-off,  the  shock  front  is  locally  normal  to  the  ground.  This  fact 
raises  a  question  about  the  -value  of  the  reflection  factor,  K|{,  which 
should  be  used  in  calculating  the  overpressure. 

Consider  first  a  weak  oblique  shock  front  in  free  air  as 
shown  in  Fig.  3^(b)« 
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Fif,  34  Shock  Front  Conflgurotlonn. 


The  pressure  Jump  across  this  shock  wave  In  free  air  Is 
given  by  P 2 -Pi  and  the  streamline  describing  the  flow  Is  turned  through 
the  angle  a  •  The  reflection  factor, Kr,  Is  defined  as  the  total  pres¬ 
sure  Jump  divided  by  the  pressure  Jump  In  free  air.  Thus,  for  the  shock 
configuration  shown  In  Fig.  3**(a): 
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If  a  boundary^  such  as  the  ground.  Is  Inserted  then  the 
flow  must  remain  parallel  to  that  boundary.  This  can  be  accomplished  by 
the  presence  of  a  second  (reflected)  shock  wave  behind  the  first  (inci¬ 
dent)  wave  as  shown  In  Fig.  3**(b).  Now,  however,  the  total  pressure 
Jump  across  this  configuration  ISP3-P1.  If  the  shock  waves  are  weak 
It  can  be  shown  (Ref.  7I-)  thatP2-Pi  “P3-P2  bo  that  In  this  case: 
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However,  for  the  normal  shock  wave  shown  in  Fig.  3^(c),  no 
reflected  wave  is  required  to  turn  the  flow  parallel  to  the  ground  so 
that  the  total  pressure  Jump  across  this  wave  is  just  P2-Pj  i-e., .the 
free  air  value.  Thus,  in  the  case  of  the  normal  shock  wave: 


It  is  obvious  that  Kj^  must  vary  between  2,0  for  the  weak 
oblique  shdck  front  and  1.0  for  the  normal  shock  front.  This  is  to  say 
that  the  reflection  x'actor  is  a  function  of  the  shock  angle.  Exactly 
how  Kjj  varies  in  net  fully  understood  at  this  tiroe>  however,  some  quali¬ 
tative  insight  may  be  obtained  from  Fig.  3^(c).  If  the  incident  shock 
angle  is  very  near  90“  and  the  shock  wave  is  weak  the  reflected  wave  can¬ 
not  turn  the  flow  parallel  to  the  ground.  This  situation  requires  that 
a  "triple  point"  exist  (Refs.  72  and  73)>  such  as  that  shown,  where  the 
shock  near  the  ground  is  normal  and  above  it  at  some  point  the  waves 
split  into  an  incident  and  a  reflected  wave.  This  situation  probably 
represents  the  transition  between  the  oblique  and  normal  shock  fronts. 
Some  general  aspects  c2  this  problem  are  considered  further  in  Refs.  Jk 
and  82. 

From  this  discussion  it  would  appear  thatKji  should  be  con¬ 
sidered  a  variable  near  cut-off.  Assuming  that  this  is  the  situation 
the  variation  of  the  overpressure  near  cut-off  and  focusing  would  take 
on  the  form  shown  in  Fig.  35* 

The  figure  shows  that  an  instrument  mounted  on  the  ground 
would  record  some  increase  in  the  overpressure  near  cut-off  but  it  would 
not  record  as  high  an  overpressure  as  predicted  using  Kr=  const.  =  2.0. 

It  would  be  desirable  to  determine  theoretically  the  varia¬ 
tion  of  Kr  with  shock  angle.  Some  of  the  preliminary  work  has  been  done 
in  Refs.  72  and  7^  but  this  work  would  have  to  be  expanded  to  obtain 
a  proper  solution.  It  would  also  be  of  considerable  interest  to  experi¬ 
mentally  confirm  the  theoretical  analysis,  both  concerning  the  variation 
of  Kp  and  the  variation  of  the  free  air  overpressure  predictions  near 
cut-off  and  focusing.  Some  possible  methods  of  accomplishing  this  are 
discussed  in  Section  V. 

(2)  High  Altitude  Turbulence  and  Shower  Clouds  -  The  problem,  that  of 
defoimation  of  the  N-wave  by  high  altitude  turbulence,  is  essentially 
that  of  finding  the  scattered  field  at  a  great  distance.  This  scattering 
problem  is  probably  one  in  which  single  sca'i'tering  occurs  since  most  high 
altitude  turbulence  layers  are  thin  (on  the  order  of  one  to  not  much  more 
than  3^000  feet  thick).  This  will  probably  result  in  deformation  of  the 
ll-wave  to  some  degree,  but  the  extent  of  this  deformation  and  the  effect 
of  propagation  of  a  deformed  N-wave  over  great  distances  is  \mknown  at 
this  time. 


MACH  NUMBER 

Fig.  35  Sonic  Boom  Noor  Cut-Off. 

The  problem  finding  the  effect  of  shower  clouds  on  the 
sonic  boom  is  much  greater.  Data  taken  during  the  Oklahoma  City  tests, 
such  as  May  29  at  O70O  and  O9OO  CST  and  May  31  at  0700  CST,  indicate 
that  profound  modi fi cat i -.ns  of  the  N-wave  can  occur  as  a  result  of  these 
clouds.  Cumulus  type  clouds  (such  as  thunderstorms)  are  essentially 
vertical  in  character  and  have  a  relatively  limited  horizontal  extent. 
The  boundaries  between  upward  and  downward  air  motions  and  between  re¬ 
gions  of  precipitation  and  no  precipitation  indicate  that  study  of  shock 
propagation  through  these  clouds  will  present  a  three  dimensional  prob¬ 
lem  in  multipath  propagation,  edge  diffraction  effects  as  well  as  the 
problems  associated  with  turbulent  modification  of  the  N-wave  within  the 
cloud  and  in  the  lower  levels.  No  significant  effort  has  yet  been  di¬ 
rected  towards  finding  a  solution  to  this  problem. 

(E)  SUMMARY  OF  RESULTS  —  The  meteorological  conditions  required 
to  produce  anomalous  propagation  such  as  focusing,  complete  cut-off,  ex¬ 
treme  lateral  spread,  and  deformation  of  the  pressure  wave  signature 
were  discussed  in  this  section.  Criteria  were  established  to  predict 
the  presence  of  these  conditions.  It  was  shown  that  realistic  varia¬ 
tions  in  both  temperature  and  wind  could  produce  focusing  or  complete 
cut-off  for  flight  at  Mach  nxaabers  below  1,3,  However,  the  focusing 
would  be  accompanied  by  cut-off  which  would  preclude  the  normal  doubling 
of  overpressure  at  the  ground  due  to  reflection  of  the  shock  wave.  Ex¬ 
treme  lateral  spread  and  distortion  of  the  pressure  wave  signature  due 
to  interactions  with  turbulence  could  occur  for  flight  at  all  Mach  num¬ 
bers.  However,  the  foimer  would  not  occur  for  flight  at  altitudes  above 
those  ■vdiere  the  maximum  wind  speeds  exist  regardless  of  the  Mach  nimber. 
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SECTION  IV  ANALYSIS  OF  OKLAHOMA  CITY  DATA 

A  limited,  amount  of  the  flight  test  data  measured  during  the  sonic 
boom  tests  at  Qhlahoma  City  has  been  analyzed.  ^Thls  analysis  has  given 
some  Insight  Into  the  problems  of  evaluating  the  data  when  comparing  It 
to  the  theoretical  results.  The  purpose  of  this  section  Is  to  present 
these  comparisons  between  theory  and  test  data.  In  general,  excellent 
agreement  was  observed  In  the  cases  considered. 

(A)  INFLUENCE  OF  TEST  AIRPLANE  GEOMETRY  -  Each  of  the  test  alrj^es 
was  considered  In  detail  and  the  shock  strength  parameter 
(see  Eq.  1)  was  determined.  The  theory  of  Bef.  was  used  In  develop¬ 
ing  the  shock  strength  parameter  for  several  angles  of  i.  The  results 
for  each  airplane  eure  briefly  reviewed  In  this  section. 

(1)  F-104A  -  Geometric  descriptions  for  the  F-104A  were  ob¬ 

tained  from  the  Iiockheed-Callfomla  Company.  These  data  were  analysed 
using  the  theory  of  Ref.  73*  It  was  found  that  for  this  particular  air¬ 
plane  more  than  two  shock  waves  would  be  produced  at  some  lift  coeffi¬ 
cients.  Special  care  was  required  Ih  determining  the  proper  value  of 
[KY.,#)]  for  the  front  shock  wave.  The  results  of  the  analysis 
are  shown  In  Fig.  36  for  several  angles  Of  S ,  IQie  figure  shows  Idiat  for 
a  rsuige  of  lift  coeff Ic*  ents,  the  shock  strength  parameter  Is  Inde¬ 

pendent  of  the  angle  6- 


I 

Fif.  36  Shock  Sfrongth  Paramotor  for  F-104A  Airflono. 
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(2)  B-58A  -  Gecmetric  descriptions  for  the  B-58a  airplane  were 

obtained  frcm  the  Fort  Worth  Division  of  General  Dynamics.  These  data 
were  analyzed  using  the  theory  of  Ref.  75  for  several  angles  of  f  with 
and  without  the  MB  pod.  Results  of  this  analysis  are  shown  in  Fig.  37* 
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W  =  AIRPLANE  WEIGHT 
S  =  AIRPLANE  WING  AREA 
Pa  =  AMBIENT  PRESSURE  AT  AIRPLANE 
M  =  AIRPLANE  MACH  NUMBER 


LIFT  PARAMETER  - 
Fig.  37  Shock  Strength  Parameter  for  B—58A  Airplan 


(3)  F-101B  -  Geosetrlc  data  for  the  P-IOIB  were  obtained  fx'om  the 
McDonnel  Aircraft  Corporation.  Hiese  data  were  eumlyzed  using  the  theory 
of  Ref.  75  for  several  angles  of  $.  Results  are  shown  In  Fig.  38. 

(B)  COMPARISON  OF  THEORY  AND  TEST  -  A  limited  number  of  comparisons 
have  beep  made  between  test  results  eind  theoretical  predictions  using 
the  data  obtained  during  the  Oklahoma  City  sonic  boom  test  series.  These 
ocmparisons  have  been  generally  broken  into  a  atxady  of  two  areas  of  the 
problem.  The  first  concerns  predictions  assuming  horizontally  stratified 
atmospheric  models.  This  has  led  to  detailed  comparisons  of  normal  pres¬ 
sure  wave  traces  with  those  predicted  by  the  theory.  The  second  concerns 
an  Investigation  of  deformed  pressure  wave  signatures.  These  have  been 
studied  both  statistically  and  In  detail. 

nie  test  data  used  for  these  computations  and  comparisons  con¬ 
sisted  of  pressure-time  traces  recorded  by  standard  NASA  Instrumenta¬ 
tion  (Ref.  84),  upper  air  winds,  temperat\ires,  and  pressures  measured 
every  two  hours  by  (M)-1A  eqvilpment,  euid  continuous  low  level  winds  and 
tes^ratures  obtained  by  Beckman-Whitley  probes  and  wire  sondes.  Ulie 
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GMD-IA  data  vere  recomputed  ahd  Interpolated  in  time  (Ref.  J6)  to  pro¬ 
vide  appropriate  data  for  computing  the  sonic  boom  strength  and  slgna* 
ture  at  the  ground  (Ref.  8)  for  each  test  case  considered.  Pressure 
signature  measurements  taken  in  the  vicinity  of  the  low  level  micro- 
meteorological  instrumentation  were  not  available,  so  that  none  of  these 
data  were  reduced  to  provide  a  description  of  the  tvurbulence.  Si^ch 
sonic  boom  measurements  covtld  be  used  in  a  study  of  the  interaction  be¬ 
tween  the  sonic  boom  and  lov-level  turbulence  using  the  techniques  of 
spectral  analysis. 


(1)  Nermoi  N-Wavt  Signoturti  -  ^Hie  detailed  pressure  T«i.ve  signature  n 

we.s  computed,  for  the  F-10Ua,  vising  the  theory  of  Ref.  75*  This  slgna-  5 

ture  was  interesting  in  that  an  intermediate  shock  wave  was  evident  be-  I 

tween  the  front  and  rear  waves.  The  pressures  were  corrected  for  the  j 

variable  atmospheric  properties  between  the  airplane  and  the  ground  by  I 

the  method  of  Refs.  6  through  8.  A  reflection  factor  of  2.0  at  the 
ground  was  assvaned.  The  predicted  pressure  wave  trace  was  then  compared  j 

to  a  number  of  imdeformed  measured  signatures  for  the  corresponding  alti¬ 
tude  and  Mach  number  of  the  test.  Some  results  for  the  F-lO^i-A  are  shown 
in  Fig.  39.  I 

In  general  the  figure  shows  excellent  agreement  with  the 
predicted  signature.  However,  it  is  interesting  to  note  that  the  slope  I! 

and  the  peak  pressure  of  the  measured  front  shock  wave  are  not  in  agree-  | 

ment  with  the  theory.  This  would  suggest  some  Influence  of  the  measur-  < 

Ing  system  Inertia.  That  is,  the  measured  wave  suggests  that  the  system  ! 

is  unable  to  react  Instantaneously  to  the  sharp  prescare  rise  at  the 
front  shock  wave.  The  response  time  lag  seems  to  result  in  an  observea 
peak  overpressure  which  is  from  7  percent  to  15  percent  lover  than  that  | 

predicted.  The  fact  that  the  remainder  of  the  measured  and  predicted  a 

pressure  wave  are  quite  close  suggests  that  the  assumptions  used  in  the  | 

theory  (including  K,j  ■  2.0)  are  correct.  S 

A  number  of  undeformed  pressure  signatures  produced  by  the  . 

F-IOIB  have  been  analyzed  in  the  same  manner  as  descrj'.bed  above.  Some  f 

typical  results  are  shown  in  Fig.  ko.  | 

These  results  are  similar  to  the  F-104a  results  in  that  r 

the  agreement  is  quite  close.  However,  the  pressure  rise  at  the  shock  ,l 

waves  is  not  Instantaneous  as  predicted  by  the  theory.  In  these  cases  , 

the  theoretical  overpressure  at  the  bow  shock  wave  is  from  5  percent  to  I 

10  percent  higher  than  the  measured  value.  This,  again,  seems  to  be  due  | 

to  the  observed  rise  rate  at  the  front  shock.  In  all  other  respects  the  | 

theory  seems  to  have  adequately  predicted  the  pressvire  signature  observed  | 

on  the  ground.  ii 

A  yeiy  limited  number  of  undefonned  B-58A  pressure  wave  ' 

signatures  were  available  for  analysis.  These  were  analyzed  In  the 
manner  described  above  and  the  results  are  shown  in  Fig.  4l. 
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Fig,  39  Comparison  of  F-104A  Thoorotieal  and  Monsvrod  Prossuro  Wavo  Signaturos, 
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The  figure  shows  that  the  agreement  of  the  obsei’ved  pres¬ 
sure  trace  with  the  theory  is  very  .close.  Sufficient  data  were  not 
available,  however,  to  estimate  the  differences  between  the  theoretical 
and  measured  front  shock  overpressure  which  resultj  the  slow  rise  of 
recorded  signature. 

Ihese  comparisons  indicate  that  when  considering  only  the 
agreement  between  the  predicted  and  observed  overpressure  of  the  front 
shock  wave  for  undeformed  signatures,  a  reflection  factor  of  2.0  may 
yield  results  which  are  higher  than  the  observed  value.  This  may  be 
compounded  by  the  fact  that  measured  data  repeatability  for  the  i  ont 
shock  overpressure  may  vary  by  as  much  as  ±12  percent.  Ihese  considera¬ 
tions  have  led  to  the  conclusion  that  when  comparing  experiment  and 
theory  the  whole  pressure  signature  should  be  analyzed.  Comparisons  of 
front  shock  overpressure  alone  could  easily  lead  to  erroneous  conclusions 
concerning  the  validity  of  the  theory. 

(2)  Distorted  H-Wovo  Si^iioturos  -  Comparisons  between  predicted  and 
measured  F-104a,  F-IOIB,  and  B-58A  pressure  signatures  were  developed 
assuming  horizontally  stratified  atmospheric  models,  =  2.0,  and  using 
the  theory  of  Ref.  75*  Comparisons  of  measured  spiked  waves  and  mea¬ 
sured  rounded  waves  for  the  F-104A  and  F-lOlB  are  shown  in  Fig.  1+2  and 
Fig.  1+3  respectively.  A  similar  comparison  of  a  spiked  wave  produced 
by  the  E-58A  is  shown  in  Fig.  1+1+.  (No  rounded  waves  were  obtained  in  che 
limited  drta  observed  for  this  airplane  during  the  Oklahoma  City  tests. ) 

Again,  these  figures  seem  to  illustrate  seme  effect  of  in¬ 
strument  response  in  that  the  pressxxre  rise  of  the  front  shock  is  not  in- 
steintaneous,  as  predicted  by  theory.  This  is  probably  caused  by  the  in¬ 
ability  of  the  measuring  system  to  respond  instantaneously  to  the  pres¬ 
sure  rise  across  the  front  shock.  It  is  interesting  to  note  that  in  the 
ca^e  of  the  spiked  waves  the  pressure  trace  is  generally  similar  to  the 
predicted  normal  signatures  with  the  exception  of  the  Impulses  attached 
to  each  shock  wave.  The  forward  portion  of  the  ro+inded  wave  is  similar 
to  the  predicted  wave  but  the  peak  is  severely  flattened.  It  would 
appear  that,  at  one  time,  both  signatures  were  of  the  type  predicted  by 
theory.  Somevdiere  in  the  atmosphere  between  the  airplane  and  the  in¬ 
strument  they  seem  to  have  been  distorted  by  the  addition  or  subtraction 
of  a  press+xre  pulse  at  each  shock  wave.  A  feature  of  the  distorted  waves 
is  that  the  portion  of  the  signature  away  from  the  shocks  agrees  reason¬ 
ably  well  with  the  theoretical  (and  observed)  undlstorted  signatures. 

In  genereuL,  the  maximum  overpressure  of  the  measured  spiked 
front  wave  exceeds  tliat  predicted  by  theory  while  the  overpressure  of 
the  rounded  front  wave  is  less  than  that  predicted  by  theory.  Because 
of  this  type  of  distortion  of  the  measured  pressure  wave  signatures, 
extreme  care  should  be  exercised  when  comparing  only  froit  shock  wave 
overpressure  data  with  the  predicted  theoretical  values.  If  this  sort 
of  comparison  is  ' t  be  made,  each  wave  should  be  classified  as  to  its 
form,  i.e.  spike  rounded,  normal,  or  a  combination  of  these.  The  most 
significant  comparisoi.  which  could  be  made  are  of  the  type  shown  in 
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Figs.  39  throu^  44,  Effort  of  this  sort  would  yield  a  much  better  un¬ 
derstanding  and  proper  Interpretation  cf  the  data. 


The  data  In  Fig,  4^  shows  how  atmospheric  turbulence  may 
affect  the  shape  of  the  sonic  boom  N-wave  signature  over  short  distances 
along  the  ground.  The  measurements  were  taken  at  O73O  hours,  in  a  flat, 
open  area  near  Lalre  Hefner  (northwest  of  Oklahcma  City).  The  overpressure 
signatiures  were  produced  by  an  F-104a  aircraft,  flying  at  Mach  1,7,  at 
an  altitude  of  28,000  ft.  The  weather  conditions  at  the  ground  during 
the  time  of  flight  were:  clea  sky,  visibility  greater  than  15  miles, 
barometric  pressure  28. l4  inches  mercury,  temperature  65  degrees  Fahren¬ 
heit,  relative  hianidlty  43  percent,  and  surface  wind  from  210‘  at  42,6 
ft/sec.  The  radio  sonde  ds,ta  at  O7OO  hours  indicate  a  stable  temperature 
lapse  rate,  but  the  strong  wind  and  the  65®F  temperatui-e  at  O73O  suggests 
that  a  strong  tiirbulent  layer  about  500  feet  deep  had  de-veloped.  The 
siui  was  12  degrees  above  the  horizon,  \diich  implies  that  the  groimd 
heating  and  consequently  the  vertical  wind  velocities  were  negligible. 

As  a  result  of  these  considerations,  the  low  level  tui’bulence  can  be  as¬ 
sumed  to  be  two-dlii!'*nsionally  isotropic  and  homogeneous  in  the  horizontal 
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(3)  Lotaral  Distribution  of  Sonic  Boom  -  Lateral  dlstrltutlon  of  the 
front  shock  wave  strength  was  investigated  for  several  flights  durii.g 
the  Oklahoma  City  tests.  Theoretical  predictions  were  obtained  using 
the  method  of  Refs.  6  through  8.  Data  for  each  fli^t  were  obtained 
from  several  sources.  Uie  airplane  shock  strength  parameter  was  ob¬ 
tained  from  Fig.  36  for  the  F-lOhA  airplane.  Meteorological  data 
throu^out  the  day  were  obtained  from  rawlnsonde  measurements  of  the 
Air  Force  Mobile  Weather  Squadron  (Tinker  AFB).  A  digital  computer 
program  (Ref.  76)  was  designed  to  compute  the  appropriate  data  for  the 
time  of  each  fll^t  from  the  rawinsonde  measurements.  This  program 
was  used  to  define  the  horizontally  stratified  model  atmospheres  for 
each  fli^t.  A  nuLber  of  comparisons  of  the  type  shown  in  Fig.  39 
were  used  to  define  the  ratio  between  the  measui'^d  and  theoretical  bow 
shock  wave  overpressure.  This  was  done  to  allow  proper  evaluation  of 
the  measorements  in  the  li^t  of  the  instrument  system  response.  Frcm 
these  comparisons  an  instrumentation  response  factor  of  O.9  vas  selec¬ 
ted  for  evaluation  of  the  F-104a  front  shock  wa->  j  overrressures  only. 
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IHieoretlcal  prediction  of  the  lateral  distrihution  of  the 
front  shock  wave  overpressures  was  calculated  for  several  flints  with 
the  above  assumptions  and  appropriate  calculated  data.  Ih!  measured 
data  was  classified  into  several  groups  of  signature  types  which  in¬ 
cluded  spiked  waves  (noted  hy  an  S),  rounded  waves  (R)^  normal  waves  (N), 
and  comhinatlons  of  these  such  as  normal  waves  with  sll^t  spikes  (NS), 
normal  waves  with  slight  rounding  of  the  peaks  (NR),  and  rounded  waves 
with  spikes  (RS).  These  data  were  compared  to  the  theoretical  predic¬ 
tions  for  the  undefonned  signature.  A  typical  set  of  these  comparisons 
Is  shown  in  Fig,  ^<^6  for  the  flints  of  April  19,  19Sk. 

These  data  show  that  with  the  exception  of  the  first  flight 
(0700)  all  the  normed  wave  traces  agree  to  within  about  ±10  percent  of 
the  predicted  vedue.  This  variation  is  within  the  experienced  how  wave 
overpressure  measuijement  error.  It  may  also  he  noted  that,  in  general, 
the  spiked  front  wave  overpressures  are  higher  theuri  that  predicted  for 
the  normal  wave,  while  the  rounded  front  shock  overpressures  are  less 
than  predicted.  This  Illustrates  that  when  the  theory  is  correctly 
evaluated  agreement  with  measured  data  is  quite  close. 

(4)  Statittical  Analytit  of  N-Wavo  Amplitudot  -  During  the  passage  of  the 
sonic  hoom  throu^  the  lower  turbulent  levels  of  the  atmosphere  its  in¬ 
ters  ‘tion  with  this  turbulence  produces  marked  effects  on  its  shape  and 
i  amplxtades.  There  is  evidence  that  turbulent  scattering  is  actively 

I  effecting  the  sonic  boom  signature.  This  evidence  may  be  deducted  from 

i  almost  any  measured  N-wave  by  the  following  reasoning;  If  there  is  tur¬ 

bulent  scattering  of  the  shock  wave,  the  changes  of  the  turbident  struc¬ 
ture  of  the  atmosphere  between  the  Instant  of  passage  between  the  bow 
shock  wave  and  the  aft  shock  wave  are  very  small  sund  the  front  and  aft 
shock  waves  should  be  effected  in  the  same  way.  Hie  following  experiment 
can  be  used  to  test  this  deduction.  First,  make  a  tracing  of  the  distorted 
N-wave  signature.  Then,  place  this  tracing  below  the  original  trace  so 
i  that  the  axis  corresponding  to  time  is  parallel  to  the  original  N-wave 

signature  (as  in  Fig.  47).  Next,  translate  the  trace  (to  tho  left  in 
this  example)  so  that  the  point  at  which  the  aft  shock  wave  begins  and 
the  front  shock  wave  begins  coincide  vertically.  Next,  perform  an  al¬ 
gebraic-graphical  subtraction  in  which  values  below  the  "zero"  AP 
line  are  treated  as  negative.  Hils  procedure,  illustrated  in  Fig.  47, 
shoiild  restore  the  original  form  of  the  N-wave.  Such  is  Indeed  the 
case  for  this  signature  and  for  a  number  of  others  that  this  experiment 
has  been  tried  upon.  It  appears,  therefore,  that  turbulence  is  scat¬ 
tering  the  energy  of  the  sonic  boom. 

An  adequate  description  of  turbulence  can  only  be  made  in 
]  statistical  teims.  For  this  reason  treatment  of  data  frem  sonic  boom 

j  tests  In  ^Ich  the  atmosphere  was  turbulent  will  require  some  sort  of 

boatlstlcal  treatment. 
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Fig.  46  Comparison  of  Latoral  Distribution  of  Shock  Front  Overpressures. 
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Fig.  47  Rtconstrucfion  of  Wove  Signaiuro  by  Graphical  Mathod. 

Usually,  the  first  attempt  to  ti-eat  data  Involves  the  use 
of  the  normal  statistical  distrihution.  The  normal  curve,  which  was 
first  developed  ’oy  deMoivre  in  1753  has  heen  studied  extensively  by 
many  mathematicians.  Its  properties  such  as  the  standard  deviation, 
which  is  a  measure  of  the  dispersion  of  the  data,  the  skewness  which  is 
a  measure  of  the  asymmetry  of  the  data  and  the  kurtosls,  which  is  a 
measure  of  the  peakness  of  the  data,  are  well  understood.  It  is  impor¬ 
tant  to  either  show  that  the  data  is  normally  distributed  or  find  an 
appropriate  transfoimatlon  uo  convert  the  observations  to  a  normal  dis¬ 
tribution.  Otherwise,  the  mathematical  tools  such  as  those  mentioned 
above  provide  a  poor  description  of  the  data  and  a  loss  in  the  under¬ 
standing  of  relative  -Importance  of  the  meaningful  physical  parameters. 
This  in  tu-Ti  leads  to  poor  or  invalid  results  when  the  statistics  are 
used  for  value  Judgements. 
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The  normal  distribution  has  a  number  of  disadvantages 
when  it  is  applied  to  the  statistical  study  of  a  set  of  measured  sonic 
bccm  overpressures.  First,  it  is  an  infinite  distribution  and,  if  used 
for  estimating  the  upper  and  lower  limits  of  overpressure,  will  predict 
negative  "peak  positive  oveiiiressui'es"  as  well  as  Infinite  peak  over¬ 
pressures  if  the  confidence  limits  are  set  hl^  enou^.  A  negative 
"peak  positive  overpressure"  is  obviously  physically  unreal,  even  a 
zero  peak  positive  overpressure  would  require  either  locally  complete 
reflection  or  refraction.  The  infinite  peak  o-verpressure  is  similarly 
physically  unrealistic. 

Since  the  action  of  turbulence  is  to  locally  redirect 
(i.e.  scatter)  the  energy  from  the  N-wave  according  to  some  angular 
pattern,  a  consequence  of  this  effect,  when  the  centers  of  scattering 
action  are  randomly  distributed  in  space  above  a  two  dimensional  samp¬ 
ling  array,  is  that  data  will  be  logarithmically  distributed.  This  is 
shown  for  the  case  of  scattering  by  a  turbulent  temperatui'e  field  in 
Appendix  IX.  It  therefore  appears  that  the  data  should  be  handled  by  a 
transformation  that  has  the  following  properties: 

•  It  should  have  an  upper  and  lower  bound. 

•  It  should  be  of  a  logarithmic  form. 

The  four  parameter  lognormal  distribution  has  these  de¬ 
sired  characteristics.  A  complete  treatment  of  the  lognormal  distri¬ 
bution  may  be  fotand  in  Refs.  63,  64,  65,  and  66.  The  general  proper¬ 
ties  of  a  lognormal  dis^ ribution  may  be  considered  in  terms  of  the 
nranber  of  parameters  involved.  /.  parameter  here  means  the  number  of 
quantities  that  are  necessary  to  describe  the  distribution.  The  min¬ 
imum  number  of  parameters  is  twoj  they  are  the  mean,  and  the  variance, 
a 2  (the  square  of  the  standard  deviation).  It  is  to  be  emphasized  that 
the  variate  cannot  assume  zero  values,  since  the  transformation  I Y=loglX 
is  not  defined  for  X  =  O . 

If  there  are  physical  reasons  to  believe  that  a  lower  or 
upper  bound  exists,  then  the  three  parameter  log-normal  distribution  in 
which  the  transformation,  Y  =  log  (X±  Y  ),  is  the  appropriate  one.  This 
results  in  either  a  threshold  or  an  upper  bound  of  the  distribution  Y  . 
Ths  three  parameters  are  then  the  mean,  the  variance,  and  either  the 
upper  Euad  lower  bound  as  appropriate.  Further,  if  there  are  physical 
reasf ns  (such  as  with  the  sonic  boom)  to  believe  that  there  is  some  upper 
bound  and  a  lower  bound  that  the  values  of  the  parameter  may  assume,  then 
the  appropriate  substitution  is  Y  =  log  (X  -  Y  )  and  we  have  the  four 

{•  ■  X) 

parameter  log-normal  distribution  in  which  the  parameters  are  the  mesua, 
the  variance,  the  lower  bound,  #,and  the  upper  bound,  Y. 

The  higher  order  log-normal  distributions  are  difficult  to 
study  analytically,  particularly  when  the  upper  and/or  lower  bounds  are 
not  known.  Fortunately,  however,  a  special  type  of  graph  paper,  i.e. 
logarithmic  probability  paper  (K  &  E  468043)  can  be  used  with  great  ease. 
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The  procedure  for  using  It  is  as  follows:  First,  the  ohser- 
vationa  are  ai’ranged  in  order  of  increasing  value.  Second,  the  percent¬ 
age  of  observations  less  than,  or  equal  to,  a  given  value  is  computed. 

This  value  is  then  plotted  on  the  abscissa  against  the  value  of  the  vari¬ 
ate  on  the  logarithmic  ordinate.  The  mean  and  the  standard  deviation  now 
can  be  easily  determined  and  the  upper  and  lower  bounds  estimated. 

This  method  of  statistical  analysis  is  Illustrated  in  Fig. 

W  for  a  specific  set  of  data.  First,  the  data  was  ordered  and  the 
corresponding  percentages  were  found.  The  data  were  plotted  on  log¬ 
normal  probability  paper  according  to  the  procedures  described  above. 

It  is  apparent  from  these  data  that  a  simple  strai^t  line  of  the  two 
parameter  lognormal  distribution  will  not  adequately  represent  the  data. 

It  appears  that  there  may  be  either  three  discrete  sets  of  data  inter¬ 
mixed  or  that  there  are  seme  sampling  fluctuations.  However,  since 
there  are  physical  reasons  to  believe  that  upper  and  lower  bounds  exist, 
any  curve  throu^  the  data  must  be  asymptotic  to  these  values.  If 
statistical  homogeneity  is  assumed,  then  a  smooth  "S”  curve  similar  to 
the  Inverse  tangent  curve  (Ref.  64)  symmetric  about  the  50  percentile 
value  should  represent  the  data.  *15118  is  characteristic  of  the  four 
parameter  type  lognormal  distribution.  The  curve  Illustrated  in  the 
Figure  represents  a  fit  of  this  type  curve  with  an  upper  bound  approxi¬ 
mately  1.8  psf  and  a  lower  bound  near  0.3  psf.  The  median  value  is  0.7 
psf  -vdille  a  one  a  standard  deviation  about  the  median  is  (+0.5,  -O.3)  psf. 


Fig.  48  Statistical  Distribution  of  Overpressure  Data. 


The  varlahlllty  of  data  obtained  from  a  series  of  sonic 
boom  measurements  is  a  function  of  a  number  of  parameters.  Among  these 
are  VEiriations  due  to  changes  in  the  large  scale  atmospheric  flow,  varia¬ 
tion  due  to  changes  in  the  type  of  aircraft  used,  changes  in  Mach  numbers, 
changes  in  altitude,  and  differences  in  intensity  of  the  low  level  tur¬ 
bulent  flow. 

To  eliminate  as  many  variables  as  possible,  data  for  the 
F-lOl^A  aircraft,  at  an  altitude  of  28,000  feet,  flying  at  Mach  1,5>  was 
chosen  for  this  initial  study.  Further,  since  it  is  apparent  tint  clouds 
will  present  a  turbulent  field  of  different  character  than  turbulence 
near  the  ground,  only  those  booms  at  times  with  less  them  3/10  clouds 
were  considered.  The  wind  records  for  Oklahoma  City  were  examined,  and 
in  no  case  was  the  surface  wind  less  than  10  knots  for  these  observations. 
This  implies  that  the  atmosphere  was  always  somewhat  turbulent  near  the 
ground.  The  data  was  grouped  at  times  of  O7OO,  O9OO,  1100,  and  13OO 
Central  Standard  Time,  in  order  to  classify  the  data  as  nearly  as  possible, 
according  to  turbulent  intensity.  Hiis  is  plotted  in  Figs.  49,  50,  51> 
and  52. 

Seve'.'al  things  are  apparent  from  these  figures. 

•  The  variability  is  greatest  in  the  afternoon  as  would  be  expected. 

•  The  variability  is  a  function  of  the  horizontal  distance  of  the  ob¬ 
servation  from  the  flight  path.  Even  at  0700,  Test  House  4  at  10 
miles  from  the  flight  path  exhibits  large  variability. 

«  Some  of  the  data  seems  to  almost  indicate  that  there  were  three  dis¬ 
crete  sample  populations,  one  at  high  overpressures,  one  near  normal, 
and  at  low  overpressure.  This  is  particularly  evident  at  Test  House  4. 

•  There  appears  to  be  an  upper  bound  near  a  value  of  two  times  the  mean 
of  the  observed  overpressure  and  a  lower  bound  near  0.3,  the  mean  of 
the  observed  overpressure. 

It  is  probable  that  the  angle  at  which  the  sonic  boom  is 
propagating  is  important.  This  is  due  to  the  fact  that  the  turbulent 
power  is  greatest  in  the  horizontal  plane  (Ref.  4l,  p.  l6lff ). 

(C)  SUIvWARY  OF  RESULTS  -  The  theory  used  in  this  report  was  com¬ 
pared  with  some  of  the  test  data  obtained  during  the  Oklahcxna  City  sonic 
boom  tests.  It  was  found  that  when  observed  undeformed  pressure  signa¬ 
tures  were  compared  to  those  predicted  by  the  theory  the  agt-eement  was 
quite  good.  The  pi^essure  rise,  trace  length,  and  slope  of  the  expansion 
region  between  the  shocks  agreed  closely  with  the  predicted  wave.  How¬ 
ever,  the  pressure  rise  at  the  shock  waves  of  the  measured  -trace  was  not 
instantaneous  as  predicted  by  the  theory.  This  results  in  some  disa¬ 
greement  betwen  the  theory  and  test  when  only  the  front  shock  over¬ 
pressures  are  compared.  Some  of  the  deformed  pressure  waves  were  ana¬ 
lyzed  statistically.  It  was  found,  frcm  this  analysis,  that  the  impor¬ 
tant  scattering  parameters  are  the  angle  of  the  path  of  propagation  of 
the  shock  wave,  and  the  time  of  day  as  related  to  the  turbulent  inten¬ 
sity  near  the  grorind. 
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SECTION  I  FREQUENCY  OF  OCCURRENCE  OF  UNUSUAL  PROPAGATION 
CONDITIONS 


Any  discussion  of  the  frequency  of  occurrenbe  of  anomalies  in  the 
sonic  Doom  overpressure  and  distribution  must  include  an  assumption  about 
the  flight  profile  of  the  airplane.  In  the  following  discussion  a  typ¬ 
ical  supersonic  transport  is  assumed  which  first  exceeds  Mach  1.0  at 
altitudes  near  UO, 000  feet  above  the  ground.  No  supersonic  flight  is 
planned  for  altitudes  below  this.  A  survey  of  wind  speed  as  a  functibn 
of  altitude  for  the  mid-latitudes  has  indicated  that  99*9  percent  of  the 
time  the  maximum  wind  speeds  occur  at  altitudes  very  near  or  below  40, 000 
feet. 

The  cut-off  and  possible  focusing  phenomena  were  shown  in  Section  III 
to  occur  at  low  Mach  numbers  for  physically  possible  wind  speeds.  For 
temperature  profiles  with  an  inversion  near  the  ground,  cut-off  and  pos¬ 
sible  focusing  will  occur  above  the  ground.  Thus,  this  phenomena  may 
occur  when  no  inversion  exist  near  the  ground  (i.e.,  approximately  50 
percent  of  the  time).  From  this  discussion,  it  would  appear  that  cut¬ 
off  and  possible  focusing  under  or  to  the  side  of  the  flight  track  would 
occur  at  least  once  on  the  ground  during  the  low  Mach  number  portion  of 
the  flight  when  no  temperature  Inversions  exist  near  the  ground  (i.e., 
approximately  50  percent  of  the  total  available  time).  Whether  or  not 
this  phenomena  constitutes  a  significant  increase  in  the  boom  strength 
on  the  ground  is  still  not  resolved  because  of  the  question  of  the  vari¬ 
ation  of  the  reflection  factor  near  cut-off. 

Extreme  lateral  spread  is  possible  only  in  the  presence  of  very  strong 
winds  with  maximum  speeds  above  150  feet  per  second.  However,  as  these 
maximum  speeds  exist  below  40,000  feet  99*9  percent  of  the  time,  flight 
above  this  altitude  would  not  produce  this  phenomena.  It  would  appear 
that  because  of  the  high  altitudes  selected  for  supersonic  portions  of 
the  flight  the  sonic  boom  will  not  extend  to  extreme  lateial  distances 
to  the  side  of  the  flight  track. 

Considering  the  flight  profiles  contemplated  for  the  supersonic  trans¬ 
port,  it  becomes  apparent  that  only  the  boom  produced  during  the  first 
and  last  few  minutes  of  supersonic  flight  may  be  significantly  influenced 
by  the  meteorological  conditions  between  the  airplane  and  the  ground. 

The  boom  produced  during  the  remainder  of  the  flight  would  be  relatively 
unaffected  by  the  meteorological  conditions  between  the  airplane  and 
the  ground. 
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SECTION  SI  EVALUATION  OF  THE  ACCURAa  OF  THE  THEORY 


The  correction  factor  presently  used  to  account  for  the  Influence  of 
the  variable  properties  of  the  atmosphere,  yPaPg,  seems  to  be  somewhat 
less  than  predicted  by  the  present  theory.  ^Comparisons  with  Okladioma 
City  sonic  boom  data  Indicate  that  the  accuracy  of  the  theory  and  method 
of  Refs.  6~6  1b  quite  good  In  accounting  for  the  influence  of  the  met» 
eorologlcal  conditions.  However,  this  data  Is  very  limited  In  scope 
because  the  test  airplane  altitudes  were  relatively  where  the  differ¬ 
ences  between  the  more  sophisticated  approach  and  yPaPg  are  smedl. 

A  limited  amount  of  high  altitude  fll^t  test  data  is  available  and  has 
been  analyzed.  These  data  (Ref.  8o)  were  ccnpared  In  Fig.  53  with  the 
theoretical  predictions  vising  the  present  theory.  Here  again,  the  theory 
Is  In  close  agreement  with  the- test  data.  However,  more  comparisons  of 
this  sort  would  be  very  deslr^le. 
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It  was  shown  in  Section  II  that  rether  large  variations  in  the  pro¬ 
perties  of  a  stratified  stmosi^ere  ceuse  Sinall  changes  (±5  percent)  in 
the  sonic  hooo  strength  generated  by  airplanes  flylrg  at  Mach  niunhera 
above  I.3.  niese  sajne  variations  produced  larger  changes  (approximately 
±20  percent)  for  flight  at  Mach  numbers  below  1.3.  In  terms  of  the  me¬ 
teorological  observations,  these  atmospheric  variations  encompass  a 
range  in  the  mean  absolute  temperatxire  (which  is  used  to  de^^ermlne  the 
pressure  altitude)  over  the  altitudes  of  interest  of  6.6  percent.  Errors 
of  this  magnitude  in  any  single  observatl.on  are  very  rare.  Further, 
when  it  is  realized  that  £iny  single  observation  is  canpared  with  nearby 
observations  during  the  analysis  of  the  atmosph(‘ric  flow  field  the 
probability  of  such  a  large  error  going  undetected  is  almost  zero.  A 
similar  comparison  holds  for  the  effect  of  errors  in  the  measurement  of 
wind.  Thus,  it  would  be  expected  that  small  errors  in  the  dally  ob¬ 
servations  (say  1  percent)  would  cause  small  changes  (±5  percent  at  the 
maximum)  in  the  theoretical  predictions  even  for  the  low  Mach  number 
range  of  fll^t.  Based  on  the  results  in  Section  II,  then,  it  would 
appear  that  data  from  present  day  upper  air  sounding  are  sufficiently 
accurate  to  be  used  in  planning  flints  for  the  sup-ersonlc  treuisport. 
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SECTION  M  CONCLUSIONS 


The  effects  of  varying  meteorological  conditions  on  the  intensity 
and  spread  of  the  sonic  boom  have  been  investigated.  This  has  been 
accomplished  by  constructing  several  stratified  atmospheric  models 
and  comparing  sonic  boom  calculatious  in  these  atmosphere  vlth  results 
in  th.  U.  S.  Standard  Atmosphere,  1962. 


It  has  been  found  that  the  influence  of  wind  and  temperature 
variation  from  standard  conditions  is  primarily  a  function  of  Mach 
number,  wh'.le  the  influence  of  pressure  variation  is  independent  of 
Mach  number.  It  was  also  determined  that  for  flight  at  Mach  numbers 
above  1.’  the  largest  influence  of  varying  the  meteorological  conditions 
fror.  those  in  the  Standard  Atmosphere  is  a  change  in  the  sonic  boom  over¬ 
pressure  of  about  ±5  percent.  Flight  at  Mach  numbers  below  1.3  may  re¬ 
sult  in  more  si.piificapt  variations. 

The  meteorological  conditions  required  to  produce  focusing,  com¬ 
plete  cutoff,  extreme  lateral  spread,  and  deformation  of  the  pressure- 
wave  signature  have  been  investigated  and  methods  for  predicting  the' 
occurence  of  these  conditions  have  been  established.  It  has  been  found 
that  realistic  variations  in  temperature  and  winds  could  produce  focus¬ 
ing  or  complete  cutoff  for  flight  at  Mach  numbers  below  1.3.  Focusing 
would  occur  simultaneously  with  cutoff  where  the  shock  waves  are  normal 
to  the  ground ,  and  the  normal  doubling  of  the  overpressure  due  to  oblique 
shock  wave  reflections  would  not  occur  in  this  region.  Extreme  lateral 
spread  and  deformation  of  the  pressure  wave  due  to  interactions  with 
turbulence  may  occur  at  all  Mach  numbers.  The  former  would  not  occur, 
however,  for  flight  at  altitudes  above  those  where  the  maximum  winds 
exist,  regardless  of  the  Mach  number. 

The  flight  path  contemplated  for  a  typical  commercial  supersonic 
transport  includes  extremely  rapid  transition  through  Mach  numbers  near 
1.0  at  altitudes  near  or  above  40,000  feet  above  the  ground.  The  time 
spent  in  accelerating  through  Mach  I.3  during  the  initial  phases  of 
supersonic  flight  and  in  decelerating  through  Mach  I.3  during  the  latter 
phases  of  supersonic  flight  amounts  to  only  a  few  minutes.  In  this 
respect,  then,  it  becomes  apparent  that  for  the  Supersonic  transport 
only  the  boom  produced  during  the  first  and  last  few  minutes  of  super¬ 
sonic  flight  may  be  significantly  influenced  by  the  meteorological  con¬ 
ditions  between  the  airplane  and  the  ground.  Furthermore,  because  the 
supersonic  flight  altitudes  arc  generally  above  those  where  the  maxi¬ 
mum  wind  speeds  exist,  the  probability  of  the  occurrence  of  extreme  lat¬ 
eral  spread  would  be  very  small. 

A  number  of  comparisons  have  been  made  with  flight  test  data  ob¬ 
tained  during  the  Oklahoma  City  flight  test  series.  Predicted  and 
measured  pressure-wave  signatures  have  been  compared  for  each  test  air¬ 
plane.  It  was  generally  found  that  the  pressure  rise,  length,  and  slope 
of  the  expansion  region  between  the  shocks  of  the  measured  wave  agreed 
closely  with  the  predicted  wave.  The  pressure  rise  of  the  observed 
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shock  va-’ss  Is  not  instantaneous  as  predicted  "by  the  theory.  Tiiis  has 
led  to  some  disagreement  betveen  theory  and  test  data  vhen  only  front 
shock  overpressures  are  considered.  !nie  presence  of  turbulence  near 
the  ground  results  in  the  deformation  of  the  incoming  pressure-wave 
signature  euid  some  of  these  deformed  signatures  have  been  analyzed.  A 
statistical  analysis  of  these  data  has  Indicated  that  the  important 
scattering  parameters  are  the  angle  of  the  path  of  the  shock  wave  and 
the  time  of  day  as  related  to  the  turbulent  Intensity  n.^ar  the  ground. 

The  effect  of  viscosity  on  the  sonic  boom  has  not  been  Included,  in 
these  studies.  Viscosity  does  not  significantly  Influence  very  low 
frequency  waves  which  would  constitute  the  fundeusental  harmonics  of  the 
sonic  boom  N-wave  (see  Refs.  05  euid  86).  Thus,  this  would  be  of  minor 
importance  when  investigating  the  propagation  of  sonic  boom  through  the 
atmosphere . 


SECTION  3QII  RECOMMENDATIONS  FOR  CONTINUED  THEORETICAL  AND 
EXPERIMENTAL  WORK 

A  muoiber  of  aspects  of  the  sonic  boom  problem  still  remain  unresolv¬ 
ed.  Although  a  substantial  amount  of  effort  has  already  been  expended 
in  obtaining  an  understanding  of  the  problem,  further  theoretical  and 
experimental  work  aimed  at  Improvement  and  extension  of  the  present 
theory  would  be  desirable.  The  inirpose  of  this  section  is  to  discuss 
some  of  the  unresolved  areas  and  to  outline  possible  theoretical  and 
experimental  approaches  which  might  be  taken  in  seeking  solutions. 

1(A)  THEORETICAL -The  theory  developed  to  date  seems  to  be  suffi¬ 
cient  for  estimating  the  shock  wave  overpressures  and  pressure  signa¬ 
tures  received  on  the  ground  for  an  airplane  in  steady  level  flight 
through  a  nonuniform  nearly  stratified  atmosphere.  Comparisons  with 
available  experiment  such  as  those  presented  in  Section  IV  have  led  to 
this  conclusion.  However,  the  theoretical  development  should  be  extended 
to  Include  the  effects  of  general  aircraft  maneuver  and  the  turbulence 
effect  near  the  ground.  Some  preliminary  work  has  been  done  in  both 
these  areas.  The  latter  was  outlined  In  Section  III.C.  Further  work 
should  be  directed  toward  seeking  a  solution  of  the  pressure  wave  his¬ 
tory  near  and  beyond  cut-off,  and  toward  an  understanding  of  the  effect 
of  Interactions  of  the  pressure  wave  with  regions  of  high  altitude  tur¬ 
bulence  such  as  exist  In  towering  cumulus  clouds.  Guidelines  for  these 
theoretical  investlga-ions  are  outlined  in  the  following  material. 

(1)  Extvnsien  to  Gtnoral  Monovvors  -  The  prediction  of  sonic  boOD 
overpressures  on  the  ground  for  eui  airplane  engaged  In  general  maneuvers 
in  a  nonuniform  atmosphere  Is  a  quite  coupllcated  task.  The  ray  tube 
area  concept,  used  In  Hefs.  6  throu£^  8,  lends  Itself  to  considerations 
of  this  sort.  The  expression  for  the  ray  tube  area  (for  instance  that 
presented  in  Appendix  II )  may  be  modified  by  the  inclusion  of  terms 
which  describe  the  distortion  of  the  ray  paths  due  to  the  maneuver. 

Some  preliminary  work  has  been  done  In  Ref.  77  toward  the  development 
of  the  required  terms. 

In  addition  to  developing  these  terms,  one  of  the  major 
problems  Involved  in  making  sonic  boom  estimates  for  a  maneuvering 
airplane  Is  the  obtaining  of  a  description  of  the  shock  front  intersection 
with  the  ground,  euid  the  location  of  the  points  at  which  this  front  may 
cusp  (or  fold  over  on  Itself)  to  form  a  focused  so  called  "super"  boom. 

The  transformations  between  the  ray  Intersections  with  the  ground  and 
the  shock  front  location  are  extremely  complicated  because  of  flight  path 
and  velocity  variation,  and  because  of  variation  in  the  atmospheric  prop¬ 
erties.  At  the  present  time,  it  is  felt  that  the  shock  front  can  be  lo¬ 
cated  by  trial  and  error  (Ref.  6).  Furthermore,  the  ray  tube  area  con¬ 
cept  cannot  account  for  the  situation  when  two  rays  generated  by  the 
airpleuie  at  largely  different  times  reach  the  ground  at  the  same  point 
emd  time  (shock  front  crossover).  This  situation  could  occur  when  the 
airplane  is  engaged  in  a  short  radius  circular  turn  (Fig.  7,  Ref.  Jd). 
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Both  the  cusp  and  the  shock  front  crossover  are  regions  which  sbovild  re¬ 
ceive  special  attention  when  investigating  the  shock  strength  distribu¬ 
tion  on  the  ground.  The  work  which  has  been  referenced  above  would 
prove  useful  in  any  future  theoretical  development, 

(2)  Sriock  History  At  and  Boyond  Cut-oH  -  The  region  near  the  cut-off 
of  the  boom  is  not  fully  understood.  In  aerodyiie.mic  terms  it  consists  of 
a  region  of  mixed  flow  where  in  an  area  slightly  before  cut-off  the  flow 
is  supersonic  and  in  an  area  beyond  cut-off  the  flow  is  subsonic.  In 
physical  terms  this  means  that  a  shock  wave  cannot  exist  beyond  cut-off. 
However,  some  sort  of  pressure  distribution  will  be  propagated  through 
the  air.  The  type  of  disturbance  is  properly  illustrated  in  Fig.  5* 

A  number  of  approaches  may  be  used  in  seeking  a  solution  to 
the  shock  wave  history  in  this  region.  One  of  these  might  be  to  extend 
the  Whitham  theory  of  Ref,.  2  to  account  for  variable  speed  of  sound  in 
the  supersonic  region.  This  would  allow  the  description  of  the  pres¬ 
sure  perturbations  generated  by  the  airplane  to  be  described  up  to  the 
cut-off.  The  subsonic  description  might  be  obtained  by  the  classical 
methods  outlined  by  Prandtl  (Ref.  79) •  The  two  solutions  could  then  be 
investigated  in  the  region  of  the  cut-off  to  describe  the  transition  of 
the  pressure  distribution  from  the  supersonic  shock  wave  to  the  sub¬ 
sonic  pressure  wave.  The  variation  of  the  reflection  factor, K^,  should 
also  be  investigated  in  the  same  region. 

(3)  )  Extantien  of  Analysis  of  Turbalant  Efftcts  -  T^e  theory  of  scattering 
of  acoustic  and  shock  waves  inherently  postulates  a  weak  interaction  be¬ 
tween  the  acoustic  or  shock  waves  and  the  turbulent  field.  This  allows 
the  assumption,  that  the  wave  length  of  the  acoustic  energy  is  much  less 
than  the  wave  lengths  of  the  turb\ilence  to  be  made  and  results  in  con¬ 
siderable  simplification  in  the  theory.  At  the  wave  lengths  associated 
with  the  sonic  booms  produced  to  date,  this  approximation  is  of  extremely 
dubious  validity  and  will  become  even  more  so  for  larger  aircraft. 

Two  possible  alternative  explanations  may  be  offered  to  ex¬ 
plain  the  spikes  observed  in  the  N-wave  traces  such  as  those  from  Okla¬ 
homa  City.  First,  they  may  be  the  result  of  coherent  scattering  in  which 
the  phases  of  the  waves  combine  at  certain  points  so  as  to  give  rein¬ 
forcement.  Coherent  scattering  can  occur  when  the  interaction  is  weak. 
Alternatively,  the  spikes  may  be  due  to  strong,  resonant  interaction  be¬ 
tween  the  shock  waves  and  turbulent  eddies.  This  type  of  interaction 
is  extremely  nonlinear  and  the  terms  usually  neglected  in  the  partial 
differential  equations  must  be  retained.  In  addition,  if  there  is  strong 
interaction  between  the  N-waves  auid  the  turbulence,  the  separation  of 
the  cross-sections  into  one  due  to  temperature  and  one  due  to  wind  is 
incorrect.  Consequently,  effects  of  a  random  temperature  field  and  a 
random  wind  field  must  be  considered  simultaneously.  The  resolution  of 
the  problem  of  the  mode  of  interaction  between  the  sonic  bocm  and”  tur¬ 
bulence  can  only  be  made  by  appeal  to  experiment,  since  the  only  "a  priori" 
reason  for  preferring  one  scattering  mode  to  the  other,  is  the  require¬ 
ment  that  the  eddies  have  an  ordered,  spatial  distribution  for  coherent 
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scattering.  However,  em  exact  physical  description  of  turbulence  is 
lacking  at  the  present  time  eind  this  requirement  for  Ein  ordered  distribu¬ 
tion  cannot  be  either  used  or  Rejected. 

Appeal  to  such  experiments  as  the  closely  spaced  microphone 
arrays  used  during  some  of  the  Oklahcna  City  sonic  booa  tests  tentatively 
indicate  that  coherent  scattering  does  not  occur.  For  instemce,  if  the 
records  frcm  Mobile  1,  Flight  2  on  2  April  1964  are  considered,  spikes 
should  occur  at  interveils  at  least  every  196  feet  (the  wave  length  of 
the  fundamentel.  harmonic  of  the  "free  air"  N-wave),  or  very  close  to 
every  microphone  for  the  200  foot  spacing  of  this  array.  This  rises  from 
the  requirement  considered  in  Section  ni.C  for  a  resonant  peak  at  in¬ 
tervals  of  (m^)  where  m  is  an  integer  and  X  is  the  wave  length  of  the 
fundamental  N-wave  harmonic.  This  is  not  the  case.  Rather,  a  regular 
progression  from  a  peak  value  of  overpressure  to  a  minimum  is  observed 
over  a  distance  of  80O  feet.  Coherent  scattering  can  thus  be  tentative¬ 
ly  ruled  out  in  this  case.  The  f\irther  observation  can  be  made  that  the 
cross-section  (i.e.  the  area  affected  by  the  eddy)  must  be  at  least  570 
feet  in  radius;  possibly  as  much  as  ll4o  feet.  For  these  sizes  it  is 
implied  that  there  is  seme  degree  of  resonance. 

For  this  reason  it  is  suggested  that  the  analysis  of  the 
interaction  of  the  sonic  boom  be  extended  as  follows: 

•  Extend  the  mathematical  analysis  to  include  the  study  of  strong 
scattering.  Green's  fxinctlons,  the  Fourier  transform  (or  possibly 
the  S'dlBtensional  LaFlace  transform)  and  of  spherical  ^rmonles  can 
be  profitably  applied  to  this  study^  since  a  more  general  type  of 
solution  can  be  generated  by  this  technique.  Currently,  initial  steps 
have  been  made  in  this  atialysls,  but  are  not  sufficiently  complete  to 
include  in  this  report. 

•  Analysis  of  a  representative  sample  of  the  Oklahoma  City  mlcrometeoro- 
loglcal  data  to  determine  the  spectral  densities  of  temperature  and  wind. 

•  Consideration  of  the  nonhomogenous  turbulent- scattering  problem  in 
which  reverberations  can  occur. 

Providing  an  analytical  solution  of  this  problem  can  be 
achieved,  reasonably  exact  statements  about  the  spatial  distribution  and 
the  maximum  values  of  the  overpressure  as  well  as  the  wave-length  of  the 
spikes,  can  be  expected  for  the  case  of  single  scattering. 

Multiple  scattering  is  a  very  difficult  problem  and  progress 
will  be  slow  in  this  case. 

The  problem  of  propagation  throiigh  convective  clouds  should 
be  studied  to  some  extent  since  einanalous  propagation  almost  surely  oc¬ 
curs..  The  internal  structure  of  convective  clouds  is  sufficiently  well 
known  that  a  preliminary  study  cem  be  made,  neglecting  turbulence. 

(B)  EXPERIMENTAL  -  AddltloneJ.  experimental  confirmation  of  the 
theories  presented  In  this  report  would  be  desirable,  especially  near  cut¬ 
off.  Controlled  experimental  data  would  also  be  helpful  In  understanding 
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more  about  the  effect  of  turbulence  on  the  pressure  signature.  A  test 
progrea  Is  briefly  outlined  In  the  foUovrtng  material  which  might  yield 
useful  data  in  both  of  these  ai'eas. 
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(1)  MastsKiMirt  sf  Sksck  Str«n|tii  At  and  Btyond  Cut>oH,  Sttcdy  Flight  > 

It  was  noted  in  Section  III  that  local  Intensification  of  the  boom  may 
occur  or.jy  simultaneously  with  cut-off.  It  was  further  observed  that 
the  reflection  factor,  Kr,  might  also  vary  in  the  same  region.  A 
series  of  field  tests  could  possibly  verify  both  of  these  postulates. 

The  theory  and  method  developed  in  Refs.  6  through  8  in¬ 
dicate  that  the  free  air  overpressure  (i.e.  accounting  for  no  reflection) 
may  increase  rapidly  in  a  very  small  region  near  cut-off.  This  is  il¬ 
lustrated  in  Fig.  $4  for  the  free  air  overpressure  under  the  flight  track 
of  an  airplane  in  steady  flight  in  a  headwind. 

The  figure  shows  that  the  theoretical  intensification  of  the 
free  air  overpressure  takes  place  over  a  very  short  distance.  If  an  air¬ 
plane  were  flown  over  an  instmaentcd  tower  at  nearly  the  cut-off  Mach 
nuaber  for  the  conditions  of  the  fli^t  the  variation  of  shock  with  dis¬ 
tance  might  be  measured.  These  meeisurements  could  be  ccmpared  with  theo¬ 
retical  predictions  for  the  same  meteorological  and  flight  conditions. 
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STANDARD  ATMOSPHERE 


M  -  1.30 


ALTITUDE -47,000  FT 
HEADWIND  AT  AIRPLANE  - 145.05  FPS 


HEADWIND  AT  GROUND 
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Uie  tower  height  should  be  in  excess  of  1000  feet  and  should 
contain  a  dense  array  of  microphones  at  regular  Intervals  aloug  Its 
height.  It  should  also  be  equipped  to  continuously  measure  wind  velo¬ 
city  and  temperature  at  regular  intervals  along  its  height.  The  surround¬ 
ing  terrain  should  be  relatively  flat  and  the  area  near  the  ground  should 
be  such  that  relatively  quiescent  conditions  exist  at  some  time  during 
the  day  or  night.  This  would  be  to  avoid  possible  interactions  with 
turbulence  near  the  ground  which  may  make  data  interpretation  difficult. 
Facilities  should  be  available  to  obtain  good  forecasts  and  measurements 
of  the  wind  and  temperatures  between  the  airplane  and  the  ground  so 
that  each  flight  may  be  programmed  to  obtain  cut-off  under  the  flight 
track  near  the  ground.  Facilities  should  also  be  available  for  tracking 
the  airplane  flight  path  and  coordinating  it  with  the  tower  instrumenta¬ 
tion.  The  flight  path  should  be  directly  over  the  tower  and  directed 
generally  into  the  wind.  One  possible  location  for  these  tests  might 
ba  the  BREN  tower  located  at  the  AEG  Nevaida  test  site. 

>  In  addition  to  the  free  air  microphones  a  set  of  microphones 
ml^t  be  mounted  on  reflecting  boemrds  at  the  same  reguleur  intervals  to 
obtain  the  variation  of  Kr  near  the  cut-off.  This  could  be  accomplished 
by  comparing  the  two  simiiltaneous  sets  of  measurements. 

The  same  type  of  facility  may  also  be  useful  in  obtaining 
measurements  of  the  turbulent  scattering  phenomena  at  veurious  levels  a- 
bove  the  ground.  Flights  could  be  made  over  the  tower  during  the  time 
of  day  when  there  is  the  greatest  probability  of  turbulent  activity  in 
the  area  near  the  ground.  Details  of  this  portion  of  the  test  program 
are  outlined  in  Section  V.B.2. 

(2)  Mcafurtmtntf  af  Prttfurt  Wav*  Dittortion  by  Low  Altitud*  Turbultnc*  - 
Present  indications  are  that  the  single  scattering  interaction  of  the 
sonic  boom  with  turbulence  is  most  important  from  an  overpressure  stand¬ 
point.  It  is  probably  limited  to  the  lower  ICXX)  to  3CXX)  feet  of  the 
atmosphere.  If  the  turbulence  extends  to  greater  depths,  the  probabil¬ 
ity  of  multiple  scattering  increases.  A  further  indication  is  that  the 
distribution  of  the  spike  amplitudes  in  single  scattering  is  trigonom¬ 
etrically  distributed.  It  would  be  very  desirable  to  determine  the 
effect  of  turbulence  at  various  levels  above  the  ground.  The  preferred 
setup  for  this  purpose  would  be  a  fixed  installation  such  as  a  hi^ 
tower  with  microphones  both  on  the  tower  and  on  the  ground  around  it, 
no  further  than  the  "free"  air  wave  length  of  the  sonic  boon  apeirt  and 
arranged  to  obtain  the  overpress\ares  over  an  area  on  the  order  of  eight 
wave  lengths  on  a  side.  The  tests  should  include  tests  at  various  angles 
of  incidence  of  the  K-wave,  and  imder  varying  wind  speeds  anfl  vf-H'ical 
temperature  lapse  retes.  It  must  be  pointed  out  that  a  series  of  nans 
will  be  necessary,  since  the  eddies  are  randomly  distributed. 

The  BREN  tower  at  the  AEG  Nevada  Test  Site  is  1527  feet 
higli,  is  instrumented  meteorologically  and  is  remotely  located  from 
inhabited  areas.  It  would  seem  to  be  an  excellent  facility  for  experi¬ 
ments  of  this  kind. 


ATMOSPHERIC  CORRECTION  FACTOR 


APPENDIX  I  VARIATION  OF  WITH  HEIGHT  OF  GROUND  ABOVE  MEAN  SEA 
LEVEL  (MSL) 


y 


The  location  of  the  ground  above  mean  sea  level  will  have  an  tffect  > 
on  the  atmospheric  correction  factor,  (see  Eq.  (2)).  The  curves  In  ii 
Pigs.  I-l,  1-2,  and  1-3  were  prepared  to  facilitate  prediction  of  this  ' 
effect  for  areas  of  the  country  which  lie  significantly  above  0  feet  i 
mean  sea  level.  | 

These  curves  were  developed  assuming  U.  S.  Standard  Atmosphere,  1962  Ij 
conditions  from  0-  feet  mean  sea  level  In  each  case.  Thus,  the  teoipera-  J 
ture  and  pressure  at  the  ground  In  each  figure  were  taken  as  the  standard  j 
values  at  2000,  UOOO,  and  6OOO  feet  above  mean  sea  level  respectively.  1 


AIRPLANE  HEIGHT  ABOVE  GROUND  -  1000  FT. 


Fig.  /-I  Affflo*p'.'«/7c  Cometion  Factor  for  tho  GrounJ  at  2000  foot  obovo  MSL. 
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ATVOSPHEfBC  CORRECTION  FACTOR 


ivMvwvxAW  ■  urn  wiM-ww  g^WtftmWVH/WWWI 


Fig.  1-2  Atmoaphtrlc  Corncflon  Focfer  for  t/m  Crovnd  at  4000  hot  abov*  MSL. 
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ATMOSPHERIC  CORRECTION  FACTOR 


16.0 


AIRPLANEHEIGHT  ABOVE  GROUND  -  1000  pT 


Fig.  1—3  Atmo$ph»ric  Correction  Factor  for  th*  GrovnJ  at  6000  foot  obov*  MSL. 
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APPENDIX  n  A  REVIEW  OF  THE  THEORY  FOR  SHOCK  WAVE  PROPAGATION 
THROUGH  A  NONUNIFORM  MEDIA  - 

TLe  purpose  of  this  section  Is  to  siaonarlze  the  results  of  the  work 
of  Refs.  6  through  8.  Specific  details  concerning  the  development  of 
this  theory  may  he  obtained  from  Ref.  6.  Basically,  the  results  describe 
second  order  perturbations  In  particle  velocity,  pressure,  emd  density 
behind  a  weak  shock  wave  In  terms  of  the  first  order  undisturbed  quan¬ 
tities.  Ihe  assungrtlons  are:  (l)  that  the  third  order  perturbations 
are  small  cctnpared  to  the  second  order  perturbations  so  that  they  may 
be  Ignored,  and  (2)  that  the  shock  wave  propagates  with  the  velocity  T 
vhlch  Is  a  function  of  Its  strength,  l.e. 


Furthermore,  It  Is  assumed  that  the  rays  which  describe  the  path  of 
propagation  of  the  shock  wave  and  the  ray  tubes,  which  are  formed  by 
bxmdles  of  the  rays,  propagate  Independently. 


Before  summarizing  the  results,  seme  nomenclature  will  be  established, 
inie  coordinate  system  shown  In  Fig.  II-l  will  be  used  In  this  and  the 
following  sections,  HI  through  VII. 


Fig.  JI-I  CoorJinatt  Systam  and  Wind  Axat. 


The  subscript  "R"  notes  quantities  In  the  (  6  ,  )  coorlnate  aystem 

rotated  throu^  the  angle  $  to  coincide  with  the  Initial  ray  direction. 
Starred  quantities  note  winds  relative  to  those  at  the  airplane  so  that 
U*  =U(  z  )  -  U(  z  =0).  The  subscripts  "a"  and  "g"  refer  to  the  value 
of  the  quantity  at  the  airplane  and  the  ground  respectively.  No  sub¬ 
script  Indicates  the  value  of  the  quantity  at  any  point  between  the  air¬ 
plane  and  the  ground. 


W3i  vw  vwvTi  «wTc  iTRsrii  tfwiniuii  viQfTVHviuncy^ir^i(yn(!yifliDGU>€J 


Ibe  following  is  a  review  of  the  equations  used  to  compute  the 
shock  wave  strength  In  nonhomogeneous  media.  Distances  are  measured 
from  the  airplane  and  the  coordinate  system  is  taken  to  move  with  the 
airplane  (and  hence  with  the  winds  at  the  airplane  altitude). 

Hie  equations  of  the  rays  in  the  (  ,  ii  ,  z  system  are  given  by: 


a-i-Us* 

dz  liR  ¥  Eq.  (Il-la) 


dri  _  Vr» 

dz 


Bq.  (Il-lb) 


IL-  1 

dz  iirI  (II-lc) 


Eq.  (II- Id) 


^R  = 


V.costf  —  Ur'* 


Eq.  (II- le) 


Eq.  (Il-lf) 
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Transformations  between  the  x,  y,  z  and  the  5  ,  t) 
ore  given  by: 


X  =  5  cos®  —  II  sinfl 


y  —  5  sin  fl  4-  I)  cos  e 


Ur*  =  U*  cos  fl  +  V*  sinfl 


Vp’'  =  -U*  sin  «  +  V*  cos  ® 

I 

I 


1, +(M=-,)cosV]-^' 

8inll  =  , 

[l+(M®— l)COe*q)]^^® 


Transformations  oe tween  a  fixed  ground  system  X, 
y>  z  system  are: 


X  =  X  +  U  (O)  t 
Y  =  y  +  V  (O)  t 
Z  =  z 


,  z  systems 


Eq.  (Il-2a) 


Eq,  (Il-2b) 


Eq.  (Il-2c) 


Eq.  (Il-2d) 


Eq.  (Il-2e) 


Eq.  (Il-2f) 


Z  and  the 


Eq.  (II- 3a) 
Eq.  (II- 3b) 
Eq.  (Il-3c) 
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where 


U  (O)  =  U(z=0),  etc. 


Some  necessary  auxiliary  relations  are  as  follows: 


U*  =  U  (2>  -  U  (O) 


Eq.  (Il-4a) 


V*  =  V  (2)-  V  (O) 


Eq.  (Il-4b) 


M  = 


a  (O) 


Eq.  (n-4c) 


a  (2)  =  a  (2)  1  + 


Jl  p_ 

t  P’(z) 


Eq.  (Il-4d) 


W*  (2)  = 


dri 

^R*  dz  Vr  dz 
ds 
dz 


Eq.  (Il-4e) 


where  a(z)  Is  the  ambient  sound  speed  In  the  atmosphere,  and 
P(z)  Is  the  ambient  pressure  In  the  atmosphere. 


The  shock  wave  may  be  located  on  the  ground  In  the  X,  Y,  Z, 
system  through  the  following  transfoiinations  for  steady  level  flight: 


=  *  +  0(O)t,- 


Eq.  (n-5a) 


If  ...Ok  =  y+v(o)t. 


Eq.  (II-51>) 


where  t  ■  t  for  the  q>  «»  -90®  ray  to  reach  the  ground. 

o  • 

nie  ray  tube  area  is  given  by: 


I 

I 

i  A(z)=S 


cos  e 


)nn 


+ 


V  * 


sin  e 


1/2 

Eq.  (II-6) 


where 


V,=  Ma(0) 


S 


The  shock  strength 


is  given  by: 


A  P 
P(z') 


(M  - 

Mz')B(z') 


(Y  +  1)1/2 


D(0)  KYo.tf) 


1/2 


'  •  dz 

dz 

(w*-l-a)*  B(z) 


1/2 


Eq.  (II-7) 


where 


B(z)  = 


A(zj  P(z) 
a(z) 


1/2 

1  Q(z) 


D(0)  = 


a(0) 


1/1 


P(0)=  P  (z  =  0),  etc. 


liumerlcal  solutions  may  be  obtained  by  an  iteration  procedure 
where  in  the  first  iteration  the  assumption  a(z)  -  a(z)  may  be  used. 
The  second  iteration  may  use  a(z)  computed  using  the  first  value  of 


A  P 


and  so  on. 


It  has  been  found  that  in  regions  away  from  cut-off 


and  focusing  the  iteration  converges  quite  rapidly. 


The  above  equations  are  developed  in  Ref.  6  (see  also  Refs.  7 
and  8).  The  equation  for  the  ray  tube  sirea,  given  here,  (Eq.  (II-6), 
is  an  improved,  more  accurate  version  of  the  aquation  in  Ref.  6  and  is 
developed  in  the  Appendix  of  Ref.  7» 


I 


Ihe  result  Is: 


AP  ■  Sonic  boom  overpressure 


Pg- 


Amblent  pressure  at  ground 
Ambient  pressure  at  airplane 
Airplane  macb  number 


Y  ■  Ratio  of  specific  heats 
I  (Yo.— SQT)  "■  Airplane  shock  strength  parameter 
gg  ■■  Sound  speed  at  ground 
Sound  speed  at  airplane 
A  ■  Ray  tube  eurea 

z  a  Distance  measured  down  from  airplane 
Sa  Distance  measured  along  ray  path 


!nhe  ray  path  and  ray  tube  area  are  primarily  functions  of  the  tem¬ 
perature  variation  In  the  atmosphere.  The  effect  of  variations  in 
ambient  pressure  can  be  obtained  by  allowing  only  the  pressure -height 
curve  to  vary,  and  by  dividing  the  above  equation  using  the  standard 


pr«Mure-hel^t  curve  by'  that  uaing  a  oonatandard  curve.  Because  the 
tel^rature-belc^t  curve  is  assuned  to  remain  constant,  the  result  is: 


P 

IlLl 

pg*/» 

<p.)Ji?. 


1 

f'hm  ^  d  Z 

/  •  dz  1 

I/I 

1 

J.  ..(4E)- 

■td. 

-  ds  j,. 

|l/2 

1 

As  variations  from  the  standard  pressure -helc^t  curve  In  the 
real- atmosphere  are  generally  stsall  It  may  be  assumed  that  the  Integrals 
vlU  remain  approximately  equal.  Ihus  the  relationship  becomes! 


Eq.  (4) 


Ihls  relationship  was  compared  with  computer  program  output  in 
Section  II.B.2  to  check  the  validity  of  the  assumptions.  It  vas  seen 
that  the  assuoptlons  are  acceptable. 
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APPENDIX  JZ  DEVELOPMENT  OF  GENERAL  EQUATION  FOR  LATERAL  CUT-OFF 
LOCATION  - 

Consider  a  general  temperature  profile  such  as  that  shown  In  Fig  IV-1, 
In  which  the  teiQ>erature  between  any  two  significant  levels  Is  assumed 
to  vary  linearly. 


F/f.  JX-1  Csnsra/  Tmnpcrotvrs  Proflh, 


For  6Lny  portion  of  this  profile  a  ray  path  Is  described  by  Eqs. 
(U-la)  and  (ll-lb).  Assuming  no  wind  the  first  of  these  becones: 


dz  Hr  \n|  / 


Further  the  lateral  (y)  displacement  of  the  ray  for  no  wind 
(Eq.  (Il-2b))  Is; 


y  =  6  sintf 


=  sintf 


i/j 


dz 


Qi^loylng  the  definition  of  Iln  frcn  Eg.  (Il»lf ),  ass\ming  that 
5=  a,  and  putting  V,  In  teims  of  Madi  number,  M»V,/a,pi  the  above 
become B : 


Integration  of  this  equation  can  be  facilitated  If  tb'  variable 
of  Integration,  z.  Is  put  In  terms  of  a/aap|.  Assuming  that  the  sound 
speed,  a,  also  varies  linearly  between  two  points  (valid  as  long  as  the 
lapse  rates  are  of  the  magnitude  normally  encountered)  the  relation 
between  a  and  z  may  be  written  as: 


where 


■  Sound  speed  at  level  z„ 

■  Sound  speed  lapse  rate  between  levels  z„  and  z^^| 

a, 

Zn+1  2b 


Using. the  above  expression 


dr= ^  d  (^) 

an  \a«Di  / 


so  that 


_  anpl  - 

y  =  07  sin  «y 


<str)  -a  . 


The  lateral  location  of  any  ray  is  given  by  the  summation  of  the 
above  integrals  for  each  significant  level  between  the  airplane  and  the 
ground.  Ibe  integration  and  sutmation  yields: 


y  =  a  I  Sin  « 


[(Mcoa  [(Mcoa 

2^ 


where  "j"  represents  the  values  at  the  last  significant  level 
(in  most  cases  the  ground).  IQie  last  ray  to  reach  the  ground  (or  last 
significant  level)  is  the  one  for  which  the  direction  cosine,  >  is 
equal  to  1.0  at  the  level  where  a>amax.  Using  Bq.  (ll-le)  the  angle6 

of  the  last  ray  is  given  by  (assuming  a  »  a,  smd  no  wind): 


^  =10  = 


a«pi 

M  cos 


Mcosfl= 


M  a^i 


:  \nrvi  YDtd  JClUr 


Using  these  expressions  and  the  definition  for  a„,  the  lateral 
location  of  the  last  ray  to  reach  the  last  significant  level  is  given 


Eq. (6) 


In  sotae  cases  the  temperature  between  two  levels  does  not  change^ 
as  in  the  stratosphere  of  the  standard  atmosphere.  Assuming  that  this 
occurs  between  levels  Zmandz^+i  (i.e.  am«am+i)  the  term  involving 
(zm+i— Zm)  would  become  indeterminate.  The  simplest  method  of  obtaining 
the  form  of  this  term  is  by  taking  the  limit  using  L' Hospital's  rule  as 
am— *am4i .  Thus,  this  term  becomes: 


To  Illustrate  the  application  of  this  relationship,  Eqs.  (6)  and 
(7)  liave  been  applied  to  the  U.S.  Standard  Atmosphere,  1962  (Ref.  1) 
whereamtx  would  occur  at  the  ground,  i.e.,ankx  ■  Ag.  Thus,  for  an  air¬ 
plane  flying  between  36,0\>0  fc'jt  and  the  groui.d.’; 


For  an  airplane  flying  between  36,000  and  65,000  feet: 


And  for  an  airplane  flying  above  65,000  feet. 
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APPENDIX  Z  DEVELOPMENTOF  FOCUSING  CRITERIA  UNDER  FLIGHT  TRACK 


The  ray  tube  area  expression,  Eq. Is  repeated  below: 


Eq.  (II-6) 


Ibe  shock  strength  will  Increase  if  this  quantity  vanishes  as 
It  appears  In  the  denominator  of  the  expression  for  the  overpressure 
E^*  (II-7)*  bube  area  can  tend  toward  zero  only  If  the  expres* 

slon  Inside  the  brackets  tends  toward  zero  as  S  Is  the  length  of  the 
ray.  SeLtlug  A  ■>  0,  speuializiag 'this 'for  uo'*wind,  and  investigating 
the  ray  under  the  airplane  (  0  ■  o)  this  becomes: 

nR=  0 


Hhe  quantity  Eq.  (ll-lf)  Is  the  direction  cosine  of  the  normal 
to  the  ohock  front  with  the  z  axis  and  Is  given  In  general  by; 


nR  = 


_ a_ 

COS0 


■=1^)1 


assuming  that  the  shock  front  propagates  at  nearly  the  local  speed  of 
sound.  When  this  quantity  Is  zero  the  shock  front  Is  normal  to  the 
ground.  Thus  for  no  wind  and  0  ■  0®,  A  can  tend  toward  zero  If  V,  ■  a . 
Putting  this  in  terms  of  the  Mach  number  and  specializing  it  further 
by  assuming  that  the  maximum  speed  of  sound  occurs  on  the  ground  (l.e., 
■  0  only  at  the  ground) 


Eq.  (9) 


Mfocuh  — 


From  this  discussion  it  can  be  seen  that  focusing  occurs  simul¬ 
taneously  with  cut-off  O). 


For  the  case  with  wind  and  #  >  0  It  is  evident  that: 


It  is  apparent  that  the  naxinum  of  (a  +  U)  must  occur  at 

the  ground  for  «  0  at  the  ground.  Hero  too,  crj^  -off  and  focusing 
occur  simultaneously. 


APPENDIX  n  DEVELOPMENT  OF  LATERAL  FOCUSING  CRITERIA 


Considering  the  entire  raj’  tube  area  expression  Eq.  (II-6)  it  is 
apparent  that  It  cam  go  to  zero  If,  and  only  If,  both  portions  of  the 
term  In  the  brackets  go  to  zero  eimultaneo\isly.  This  requires  that: 


V  * 

^  —  sin  tf  =  0 

V  a 


at  the  same  location.  Here  again  two  possibilities  are  open.  These  are 


U.?  V  *  V  * 


tnat  — ^  -  cos  d  B  0  when  _J1  -  sin  6  •  0,  or  n_  ■  0  when  _E.  -  sin  $  m  o, 

V,  V,  "  V, 


It  can  be  shown  that  the  first  of  these  leads  to 


^U**+V**  =  V, 


which  Is  physically  unlikely.  Investigating  the  second  a  simultaneous 
solution  must  be  sought  for: 


M  cos  9 - ^ 

.  a. 


—  sin  9  =  0 


A  solution  for  the  value  of  V*  required  In  terms  of  the  Mach  number, 
U*,  and  sound  speeds  can  be  obtained  1*"  the  following  transfonaatlons 
are  used: 


coao  = 


ainO  — 


r*cosfl  +  V*sinfl 

Eq.  (Il-2c) 

-0*^8111  4-V*co8fl 

Bq.  (Il-2d) 

1 

[l  +  (M*-l)C08*q)] 

|I/2 

e 

CVl 

1 

M 

M 

s 

(M'  — 1)‘^*C089 

[i  +  (M*-1)c08^9 

1/2 

Bq.  (n-2f) 

After  substituting  these  In  the  above  expressions,  equating  them, 
and  after  sane  algebraic  memlpulatlon  It  can  be  shown  that; 


a. 


1/1 

Eq.  (VI-1) 


This  relationship  was  used  In  developing  the  curves  In  Fig.  25*  Here 
again  it  Is  evident  fron  the  derivation  that  focusing  and  lateral  cut-off 
(in  this  case)  occur  slnriltaneously. 
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APPENDIX  m  DEVELOPMENT  OF  EXTREME  LATERAL  SPREAD  CRITERIA 


HSie  sonic  boar,  distribution  cein  extend  over  large  lateral  distances 
if  the  ray  vhich  leaves  the  airplane  in  the  x-y  plane  (Pig.  II-4) 
reaches  the  ground.  HSie  conditions  inquired  to  produce  this  effect  can 
be  obtained  from  an  investigation  of  the  expression  fori„,  Eq.  (ll-le). 
For  any  ray  to  reach  the  ground  I  <0  throu^out  the  path  of  propa¬ 
gation.  ^ 

^us : 


*■  ■'  Mcos»-^ 


The  two  rays  which  leave  the  aiirlane  in  the  x-y  plane  are  the 
IP  =  0"  and  V  =  180*  rays.  Using  the  relationships  between  $  and  ^  , 
Eqs.(ll-2e)  and  (ll-?f)  for  t  =  0*  and  l80* : 


COStf  =  j^ 


—  X  M 


Thus,  the  ray  in  the  x-y  plane  will  reach  the  groiind  if 


I.  ^ 

- m 


Rearranging  this,  the  wind  required  to  produce  this  effect  is  given 


Ur  <  a.  —  a 


Ihe  quantity  Ur  can  be  'Titten  as 


US  =  U*  cose  +  V*sine 


Eq.  (Il-2c) 


which  for  the  %  »  0“  or  l8o®  ray  becoines: 


R~  mi 


U  :t  (M^  -  1)^^^  V  U,  d:(M^- V, 
M  M 


Uslag  this  In  the  above  Inequality  and  rearranging  slightly 

.  .  U.  :t  V.  ^  ^  .  U  d:  (M*- 1)^^*  V 

*•+ - M - >•+ - M - 


or 

A.  >  A  Eq.  (12) 

This  expression  shows  that  for  the  ray  In  the  X-y  plane  to  reach  the 
ground  the  value  of  the  expression  for  A  at  the  airplane  (A,)  must  be 
greater  than  at  any  point  between  it  and  the  ground. 


1*. 
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APPENDIX  Xm  REVIEW  OF  ROUTINE  CALCULATION  OF  SONIC  BOOM  IN 
STANDARD  AND  NON-STANDARD  ATMOSPHERES. 

In  order  to  aid  in  obtaining  rapid  estimates  of  the  sonic  boan  pro¬ 
duced  in  a  real  atmosjAiere  the  information  of  Sections  II. A. 4,  and  II.B.7 
is  repeated  here.  These  methods  must  be  used  with  extreme  care  for  making 
estimates  at  Mach  numbers  less  than  about  1.3>  criterion  developed 

in  Section  III  should  be  checked  for  fll^t  at  these  low  Mach  numbers  to 
determine  if  anamalous  propagation  may  occur.  The  conditions  \^lch  can 
caxise  complete  cut-off  are  given  in  Sections  III.A.l  and  III.B.l.  The 
conditions  \^lch  can  cause  focusing  at  the  ground  are  given  in  Sections 
III.A.2,  III.B.2,  and  III.B.3>  vhlle  the  conditions  \dilch  can  cause  ex¬ 
treme  lateral  spread  are  given  in  Section  III.B.U. 


A.  CALCULATIONS  IN  STANDARD  ATMOSPHERE  -  Once  the  airplane  geometiy 
Inputs,  KYg  ,0),  (see  Eq.  (l))  have  been  established  for  each  altitude 
and  Mach  number  of  Interest,  routine  calculations  of  sonic  boom  strength, 
lateral  dlstrlbtitlon,  and  lateral  extent  for  steady  level  fll^t  in  the 
U.S.  Standard  Atmosphere,  1962  may  be  obtained  in  the  following  manner: 

(a)  Compute  AP  nt.  track  (2)  and  Fig.  2  for  each 

Mach  nmber  and  altitude. 

(b)  Compute  lateral  dlstrlbtition  of  the  boom  strength  from  Eq.  (3) 
for  each  Mach  number  and  altitude. 

(c)  Obtain  the  location  of  lateral  cut-off  from  the  curves  in  Fig. 

L  for  each  Mach  number  and  altitude,  and  terminate  the  lateral  distribu¬ 
tion  of  boom  strength  at  this  point. 


B.  CALCULATIDNS  IN  NON-STANDARD  ATMDSPHERE  The  method  for  calcu¬ 

lating  sonic  boom  distribution  on  the  ground  was  outlined  above,  for  the 
U.S.  Standard  Atmosphere,  1962,  with  no  wind.  A  similar  procedure  would 
be  used  in  making  routine  estliaates  for  a  general  atmosphere  with  wind. 
The  following  method  should  be  used  with  extreme  care  for  Mach  numbers 
between  1.0  and  1.3,  especially  In  cases  when  wind  shears  are  to  be 
considered.  Criteria  for  meteorological  conditions  which  can  cause 
anomalies  in  the  overpressvire  and  distribution  of  the  sonic  boom  are 
developed  in  Section  III.  Ihese  should  be  checked  when  making  calcu¬ 
lations  in  the  above  Mach  number  range. 

Once  the  airplane  geometry  Influence,  I(Yo,ff),  (see  Eq.  (l))  has 
been  determined  for  each  Mach  number  and  altitude  of  Interest  the  cal¬ 
culation  of  sonic  boom  strength  and  distribution  may  proceed  as  follows: 

(a)  Calculate  AP  un^r /n,  ,„ck  ^9*  (5)  an^  Fig.  2  (or  Appendix 
Figs.  I-l,  1-2,  or  1-3)  for  ground  located  above  0  ft.  MSL  for  each  alti¬ 
tude  and  Mach  number. 
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•  The  variation  of  the  tempei'ature  correction  factor,  ,  may  he  esti- 
oated  free  Fig.  8,  with  the  value  at  Mach  1.2  estimated  from  Fig.  9* 

•  The  variation  of  the  correction  for  wind  shears  with  Mfich  number  may 
be  estimated  from  Fig.  12  with  the  value  at  Mach  1.3  taken  from 
Fig.  13. 

(b)  Compute  the  lateral  distribution  of  sonic  boom  strength  from 
Eq.  (3)  for  each  altitude  and  Mach  number.  (Caution  must  be  exercised 
when  estimating  the  lateral  distribution  for  low  Mach  numbers  and  high 
winds,  Fig.  15 ). 

(c)  Obtain  the  location  of  lateral  cut-off  from  Eq.  (6)  and  Eq. 

(7)  for  each  Mach  number  and  altitude,  and  terminate  the  lateral  distri¬ 
bution  at  that  point.  (Effect  of  moderate  winds  on  the  lateral  cut-off 
location  may  be  estimated  from  Fig.  I9. ) 

The  above  procedure  is  illustrated  schematically  in  Fig.  VIII-1. 


LATERAL  DISTANCE 


Fig.  XIU-l  Routine  Calculation  of  Sonic  Boom. 
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APPENDIX  H  DATA  TRANSFORMATION  FOR  TURBULENT  SCATTERING 


The  statistical  study  of  sonic  boon  overpressures  requires  an  ap> 
proprlate  tranefomatlon  of  the  test  data.  The  form  of  this  transforma¬ 
tion  may  be  determined  by  finding  the  probability  distribution  of  a  ran¬ 
domly  located  set  of  measurements  In  an  arbitrary  plane  from  the  mathe¬ 
matical  description  of  the  scattering  of  the  shock  wave  by  the  turbulence. 
Tbe  geometry  of  finding  the  data  treinsfoxmatlon  Is  shown  In  Fig.  IX-1. 


Fig.  IZ-1  C»om»try  for  Dotormlning  Probability  Distribution 


In  this  figure  a  -^urbulent  tenperature  or  wind  element  at  (A)  has 
scattered  the  energy  symmetrically  In  the  angle  <p  about  the  direction 
of  the  shock  propagation  IT.  A  number  of  Intensity  measurements  of  the 
scattered  field  are  made  by  microphones  placed  randomly  In  the  ( x,y ) 


or  ground  plane  and  It  Is  desired  to  treat  the  measured  data  statistically. 
The  actual,  form  of  the  probability  distribution  of  the  measurements  can  be 
fo\?nd  from  the  probability  theorem, 


P(x,y)=  yyf(x,y)dxdy  Eq.  (EC-l) 

R 

vhlch  states  that  the  Joint  probability  P(x,y)of  a  given  value  of  a 
variate  f(x,y)  occurring  In  the  region  R  In  the  (x,y)  plane  Is  egual  to 
the  Integral  of  f(x,y)over  the  region  R.  The  cross-section  «U) 
discussed  In  Section  III-C-3  gives  the  distribution  of  the  scattering 
Intensity  In  terms  of  the  solid  emgle.  When  this  dlstrlbutl^  Is  pro¬ 
jected  on  the  (x,y)  plane,  It  gives  the  desired  function  f(x  ,y) ,  thus, 
the  probability  dlstrlbxxtlon  becomes: 

p(4-]=  /o  (7)  •  dfA  Bq.  (IX-2) 

'  *o'  R 

vhere:  I  Is  the  observed  Intensity  of  the  scattered  wave 

Iq  Is  the  Incident  Intensity  of  the  original  wave 

dA  Is  the  vector  differential  area  In  the  (x,y)  plane 

0(/)ls  the  scattering  cr.  is-sectlon  discussed  In  Section  111-0-3 

If  the  definition  of  £q.  (lX-2)  Is  extended  to  a  number  of  scattering 
centers,  randomly  located  In  space  above  a  single  fixed  microphone,  the 
probability  distribution  becomes: 


P(jj=  E  f  a  Eq.  (lX-3) 

'O'  I  R 

I 


where  the  "i"  subscript  refers  to  one  measurement  taken  of  a  single 
scattering  event  of  an  ensemble  of  such  events. 

Figure  IX-1  shows  that  If  the  Incident  wave  vector  K,  Is  not  per¬ 
pendicular  to  the  (x,y)  plane  of  observation,  the  problem  Is  geometri¬ 
cally  complex,  requiring  the  Integration  of  elliptic  functions.  To 
Illustrate  how  the  form  of  the  probability  distribution  of  the  measure 
ments  may  be  obtained,  a  sljiq>llfled  case  of  scattering  by  a  turbulent 
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temperature  field  Is  developed  vihere  the  N-wave  is  props, gating  In  a  di¬ 
rection  perpendicular  to  the  plane  of  observation  (the  ground).  Sub¬ 
stitution  of  the  equation  for  the  txirbulent  temperatare  cross-section 
from  Section  UI-C 


a  ^d)  =  f  K*  cos2  q>  «  (IT) 


into  Eq.  (IX-3)  gives 


'*0/  ^  J=1  i  =  i  O  o  r,  z, 


Eq.  (n-U) 


where  z^  Is  the  height  of  the  ''l-th''Bcatterer  above  the  ground 

(r|,0i)are  polar  coordinates  in  the  horizontal,  ground  plane 
N  Is  the  total  number  of  observations 

Kj  Is  the  "j  -th"  wave  number,  which  contributes  In  an  unspeci¬ 
fied  manner  to  the  Incident  wave. 

After  performing  the  Integrations,  it  follows  that  the  measurements 
may  be  represented  by  the  form  of  the  probability  distribution: 


Eq.  (lX-5) 


The  summation  of  logarithms  In  general  cen  be  written  as: 


Eoj  1  nb^  =  ln[^n(b/‘j 


Eq.  (n-6) 
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vhlch  states  that  the  sum  of  a  set  of  logarithms  of  terms  in(ai,bj)i6 
equal  to  the  logarithm  of  the  product  of  terms  in  (bi)a  .  A  refomula- 
tlon  of  Bq,  (IX-5)  In  these  terms  leads  to  the  form  of  the  probability 
distribution  of  the  sonic  boom  measurements; 


K.  4>  (k)  I 


Eq.  (lX-7) 


Equation  IX-7  shows  that  the  form  of  the  probability  distribution 
is  logarithmic  and  consequently  a  logarithmic  transformation  of  the 
overpressure  data  would  be  appropriate  for  statistical  purposes  for 
this  case. 
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